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NONLINEAR  DIFFRACTION  OF  OCEAN  GRAVITY  HAVES 
ABSTRACT 

In  an  Irregular  sea,  waves  having  different  wave  numbers  Interact 
nonllnearly,  giving  rise  to  long  waves  at  the  difference  frequency  and 
wavenumber.  The  long  waves  are  associated  with  motions  that  have  large 
characteristic  time  and  space  scales. 

The  method  of  multiple  scales  in  time  and  in  space,  in  conjunction 
with  perturbation  expansions,  enables  us  to  separate  the  flow  into  com¬ 
ponents  for  the  general  case  when  there  are  wave-trains  propagating  in 
different  directions.  In  particular,  it  is  of  great  interest  to  study 
the  effect  of  modification  of  the  short  waves  by  diffraction,  refrac¬ 
tion,  reflection  and  radiation.  Using  the  method  of  matched  asymp¬ 
totics  we  determine  the  long-wave,  that  consists  of  forced  waves  tra¬ 
velling  with  the  short-wave  groups  and  of  an  additional  wave  that  pro¬ 
pagates  away  from  the  "zone  of  modification"  at  the  long-wave  velocity 
(gh)^/^.  The  resulting  theory  has  a  wide  range  of  applications.  We 
have  studied  the  following  problems: 

a)  Slow  Sway  of  a  Moored  Floating  Body  in  Hater  of  Finite  Depth 

Both  large  and  small  amplitude  cases  were  studied,  as  well  as  the 
effect  of  different  geometries  and  mooring  stiffness.  New 
analytical  results  for  a  block  with  small  clearance  underneath 
are  Included  and  the  existence  of  resonance  is  displayed. 
Periodic  and  different  transient  Inputs  were  studied. 
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b)  Wave  Trapping  on  a  Shelf 

Normal  and  obliquely  Incident  slowly  varying  waves  were  studied, 
retaining  short  wave  reflection.  For  the  obliquely  Incident  waves 
trapping  can  occur.  This  trapped  wave  cannot  be  excited  from  off 
the  shelf  by  any  linear  mechanism.  The  shelf  response  Is  computed 
and  singular  resonance  Is  found.  The  transient  response,  however. 
Is  finite.  It  too  Is  determined. 

c)  Excitation  of  Interfaclal  Waves  In  the  Lee  of  a  Breakwater 

The  method  of  multiple  scales  Is  used  to  determine  the  non-reso¬ 
nant  Interaction  of  surface  waves  with  Internal  waves  In  a  two 
layer  fluid.  We  study  the  effect  of  diffraction  by  a  long  break¬ 
water.  While  the  short  waves  are  decaying  behind  the  breakwater, 
free  long  Internal  waves  are  propagating  directly  Into  that 
shadow  zone. 

The  present  method  of  analysis,  and  In  particular  the  multiple 
scales  expansion,  proves  to  be  a  useful  tool  In  studying  modulation  and 
nonlinearity  In  several  aspects  at  wave  propagation,  and  can  be 
extended  to  the  study  of  a  heretofore  unexplored  aspect  of  harbour 
resonance.  The  Ideas  will  be  sketched  In  the  Epilogue. 


Thesis  Supervisor:  Chlang  C.  Mel 

Title:  Professor  of  Civil  Engineering 
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INTRODDCTION 


1.1  Survey  of  Literature 

During  the  last  fifteen  years  there  has  been  an  Increasing  inter¬ 
est  in  the  low  frequency  effects  of  irregular  waves  interacting  with 
floating  bodies.  The  low  frequency  force  and  motion  can  be  suffic¬ 
iently  large  and  are  Important  to  the  design  of  moored  and  dynamically 
positioned  offshore  structures.  They  can  also  affect  the  behaviour  of 
ships  as  well  as  that  of  Icebergs. 

Moored  offshore  platforms  are  subject  to  seas  with  narrow-banded 
spectra.  The  moored  vessel  may  have  natural  frequencies  of  horizontal 
plane  motions  (sway,  surge,  yaw)  at  the  order  of  0.01  Hz,  which  are 
much  lower  than  the  dominant  frequencies  of  the  Incident  sea,  typically 
0.1  Hz.  It  has  been  observed  that  the  horizontal  response  of  a  moored 
vessel  is  Indeed  the  largest  at  the  low  frequency  range  around  its 
natural  period  (Verhagen  &  Van  Sluijs  (1970),  Hsu  &  Blenkarn  (1970), 
Remery  &  Herman  (1971);  see  Fig.  1.1).  At  these  frequencies  the  prime 
source  of  forcing  is  associated  with  the  modulation  periods  of  short¬ 
wave  groups.  Nonlinear  interaction  of  the  high  frequency  waves  forces 
waves  at  the  difference  (beat)  frequencies  and  wavenumbers. 

Several  authors  have  approached  the  long  period  motion  within  the 
potential  theory.  A  comprehensive  review  of  this  topic  has  been  given 
by  Ogllvie  (1983),  who  also  cites  a  few  references  that  treat  viscous 
effects.  In  a  regular  wave  train  the  steady  drift  force,  which  is 
second  order  in  wave  slope,  can  be  computed  from  the  first  order 
(linear)  solution.  In  Irregular  waves,  Newman  (1974)  has  found  that 
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y  Response  to  a  Narrow  Banded  Wave 


the  slow  drift  force  can  be  written  as  a  quadratic  "transfer"  function 
of  the  wave  components.  The  coefficients  of  this  function  can  be  expan¬ 
ded  as  functions  of  the  difference  frequencies.  He  suggests  that  for 
small  frequency  differences  the  coefficients  can  be  approximated  by 
their  values  at  zero  difference  -  thus  the  slowly  varying  drift  force  is 
almost  as  simple  as  the  steady  drift  force. 

Newman's  result  has  been  supported  by  Martinsen  (1983)  following  an 
Idea  suggested  by  Hsu  &  Blenkarn  (1970).  In  the  time  domain,  the  slowly 
varying  wave  Is  viewed  as  a  regular  wave  train  for  any  given  Instant, 
with  amplitude  and  frequency  of  the  Irregular  wave  at  that  Instant  and 
place.  The  slowly  varying  force  Is  then  assumed  to  be  given  by  the 
steady  drift  force  that  would  correspond  to  a  regular  wave  train  that 
has  those  Instantanous  and  local  characteristics.  Martinsen  points  out 
In  addition  that  for  certain  geometries,  small  frequency  changes  may 
result  In  large  changes  In  the  steady  drift  force;  this  was  not  Included 
In  the  analysis  of  Roberts  (1981)  who  used  the  same  approach  to  study 
drift  motion  In  random  seas  with  nonlinear  restoring  forces. 

In  computing  the  slow  drift  force,  the  pressure  Is  Integrated  over 
the  wetted  surface  of  the  body  and  the  low  frequency  terms  are  collect¬ 
ed.  Nearly  all  of  these  terms  have  corresponding  terms  In  the  expres¬ 
sion  for  the  steady  drift  force  due  to  a  regular  wave  train,  with  the 
exception  of  a  term  that,  includes  a  slowly  varying  potential.  This  pot¬ 
ential  has  not  been  Included  by  Newman,  Martinsen  or  Roberts.  Bowers 
(1976)  computed  the  low  frequency  potential  forced  by  two  regular  wave 
trains  without  diffraction  effects.  A  similar  approach  was  followed  by 
Pinkster  (1980).  Vlhen  the  water  depth  Is  small  compared  with  the  wave- 


9 


length  of  the  slow  wave,  this  potential  becomes  large  and  needs  to  be 
considered . 

Trlantafyllou  (1982)  has  observed  that  for  finite  depth  the  slow 
potential  is  of  first  order.  He  used  a  multiple  time  expansion  to 
study  large  amplitude  (0(1))  motion;  this  technique  was  also  employed 
by  Molin  &  Bureau  (1980).  Using  more  traditional  arguments,  Faltinsen 
&  Loken  (1979)  have  calculated  the  contribution  of  the  slowly  varying 
potential  by  introducing  a  related  boundary  value  problem  and  using 
Green’s  Theorem  for  two  dimensional  bodies  in  water  of  infinite  depth. 
They  have  concluded  from  computed  results  that  the  contribution  of  the 
term  in  question  is  small  when  the  wave  spectrum  is  narrow  banded. 

The  related  problem  of  diffraction  by  fixed  coastal  structures  is 
also  of  engineering  importance.  The  linear  problem  has  been  exten¬ 
sively  studied  in  the  literature  (see  Mei  (1983)).  In  particular,  the 
propagation  of  waves  which  are  short  compared  with  the  structure's  dim¬ 
ensions  can  be  approximated  by  the  theory  of  geometric  optics  and  the 
parabolic  approximation.  Not  much  theory  is  available  on  the  nonlinear 
aspects . 

Interaction  of  the  long  waves  that  accompany  groups  of  short  waves 
can  also  occur  when  the  wave  motion  is  modified  by  topography  (refrac¬ 
tion)  or  by  coastal  structures  such  as  islands  or  harbours  (diffrac¬ 
tion).  Mei  &  Benrooussa  (1984)  have  discussed  the  refraction  of  slowly 
varying  waves  due  to  slowly  varying  depth  where  the  length  scales  of 
the  variations  are  comparable.  They  have  found  a  nonlinear  mechanism 
for  producing  free  long  waves  which  do  not  propagate  in  the  same  direc¬ 
tion  as  the  wave  group.  In  particular  long  waves  can  be  trapped  over  a 
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ridge.  Foda  and  Mel  (1981)  have  studied  such  trapping  using  high  order 
theory.  Symonds  and  Bowen  (1984)  have  studied  trapping  of  long-waves 
by  breaking  short-waves. 

Bowers  (1977)  has  studied  a  problem  that  Is  essentially  one  dimen¬ 
sional,  where  a  narrow  rectangular  wave  flume  resonates  a  very  narrow 
rectangular  basin.  He  has  shown  that  the  free  long  waves  Induced  by 
Incident  short  wave  groups  can  resonate  such  a  harbour.  The  resonant 
periods  are  much  longer  than  those  of  the  Incident  short  waves,  but  may 
be  close  to  their  modulation  period. 

1.2  Scope  of  Thesis 

In  this  thesis,  both  floating  bodies  and  large  fixed  structures 
will  be  treated. 

It  Is  well  known  that  wave  motion  with  long  periods  Is  associated 
with  long  length  scales.  The  modulation  length  scale  and  the  long 
waves  associated  with  it  are  expected  to  be  much  larger  than  the  short 
wave-length.  This  suggests  the  introduction  of  multiple  scales  In 
space  as  well  as  In  time.  In  Chapter  2  a  general  theory  Is  presented 
for  the  interaction  of  short  waves  with  long  waves  in  water  of  finite 
depth.  The  motion  of  a  floating  body  that  performs  slow  drift  motion 
in  water  of  finite  depth  Is  our  first  application.  Both  large  and 
small  drift  displacements  are  studied  in  Chapter  3.  The  effects  of  the 
determining  factors,  the  mooring  stiffness,  and  the  body  geometry  are 
examined.  Both  two  and  three  dimensional  geometries  are  discussed  In 
general  terms.  A  special  two-dimensional  configuration  which  exhibits 
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damped  resonance  Is  that  of  a  horizontal  rectangular  block  with  a 
narrow  gap  underneath  it.  This  geometry  Is  studied  analytically,  and 
compared  to  a  block  with  no  gap.  Solutions  for  transient  evolution 
that  feature  initial  growth,  approach  to  quasi-steady  state  and  decay 
are  also  given. 

Slowly  varying  wavetralns  may  interact  with  the  sea  bottom  topo¬ 
graphy  as  well  as  with  coastal  structures  such  as  Islands  and  harbours. 
The  characteristic  length  of  these  structures  is  much  larger  than  that 
of  the  floating  bodies.  To  illustrate  the  physics,  we  study  a  few 
simple  geometries.  In  Chapter  4  the  interaction  of  slowly  varying 
waves  with  a  shelf/ridge  Is  studied.  For  oblique  Incidence  trapping  of 
waves  and  resonance  are  found.  The  transient  response  of  the  shelf  Is 
also  studied,  displaying  a  'wake*  effect  after  the  Incident  waves  pass 
through. 

In  Chapter  5  we  study  the  diffraction  of  slowly  varying  surface 
waves  by  a  long  breakwater  In  a  two  layer  fluid.  In  particular,  long 
Interfaclal  waves  are  forced  by  the  wave  groups.  We  show  that  free 
long  Internal  waves  are  emitted  from  the  shadow  boundary  Into  the  lee 
of  the  breakwater  and  compute  their  magnitude.  In  Chapter  6  the  work 
Is  concluded.  Limitations  and  possible  extensions.  Including  harbour 
resonance,  are  discussed. 
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THE  INTERACTION  OF  SHORT  PERIOD  HAVES  HITH  UNIG  PERIOD  HAVES 


2.1  Introduction  and  Definitions 

It  Is  well  known  that  progressive  waves  Induce  a  depression  In  the  mean 
water  surface  level,  which  Includes  the  "set  down",  as  well  as  a  drift 
current,  which  includes  "Stokes'  drift".  These  quantities  are  averages  over 
the  short  period  of  the  surface  wave.  They  are  second  order,  proportional  to 
the  square  of  the  wave  steepness,  Longuet-Higglns  and  Stewart  (1962,1964) 
have  obtained  expressions  that  relate  these  mean  quantities  to  the  'radiation 
stress'  which  stands  for  the  momentum  flux  due  to  the  wave  propagation. 

When  the  amplitude  of  the  short  waves  varies  slowly  in  time  and  space, 
the  'set  down'  and  the  drift  vary  correspondingly,  at  the  same  spatial  and 
temporal  rates.  They  are  proportional  to  the  square  of  the  amplitude  of  the 
primary  waves.  This  slow  variation  has  the  nature  of  a  wave,  and  will  be 
referred  to  as  a  long,  or  low  frequency,  wave.  By  carrying  out  a 
perturbation  expansion  of  the  wave  potential  and  using  the  technique  of 
multiple  scales,  we  study  the  interaction  between  a  rigid  structure,  groups 
of  slowly  varying  waves  that  propagate  in  different  directions,  and  the  long 
period  waves  that  are  Induced  by  them.  Specific  examples  are:  floating 
bodies,  depth  variations,  and  breakwaters. 

Under  the  usual  assumptions  of  potential  theory  (ideal  Incompressible 
fluid,  Irrotatlonallty,  no  separation  etc,)  we  may  write  Laplace  Equation  for 
the  velocity  potential  ♦(x,  y,  z,  t): 

A*  ■  0  in  the  fluid  (2,1,1) 

where  (x,y,z)  are  Cartesian  coordinates,  with  the  positive  z  axis  pointing 
vertically  upwards,  and  t  denotes  time.  Writing  g  for  gravity  acceleration. 
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P  for  pressure  and  p  for  the  fluid  density,  we  may  write  Bernoulli's  Equa¬ 
tion: 

1  o 

-■^“gz  +  *t'^2  h*|  (2.1.2) 

We  let  the  pressure  at  the  free  surface  be  zero.  Assuming  small  wave  steep¬ 
ness,  and  Taylor  expanding  the  free  surface  boundary  condition  for  ♦  around 
the  rest  position  of  the  free  surface,  we  get: 

*tt  8‘z  ■  [■  i  *  i  ‘t  ‘ztlt  ■  V  ‘t> 

(2.1.3) 

+  0(7j^  *^,  («^)^)  (z  =  0) 

(cf.  Benney  (1962)). 

At  the  rigid  bottom,  the  kinematic  boundary  condition  is: 

*z  “  ^h  ^  (2.1.4) 

where  Vh  =  (3/3x»  9/3y)  is  the  horizontal  gradient,  and  h  *  h{x,  y)  is  the 
bottom  depth. 

The  water  is  assumed  to  be  of  intermediate  depth,  comparable  to  the  wave 
length: 

kh  =  0(1)  (2.1.5) 

We  shall  study  only  bathymetrles  that  are  piecewise  horizontal  (Vh  is  zero) 

for  which  (2.1.4)  reduces  to: 

♦  “  0  (z  -  -h)  (2.1.6) 

z 

There  are  two  small  parameters  that  are  associated  with  the  peturbation 
expansion  for  slowly  varying  small  amplitude  waves.  The  first  one  is  the 
wave  steepness: 

e  -  0(k  eA)  «  1  (2.1.7) 

where  k  is  the  central  wavenumber  and  eA  is  the  free  surface  amplitude  of  the 
short  wave. 


14 


The  second  parameter  Is  the  modulation  ratio: 

W  -  0(pft/a))  «  1  (2.1.8) 

where  ui  “  0(g/h)^/^  Is  the  central  frequency  of  the  short  wave,  and  yfi  Is 
the  frequency  of  modulation  of  the  short  wave,  or  -  equivalently.  Its 
frequency  bandwidth.  For  simplicity,  we  shall  choose  e  equal  to  y,  so  as  to 
render  the  effects  of  dispersion  and  nonlinearity  comparable. 

The  vastly  different  spatial  and  temporal  scales  of  the  short  waves  and 
their  envelopes,  suggest  the  Introduction  of  multiple  scales.  This  will 
enable  us  to  distinguish  between  quantities  of  different  magnitudes  and 
obtain  a  consistent  analysis  of  the  different  orders. 


2.2  The  Multiple-Scales  Analysis 

Since  the  wave  envelopes  have  typical  temporal  scales  of  O(eto)”^  and 
spatial  scales  of  0(ek)“^,  we  define  the  stretched  coordinates: 


(Xj,  yj,  tp  -  e(x,  y,  t)  (2.2.1) 

Since  the  fluid  depth,  h.  Is  assumed  to  be  much  smaller  than  the  long  scale, 
there  Is  no  need  to  stretch  the  vertical  coordinate. 

The  potential,  $,  and  the  free  surface  displacement,  q,  are  expanded  In 
powers  of  the  small  parameter  e: 

2 

(♦,  ?)  “  g(4j,  Cj)  E  (^2*  ^2^  ***  (2.2.2) 


where  *11 


and  functions  of  the  fast  and  the  slow  coordinates: 

(X,  y,  X,  t,  X,,  y,,  tj).  -  c^(x,  y,  t,  Xj .  y,,  tj)}  - 


0(1) 

(2.2.3) 


All  quantities  In  this  work  are  0(1)  with  powers  of  c  signifying  the  order. 
The  partial  derivatives  become: 
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(2.2.4) 


a/9t  (e"  ♦)-£”{♦  +  e  ♦  ] 


Expanding  (2.1.1)  In  this  fashion,  we  obtain,  after  separating  the 


orders: 


-  0 


(2.2.5a) 


^*2  “  ■^♦ixxj  ■  ^♦lyyj 


(2.2.5b) 


In  the  same  way,  (2.1.3)  gives  at  (z  *  0): 

+  gi,  **  0 
Itt  ®  Iz 


(2.2.6a) 


»2tt  *  «*2z  ■  -  i  Vl>  <2.2.1 


and  from  (2.1.6),  at  (z  ■  -h): 


*lz  *  “ 


(2.2.7a) 


‘2z  ■  ° 


(2.2.7b) 


The  dependence  on  fast  time,  t,  is  in  the  form  where  m  is  an 

integer.  We  further  expand  terms  of  these  time  harmonics. 


defined  by: 


<*n>  'a'  ■  I  <*nm-  'am>' 


-im<i)t 


(2.2.8) 


m*  -n 


where  are  complex,  containing  the  phase  information.  We  require 
♦nm  *  ♦n,  -m  (where  (  )*  is  the  complex  conjugate)  so  that  are 
real. 


In  the  following  sections,  we  shall  solve  the  perturbation  equations  at 
the  different  orders  and  harmonics  that  are  required  to  determine  the  leading 
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order  slow  motion  and  long  waves.  The  first  order  first  harmonic  potential 
represents  the  short  waves.  Because  of  the  use  of  multiple  scales,  the 
long  wave  potential,  too,  is  first  order.  It  Is  the  zeroth  harmonic,  ♦jg* 

By  focusing  on  the  long  spatial  scales  we  shall  obtain  a  governing  equation 
for  t]ig.  In  terms  of  the  solution  for  4^^.  In  this  chapter,  only  the  general 
form  of  the  solution  will  be  studied.  The  specifics  of  the  near  field  solu¬ 
tion  and  matching  to  the  far  field  are  left  to  the  following  chapters.  We 
begin  by  studying  the  short  waves  potential, 

2.3  The  Short  Haves  Away  froa  the  Body 


Extracting  the  first  order  first  harmonic  equations,  for  the  primary 
short  wave  potential,  from  (2.2.5a,  2.2.6a,  2.2.7a),  we  obtain: 

(2.3.1) 


“a  ■  “ 


»uz  -  ‘>‘a  ■  “ 


(-h  <  z  <  0) 
(z  -  0) 


(2.3.2) 


where  a  *  — 
g 


(z  -  -h) 


(2.3.3) 


These  equations  are  formally  the  same  as  the  linearized  equations  for 
purely  periodic  waves  with  frequency 

He  focus  our  attention  on  a  form  of  the  solution  that  Is  general  enough 
for  the  classes  of  problems  studied  in  the  present  work.  This  form  assumes 
periodicity  of  the  short  wave  In  the  y  direction: 

*11  “  ^0^^^  ^1'  +  ®0^*1»  *^1^® 


(2.3.4) 


+  I  f^(z)  ls^(xj,  yj,  tj)  e 
n“l 


a  x-Hyy  _  -a  x+lyy 

"  +  Sq(xj,  yj,  tj)e  "  1 
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w* 


;  n  -  0,  1,  2 


are  the  eigenfunctions: 


(2.3.5a) 


(2.3.5b) 


(2.3.6) 


We  may  also  write;  a  •  k  sin  ©,  y  »  k  cos  9,  where  9  is  the  direction  of  pro¬ 
pagation  of  the  propagating  modes.  In  any  given  region  either  or  {s-} 

n  *  1,2...  may  be  present,  and  their  contribution  Is  decaying  exponentially. 
These  modes  are  called  the  evanescent  modes.  When  the  water  depth,  h,  takes 
on  different  but  constant  values,  the  wavenumbers  end  eigenfunctions  must 
change  according  to  the  dispersion  relation  (2.3.5b). 

The  propagation  of  the  wave  envelopes  for  the  propagating  modes,  Sq  and 
Sq,  can  be  found  from  solvability  conditions  for  the  second  order  first  har¬ 
monic  potential  ♦21*  view  of  (2.3.4),  the  governing  equations  for  the 
propagating  modes,  e.g.,  the  component  associated  with  Sq,  can  be  separated 
from  (2.2.5b,  2.2.6b,  2.2.7b),  to  get: 

A*2i  -  ®0x  ®0y  ^  (-h  <  z  <  0)  (2.3.7) 
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(z  -  0) 


(2.3.8) 


21(1)  +  , 

*21  ■  ®*21  “  g  ®0t,  *0 

*21  *  0 
z 

For  solvability  of  ♦2i»  w®  must  have: 

®0  “  ®o  ~  *1  ~ 


where 


M 


1“  ;  V  - 
3(1)  3(1) 


(cf.  Mel  (1983,  p.  52)). 

The  propagation  of  Sq  Is  found  In  the  same  way.  It  Is: 


(2.3.9) 


(2.3.10) 


(2.3.11) 


Sq  -  Sq  (tj+  yXj  -  vyp  (2.3.16) 

We  stress  that  the  details  of  Sq  are  still  unknown,  only  the  dependence 

on  x^,  y^,  and  t^  Is  fixed.  Further  Information  must  await  the  study  of  the 
neighborhood  of  the  body  where  evanscent  modes  are  Important. 


2.4  The  IfOng  Waves  Away  from  the  Body 


We  now  turn  to  the  long  waves.  They  are  associated  with  the  zeroth 
harmonic  of  the  potential,  and  of  the  surface  elevation.  Due  to  the  Intro¬ 
duction  of  stretched  coordinates.  It  is  the  first  order  slow  potential  that 
gives  rise  to  second  order  free  surface  displacement 


ii_  * 

■  g  3t  ^*10 


1  2 

i  "  *ioti 


O(e^) 


(2.4.1) 


We  shall  proceed  to  find  the  governing  equation  for  *^g. 

The  Interaction  of  short  period  waves  and  long  period  waves  Is  different 
in  open  water,  over  large  horizontal  scales,  from  that  in  regions  of  wave 
generation,  such  as  the  nelgborhood  of  a  floating  body,  a  depth 
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discontinuity,  or  other  structures,  where  evanescent  short  waves  exist.  Each 
region  will  be  studied  separately  and  the  resulting  solutions  will  be  matched 
asymptotically  to  determine  their  values.  The  region  where  the  evanescent 
modes  are  Important  -  that  Is  -  within  a  few  short  waves  from  the  wave 
generation  area  -  Is  referred  to  as  the  'near  field'.  We  denote  the 
potential  In  the  near  field  by 

1(1  -  ♦  kx  -  0(1)  kxj  «  1  (2.4.2) 

The  region  In  which  only  propagating  modes  need  to  be  considered,  a  few  wave 
groups  away  from  the  wave  generation  area.  Is  called  'the  far  field'  and  the 
potential  there  Is  denoted  by 

*  -  ♦  kx  »  I  kxj  -  0(1)  (2.4.3) 

see  Figure  2.1 . 

With  this  notatloil  we  shall  study  separately  the  specific  characteris¬ 
tics  of  the  potential  In  each  region. 

The  short  scale  variation  of  ♦jq  Is  found  by  collecting  the  zeroth  har¬ 
monic  from  (2.2.5a),  (2.2.6a)  and  (2.2.7a): 


> 

o 

1 

o 

(-h  <  z  <  0) 

(2.4.4) 

A  «  0 

*10z 

(z  -  0,  -h) 

(2.4.5) 

The  eigenfunctions  for  these  equations  are: 

g^(x,  y)  cos  n  -  0,  1,...  (2.4.6a) 

where 

(Ajj  -  n^x^/h^)  ••  0  (2.4.6b) 

For  all  n  #  0,  the  solutions  decays  exponentlaly  away  from  the  near 
field.  We  conclude  that  the  far  field  component  of  Is  Independent 
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Figure  2.1  Definition  Sketch 
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of  the  short  scale  (x,  y,  z): 

“  ^10^*1*  ^1*  ^1^  (2»A,7) 


Let  us  now  seek  the  governing  equations  for  mainly  from  (2.1.3). 

To  second  order  In  e,  ^iQz  zero  at  (z=0).  However,  by  Integrating 
(2.1.1)  from  -h  to  0  with  respect  to  z,  and  using  (2.1.6),  we  find: 


0 


0 


♦  n  I  dz  =  -  f  V  t  dz 
z  |z=0  J.  ^zz  J  h  * 

~h 


also 


^^h  *^30 


-h 


*10  ^^h‘  ^1  ^  ^1*  ^h^*20  ^  \  *30 


(2.4.8) 

(2.4.9) 


where  =  (3/3Xj,  3/3yi) 

The  first  term  on  the  right  hand  side  of  (2.4.9)  is  only  a  function  of 
the  long  space  scales  (xj,  y^)  while  the  other  two  terms  are  functions  of 
short  (x,  y)  as  well.  The  third  order  zeroth  harmonic  of  (2.1.3)  is 


^8  ♦z^30 


^lOtit^ 


(z“0) 


We  substitute  (2.4.9)  into  (2.4.8)  and  the  resulting  equation  into  (2.4.10). 
The  long  spatial  scale  component  of  (2.4.10)  is  thus  given  by: 


•^lOt^tj  ^1  *10 


(2.4.11) 


On  the  right  hand  side  of  (2.1.3),  the  long  scale  component  of  {0[(^^)5j], 
0l(^^)y]}  is  of  order  0(6**)  while  the  zeroth  harmonic  of  {0[(^3)^]j 
0(6**).  The  third  order  zeroth  harmonic  of  the  other  terms  on  the  right  hand 
side  of  (2.1.3)  is 
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^  ♦  (O  tn  *  *»t, 

-  *11  ♦u>  ’h*’! 


(2.4.12) 

'^*11  '*’ 


We  note  that  In  the  far  field.  In  the  absence  of  evanescent  modes  In  (2.3.4), 
the  terms  of  (2.4.12)  have  components  of  the  form: 

(s-  s;*e-2^“*  *)  , 


+  2 

«0 

Clearly,  only  the  last  two  terms  contribute  to  the  long  spatial  scale  compo¬ 
nent.  This  means  that  (jijg  is  Induced  only  by  the  self  Interaction  of  each 
propagating  short  wave  with  itself.  Cross  Interactions  between  the  propagat¬ 
ing  modes,  or  with  the  evanescent  modes,  give  rise  to  higher  order,  short 
scale  waves.  In  the  second  square  brackets  of  (2.4.12)  only  the  first  term 
is  purely  long  scale  in  space,  and  the  last  term  is  discarded  since  it 
depends  on  short  space  scale. 

When  only  the  self  interaction  terms  are  taken  on  the  right  hand  side 
of  (2.1.3)  the  result  is: 


^  2  ^2.^.  r,  2  ,  2.  3 

♦lOtjtj  ■  **•  'l  *10  ■ 


*  -  2i»a  Tr- -  2»ir  |»J|' 


3x 


1 


ay] 

3 


+  fo(0)  [(a  -  k  )  +  2a,a  -  2a,Y  |s-| 


(2.4.13) 

-|2 


which  is  the  differential  equation  governing  the  far  field  slow  potential 
♦io(*i»  yp 
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where  the  overbar  stands  for  the  long  scale  component,  the  terms  make  up 
(2.4.11)  while  the  terms  due  to  the  short  waves  make  up  (2.4.12).  Hence,  the 
wave  equation  for  (2«4.13)  can  be  retrieved.  Finally,  a  third  and  more 

formal  way  to  obtain  (2.4.13)  is  by  studying  solvability  conditions  for  ♦30* 
this  is  lengthy  and  will  be  omitted. 

The  form  of  (2.4.13)  can  be  further  simplified.  By  making  use  of  the 
relations  for  the  propagation  of  the  wave  envelopes  (2.3.10)  and  (2.3.11), 
(2.4.13)  becomes: 

'•’lOtjt  "  ®^^1  '•*10  “  ^0^®^  ~  [jsj  ^  +  |sq|^1^.^  (2.4.18) 


24 


which  is  hyperbolic.  Finally  we  write: 


ay  +  ykv 


g  g 


(2.4.19) 


(2.4.18)  then  becomes: 


hot  t  -  eK  ho  -  -  2“  ^ 


'lOtjt, 


(2.4.20) 


This  is  a  generalization  of  the  theory  of  Longuet-Higgins  and  Stewart 
(1962,  1964)  for  a  single  wavetrain. 

The  wave  Equation  (2.4.13)  or  (2,4.20)  governs  the  long  wave  in  the  far 
field.  The  specific  solutions  will  depend  on  the  particular  problem  at  hand 
and  will  be  determined  by  matching  to  the  ne?  ,x>sid.  In  the  following  chap¬ 
ters  we  shall  study  a  few  examples  ol  the  application  of  the  present  theory 
to  problems  that  Involve  wave  propagating  in  different  direction  at  the  same 
location. 
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3.  SLOW  SWAT  OF  A  FLOATING  BODY 

3.1  Introduction  and  General  Considerations 

We  first  deal  with  the  slow  motion  of  moored  floating  bodies  in  slowly 
varying  waves.  The  linearized  theory  does  not  account  for  low  frequency 
responses  due  to  narrow  banded  incident  swells.  It  has  been  observed,  how¬ 
ever,  that  such  responses  exist,  and  can  be  very  large.  It  is  also  known 
that  the  resonant  frequency  of  a  mooring  system  usually  lies  below  the  swell 
frequencies,  while  it  can  fall  Inside  the  range  of  the  lower,  modulational 
frequencies . 

In  general,  a  floating  body  has  six  degrees  of  freedom.  Among  them, 
three  modes  of  motion  Induce  buoyancy  restoring  forces  and  moments.  These 
are:  heave,  pitch  and  roll.  To  leading  order,  the  buoyancy  forces  and 
moments  are  proportional  to  the  displacements.  For  example,  a  heave  motion 
with  amplitude  H  will  Induce  a  buoyancy  force  with  amplitude  -pgHA^,  where 
denotes  the  water-plane  area.  H  can  be  large  only  for  a  bottle  shaped 
body,  on  which  we  shall  comment  later. 

For  a  body  that  is  not  bottle  shaped,  the  only  Interesting  slow  motions 
are  horizontal:  surge,  sway  and  yaw.  We  shall  examine  a  sway  problem  to 
Illustrate  our  approach.  The  geometry  chosen  is  that  of  a  rectangular  hori¬ 
zontal  cylinder  in  beam  seas.  The  cylinder  axis  lies  on  the  y  axis  and  the 
waves  are  incident  from  (x  -»•  -<»)  (Figure  3.1).  The  problem  is  then  two 
dimensional,  with 

3/ay“0  (3.1.1) 

for  all  quantities. 
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We  assume  that  the  body  performs  sway  motion  only,  to  the  leading  order, 
at  both  fast  and  slow  time  scales.  This  assumption  simplifies  the  presenta¬ 
tion  and  will  be  shown  to  be  valid  for  some  cases.  Computed  results  for  the 
general  motion  will  be  presented  in  the  other  cases.  The  mooring  Is  modelled 
by  a  linear  spring. 

The  wave  field  is  coupled  with  the  motion  of  the  body  through  boundary 
conditions  on  the  body  surface.  The  Incoming  waves  are  diffracted  and 

reflected  by  the  presence  of  the  body.  Waves  are  further  radiated  by  the 

body  motion  caused  by  waves.  The  effect  of  sway  Is  expressed  by  the  kine¬ 
matic  boundary  condition.  Denoting  the  sway  displacement  by  X: 

*x  “  ^t  *  {x  -  X(t)  ±  B;  -D  _<  z  £  5}  (3.1.2) 

$  =  0  on  b  =  (z  =  -D;  -B  <  x  -  X  <  B I  (3.1.3) 

S-  are  the  sides  of  the  body,  which  has  breadth  2B  and  draught  D  (see 
Figure  3.1).  B  and  D  are  assumed  to  be  comparable  to  the  depth  h  which  is  of 
the  order  of  the  wavelength  2Tr/k.  The  boundaries  of  the  body  at  rest  are 
given  by 

Sq  =•  {x  =  ±B;  -0  <  z  <  0},  Sq  =  {z  =  -D;  -B  £  x  _<  B} 

The  effect  of  the  waves  on  the  body  motion  is  given  by  the  dynamic 
boundary  condition; 

MX  +  Kn*  “  /  P  dz  -  /  P  dz  (3.1. A) 

S  S 

where  M  is  the  mass  of  the  body  and  is  the  elastic  constant  of  the  mooring 
system.  The  right  hand  side  is  the  hydrodynamic  force  where  the  pressure  P 
is  given  by  (2.1.2).  We  shall  see  that  the  driving  force  for  the  slow  motion 
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(zeroth  harmonic)  Is  O(e^).  The  magnitude  of  X  depends  on  the  mooring  stiff¬ 
ness:  If  Kq  Is  0(e),  X  must  be  0(e)  (small  displacement)  in  order  to  balance 
the  hydrodynamic  force.  This  case  Is  studied  In  Section  3.3.  The  mass  of 
the  body  is  0(1)  so  that  the  Inertia  of  the  body  for  the  slow  motion  Is  O(e^) 
since, 


MX 


tt 


2 

e 


MX. 


(3.1.5) 


It  Is  negligible  if  X  Is  0(e)  (small  displacement). 

On  the  other  hand.  If  Kq  Is  O(e^),  X  must  be  0(1)  In  order  to  balance 
the  slow  drift  force.  The  body  Inertia  Is  then  O(e^).  This  case  will  be 
studied  In  Section  3.4. 

As  a  special  but  Interesting  case,  we  shall  see  that  large  blockage  by 
the  cylinder  may  give  rise  to  large  added  mass  and  damping  coefficients, 
rendering  X  at  0(e)  regardless  of  the  mooring. 

The  problem  will  be  studied  separately  in  the  near  field  ,  within  a  few 
short  waves  from  the  body,  and  In  the  far  field,  a  few  wave  groups  away  from 
the  body  on  either  side.  We  examine  the  form  of  i  In  the  far  field  where  it 
Is  denoted  (^.  Because  of  radiation  conditions,  the  only  waves  in  addition 
to  the  Incident  waves  must  be  propagating  away  from  the  body.  Since  the 
equations  for  are  linear,  (2.3.10)  and  (2.3.16)  for  Sg  and  Sg  may  be 
written  as: 


sit  =  Ag(tj  -  Xj/Cg)  ;  Sn  -  RA^Ct,  +  tt^)  (x,  <  0) 


’0 


O'  1 


g 


(3.1.6) 


So  -  TAg(tj  -  Xj/Cg)  ;  Sg  -  0 


(xj  >  0) 


where  R  and  T  are  complex  reflection  and  transmission  coefficients. 
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respectively.  We  shall  see  later  that  they  Include  both  diffraction  and 
radiation  effects.  The  equation  for  <^jq,  (2.4.20),  becomes: 


*^10tjtj  "  *10xjXj 

f^CO)  [(o^  -  k^)  Cg  +  2u)k] 


to  the  left  and  the  right  of  the  body,  respectively. 

The  forcing  terms  produce  waves  which  propagate  at  velocities  Cg  and 
-Cg.  There  are  also  free  waves  which  propagate  away  from  the  origin  at 
velocities  (gh)^/^  and  -(gh)^/^.  These  are  solutions  to  the  correspond¬ 
ing  homogeneous  equation.  Formally,  the  solution  can  be  written  as: 


^in  *  { 


^10  ^*1  ~  ^g  ^  *^10  ^*1  *  ^g  ^  *^10  ^*1  ^  *"1^  0) 


<t*Io  ^*1  ~  ^lO  ^*1  ■ 


(x^  >  0) 

(3.1.8) 


The  values  of  (a  =  I,  R,  or  T)  at  (xj  =  0)  can  be  obtained  from  (3.1.7) 

and  (3.1.8) 


(a)  _  (a)u2|  ^0^^^ 


2  2.  2(0  k 


[x.  ]  -  [1] 


gh/Cg  -  1 


'J  r/i  *  1  r  vx;.  r  ^ 

■  *  “I  1  ■  ’  ’  57“’)  'Vi’  ■  ‘  ’ 

^20-  will  be  determined  through  matching  to  the  near  field  i|)jo*  order 
to  perform  the  matching,  we  shall  need  the  inner  expansion  of  (3.1.8)  as 
(xj  -»■  0)  for  (xj  >  0)  and  (x^  <  0),  which  is 


♦10  -  [♦fo  *  *10  *  ♦mJ,  .0  *  -1  ^  It, 

1  g  g  /gh  1 


(xj  <  0) 


(3.1.10a) 
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♦lO  -  [4  *  ♦lolx,  -  »  *  »1  *  ^It, 

1  g  -/gh  1 


(xj  >  0) 


(3.1.10b) 

where  we  have  replaced  derivatives  with  respect  to  Xj  by  derivatives  with 
respect  to  t^,  expanding  (3.1.8)  In  the  form: 


<=  ■  «<»)  "’‘if:  -  <=«» "  -I  f:  "g 


(3.1.11) 


'1  *  ""I 

We  now  turn  to  specific  problems.  The  first  case  examined  will  be  that 
of  slow  sway  of  small  amplitude  (0(e) )• 

3.2  The  Solution  for  the  Fast  Motion 


When  the  mooring  Is  somewhat  stiff, 

Kq  “  =  0(e)  (3.2.1) 

the  slow  sway  displacement  (and  the  fast  one)  Is 

X  -  0(e)  (3.2.2) 

which  can  be  expanded  Into  harmonics: 

CO  n 

^  “  I  e”  I  X  exp(-lmu)t)  (3.2.3) 

n“l  m*-n 

The  boundary  conditions  on  the  body  (3.1.2  -  3.1.4)  can  be  Taylor-expanded 
about  the  rest  position,  Sq. 

At  the  first  order,  first  harmonic,  we  get  for  the  short  wave  potential, 

♦li¬ 


on 


(3.2.4) 
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0 


(3.2.5), 


-oj  M  /_  fjjdz] 


(3.2.6) 


Together  with  (2.3.1  -  2.3.3)  and  the  radiation  condition,  these  equa¬ 
tions  define  the  problem  for  which  Is  the  near  field  of  Formally 

these  equations  are  Identical  to  the  equations  for  the  linear,  time  harmonic 
problem  of  a  rectangular  cylinder  swaying  freely  In  regular  waves.  This 
problem  was  solved  by  Mel  and  Black  (1969),  using  a  variational  method.  We 
shall  employ  the  equivalent  Galerkln  method  In  order  to  solve  for  the  dynamics 
of  the  sway  motion  as  well  as  Its  kinematics.  To  start  with,  we  break  the 
problem  down  Into  a  diffraction  problem  and  a  radiation  problem. 

The  diffraction  problem  deals  with  a  fixed  body.  We  break  It  down 
further  into  a  symmetric  and  an  antisymmetric  problem.  The  analysis  Is  con¬ 
fined  to  (x  <  0).  Starting  with  the  symmetric  (S)  problem,  (2.3.4)  becomes: 


S  ,  .  lk(x+B)  ^  _  -lk(x+B)  ^  ^  , 

♦u  ■  ■'s  '  '*  I,  Vn  ' 

n“l 


k  (x+B) 
n 


,  (x  <  -B)  (3.2.7) 


under  the  body,  the  potential  has  the  form 


i|», ,  *  A-  F„  +  y  A  F  cosh  K  x 
0  0  n  n  n 

n“l 


(x  >  -B)  (3.2.8) 


For  the  antisymmetric  (A)  problem: 


-  f  K  /  lk(x+B)  -lk(x+B) 

'*'11  '  Vo^® 


OB  k  (x+B) 

)  +  y  b  e  "  (x  <  -B)  (3.2.9) 


lb,  ,  *  B^  F»  x  +  y  B  F  slnh  K  x 
’^11  0  0  n  n  n 

n=l 


(x  >  -B)  (3.2.10) 
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The  eigenfunctions  under  the  body  satisfy  (2.3.3,  3.2.5)  and  are  given  by: 


1/2  - 

Fq  -  (h-D)  '  ;  -  /2/(h-D)  cos  (z+h),  n  -  1,  2,  ...  (3.2.11) 

Rs  and  are  the  reflection  coefficients* 

There  are  two  continuity  conditions  that  apply  at  the  boundary  between 
the  fluid  region  under  the  body  and  the  region  to  its  left:  both  the  poten¬ 
tial  and  must  be  continuous.  Let  us  define: 

(3.2.12a) 


u,  - 


A  „  .S 
llx  ’  ”  "^llx 


”  'I'T,™  (x  *  -B)  (-h  <  z  <  -D) 


=0  -  0  (-D  <  z  <  0)  (3.2.12b) 

and  Ug  are  the  complex  amplitudes  of  the  horizontal  velocity.  Using 
the  orthonormality  of  (fn)  fil'd  that  of  we  can  express  {fin}»  {A^} 

in  terms  of  Ug,  and  {b^}.  (Bn)  terms  of  U^: 


lk( 


1-Rg)  /jj  ^0  ®0  *  ik(l-R^T  /jj  ^0 


1 

r-  /  u.  f 


(3.2.13a) 


n  k  •*„  S  n 
n  -H 


1 

t)  =  ^  /  U.  F„ 
n  k  „  A  0 
n  -H 


dz  (n  >  1) 


^0 


-H 


K  slnhK  B 
n  n  -H 


(3.2.13b) 

1 

/  U  F  dz  ;  B  =•  - - ir-p-r-  r  U.  F  dz  (n  >  1) 

iu  s  n  ’  n  K_coshK  B  /  A  n  — 

n  — H 


-D 


Ag  does  not  affect,  hence,  cannot  be  determined  from  U3. 

S 

We  may  now  use  (3.2.13)  to  write  down  the  continuity  condition  on 
i|/jj  at  (x  ■  -B).  Eliminating  Rs  (or  R^)  and  keeping  ag  and  bg,  we 
get: 
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(-h  <  z  <  -D) 


(3.2.14) 


Multiplying  (3.2.14)  by  {Fq,...Fjj}  and  integrating  from  -h  to  -D  with 
respect  to  z,  we  get  N+1  equations  for  {BQ.BjKjCoshKj,  . . .  ,BNKf}CoshKNB}. 
Multiplying  (3.2.15)  by  {Fj,...,Ffj}  we  get  N  linear  algebric  equations  for 
{-AjKjSinhKj^B, . . . ,  -ApfKpjsinhK(jB} .  Details  are  given  in  Appendix  A. 

Having  obtained  an  approximation  for  the  velocity  (x  =*  -B)  we  can  determine 
R5,  and  the  force  on  Sg  from  (3.2.13). 

The  diffraction  problem  for  a  wave  incident  from  the  left  is  solved  by 
superposition  of  the  symmetric  and  the  antisymmetric  solutions.  We  have: 
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'd  ■  7  «S  *  "a*  =  ’'d  ■  I  ('s  -  V  =  ^  (3.2.18) 

here  Rj,  T^j,  and  F  are  respectively  the  reflected  wave,  the  transmitted 
wave  and  the  force  due  to  diffraction,  per  unit  incident  wave  amplitude,  and 

=  pirn  /  dz  (3.2.19) 

"o 

The  solution  of  the  radiation  problem  for  sway  is  done  in  the  same  way, 
with  (3.2.12b)  replaced  by: 

“  V  (x  =  -B)  (-D  <  2  £  0)  (3.2.20) 

where  V  is  the  complex  amplltiide  of  the  sway  velocity  (the  sway  potential  Is 
antisymmetric),  see  App.  .''-x  A  for  details. 

Having  found  the  added  mass,  y,  and  the  damping  coefficient,  X,  for  the 
swaying  cylinder,  (3.2.6)  becomes: 

[-0)^  (M+y)  -  luXl  Xjj  *  A  F^  (3.2.21) 

for  an  Incident  wave  amplitude  A.  If  the  left-going  radiated  wave  has  complex 
amplitude  Rj.  for  a  unit  amplitude  sway  we  find  that 

R  =  Rj  +  R^  Xj^/A  =  {R^  +  Rj..F^/[-a>^(M+y)  -iuX]  }  (3.2.22) 

T  -  -  R^  Xjj/A  =  -  R^-F^/[-a)^(M+y)  -  lo)X]  } 

As  a  special  case,  we  shall  be  Interested  In  a  body  which  extends  the 
entire  sea  depth  and  slides  on  the  bottom,  that  is: 

D  -  h  (3.2.24) 

For  this  case,  we  obtained  an  analytical  solution  for  the  combined  diffrac¬ 
tion  radiation  problem  of  (Agnon  &  Mel  1985a).  This  solution  Is  given  In 
Appendix  B. 
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As  another  Interesting  special  case,  we  consider  a  narrow  gap  between 
the  bottom  and  the  body: 

h  -  D  =■  0(£h)  (3.2.25) 

In  view  of  the  narrowness,  in  the  gap  represents  a  uniform  flow  which 

depends  on  the  pressure  gradient  between  its  two  ends: 

-9P/8X  =  -imp  (B,  -h)  -  (-B,  -h)]/2B  =  0(e)  (3.2.26) 

The  effect  of  this  flow,  when  multiplied  by  the  gap  width,  amounts  to  an 
O(e^)  flux.  It  is  negligible  outside  the  gap  since  it  is  equivalent  to  an 
O(e^)  sink-source  pair.  The  potential  outside  the  gap,  is  then  given  to 

order  0(e)  by  the  solution  in  Appendix  B  for  a  sliding  body,  as  if  the  gap 
did  not  exist. 

We  now  turn  to  the  slow  sway. 

3.3  Slow  Sway  of  Small  Aaplltude 

The  slow  potential  fjo,  which  in  the  near  field  is  denoted  by  is 

governed  by  (2.4.4)  and  (2.4.5).  For  the  boundary  condition  on  the  body,  we 
derive  the  second  order  zeroth  harmonic  from  (3.1.2)  -  (3.1.4).  Expanding 
(3.1.2)  about  Sg  we  get: 

I'x  “  "^x  ^^0^  *xx  ^  \  (3.3.1) 

where  <1/^^,  ij/^x  on  the  right  hand  side  are  evaluated  at  Sg  which  are  the 
mean  positions  of  the  side  walls. 

The  second  order  zeroth  harmonic  component  is: 

'*'l0x^  '*'20x  ^'•'llxx  *11  "  ^lOt^  ^0^  (3.3.2) 
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where  ^20^  is  zero,  since 


^20  “ 


To  find  the  dynamic  boundary  condition,  we  note  that: 


(3.3.3) 


d  Y  3  Y 

2  ^10  ^  ^lOt.t,  “  ^ 

at  11 


(3.3.4) 


iJ)^(S  )  =■  4)^(Sq)  +  X  +  0(e  ) 


(3.3.5) 


I  i|)  dz  =  ij)  (z=0)  c  +  0(e  ) 
0 


(3.3.6) 


/  z  dz  -  Y  c 
0 


(3.3.7) 


The  Integral  of  the  pressure  on  S-  will  have  contributions  from  Integration 
between  -h  and  zero  and  from  Integration  between  zero  and  The  contri¬ 

bution  to  the  latter  due  to  the  quadratic  terms  in  the  expression  for  the 
pressure  (2.1.2)  is  0(e^)«  The  only  contribution  at  O(e^)  is  that  evaluated 
in  (3.3.6),  (3.3.7). 

Hence,  the  second-order,  zeroth  harmonic  of  (3.1.4)  is: 

KX^q  =  -p  [/_  P2Qdz  -  P2Qdz  +  6(lu„|,jj  Cj*  +  *)  +  6(gqj  (3.3.8) 


where 


^20  "  '^lOtj  h’^'lll^  (^“'I'tl^^ll  +  *)  +  gz 

(from  the  Bernoulli  equation  (2.1.2)  and  from  (3.3.5)), 

6G  =  G(-B,  0)  -  G(B,  0) 


(3.3.9) 


(3.3.10) 
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Equation  (3.1.3)  simply  gives: 


i1)2q  “0  on  and  (z  =  -h)  (3.3.11) 

z 

We  note  that  the  kinematic  boundary  condition  (3.3.2)  couples  4»jq,  which 
is  a  function  of  long  spatial  scale,  with  ij)20>  which  is  a  function  of  short 
spatial  scale.  The  right  hand  side  of  (3.3.8)  (see  also  (3.3.9))  can  be 
split  into  contributions  due  to  iJ;jg  and  Since  4;^^  is  formally  the  same 

as  4)j2  for  a  regular  wave  train,  the  terms  Involving  have  the  same  form 
as  the  terms  in  the  expression  for  the  steady  drift  force  in  a  regular  wave- 
train.  In  Appendix  D  we  evaluate  these  terms  and  obtain  an  expression  which 
agrees  with  that  obtained  by  Maruo  (1960),  by  Newman  (1967)  and  by  Longuet- 
Higglns  (1977),  using  conservation  of  momentum  and  maintaining  that  there  is 
no  energy  absorbed  by  the  body.  The  contribution  of  the  short  waves  is: 

Q 

pg  |a|^  |r|^  ^  (3.3.12) 

where  A,  the  amplitude  of  the  incident  short  wave  given  by: 

2ioj  fn(0) 

A  =  - - - A_  (3.3.13) 

g  0 

Aq  Is  the  amplitude  of  the  incident  short  wave  potential  $ j ^  and 

(C  =  u/k)  is  its  phase  speed.  Equation  (3.3.8)  can  now  be  written  as: 

Q 

KXio  *  -p[  /_  “  /+  ^  ^  (3.3.14) 

The  form  of  in  the  near  field  will  depend  on  the  specific  geometry 
of  the  body.  Denoting  the  gap  width  by  H,  we  distinguish  among  three  cases, 
according  to  the  width  of  the  gap  between  the  body  and  the  bottom: 
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a)  large  gap:  H  /  h  =  0(1), 

b)  small  gap:  H  /  h  “  0(e), 

c)  no  gap:  H  =  0. 

These  three  cases  are  studied  In  the  following  subsections.  Each  case 
represents  a  different  asymptotic  limit.  Together  they  give  a  picture  of  the 
dependence  of  the  body's  response  on  the  size  of  the  gap. 

3.3.1  The  Large  Gap  Case 

A  geometry  common  In  practice  Is  one  In  which  the  gap  between  the  body 
and  the  bottom  Is  comparable  to  the  water  depth. 

The  variation  of  ij/^g  with  respect  to  (x,  z)  Is  governed  by  Laplace  Equa¬ 
tion  (2.4.4)  and  by  homogenous  Newmann  conditions  on  the  mean  free  surface 
and  the  bottom  (2.4.5).  The  flux  through  Sq,  Sq  Is  O(c^),  because  of 

(3.3.2),  and  the  flux  at  the  matching  boundaries  to  the  far  field 

S*  =  {X,  z|  1  «  kx  «  -h  <  z  <  0}  (3.3.15) 

(see  Figure  3.2),  Is  O(e^)  because  of  (3.1.10). 

Since  the  cross  section  of  the  gap  Is  0(1),  the  flux  through  It,  which 
Is  the  horizontal  gradient  of  the  slow  potential  Is  Vhc^lO*  H  must  be 
matched  to  the  far  field  according  to 

\  e+lO  ”  /H)  “  O(g^)  (3.3.16a) 

Hence,  the  spatial  variation  of  the  slow  potential  In  the  near  field  with 
respect  to  the  short  scale  Is  second  order  In  e,  except  In  the  vicinity  of 
sharp  corners.  This  Is  In  agreement  with  the  0(e)  kinematic  boundary 
condition: 
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'I'lOx  ”  ®  (3.3.16b) 

Thus,  Is  only  a  function  of  time  almost  everywhere  in  the  near  field: 

i|>10  -  (3.3.17) 

Integrating  across  the  hull,  we  get 

/  'I’lOf-  dz  “  j  dz  “  0  (3. ^.18) 

lUtj  lUtj 

^0  0 

The  dynamic  boundary  condition  (3.3.14)  becomes: 

Q 

K  -  pg  |a|^  |r|^  ^  (3.3.19) 

Formally  this  result  is  the  same  as  that  for  the  steady  drift  motion  in  a 
regular  (unmodulated)  wave  train,  and  agrees  with  Newman's  (1974)  theory,  in 
which  the  slow  drift  force  is  approximated  by  the  steady  drift-force  form¬ 
ula. 1 

The  present  theory  allows  us,  however,  to  find  the  long  wave  field,  and 
the  vertical  force  as  well.  To  determine  the  wave  potential  we  match  the 
slow  potential  and  its  gradient  -  on  both  sides  of  the  near  field  -  with  the 
corresponding  quantities  in  the  far  field.  Matching  of  the  potential  gives: 

4ijq(0»  tj)  *^10^^’  ^1^  (3.3.20) 

To  the  far  field  observer,  the  near  field  shrinks  to  a  line  at  (x^  »  0) 
across  which  the  potential  is  connected  via  ij^^g.  Combining  (3.3.10)  and 
(3.3.20),  we  have: 

**  (3.3.21) 

1  In  the  large  gap  case,  large  slow  sway  is  possible  (for  weaker  mooring). 
This  will  be  studied  in  Section  3.4.  The  results  presented  there 
include  (3.3.19). 
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In  order  to  determine  the  unknovm  qjjQ,  we  need  a  further  condition  on  the 
jump  in  3(^jq/3xj.  This  condition  is  related  to  the  horizontal  flux.  At  the 
second  order  zeroth  harmonic,  the  flux  across  S^,  where  kx  >  1  and  x,  <  1  is 

CD  ^ 


+ 

z=0 

S+ 


This  flux  should  be  matched  to  the  far-fleld  flux: 


-h 


3X  |g  + 


dz 


(3.3.22) 


(3.3.23) 


The  overbar  stands  for  averaging  over  a  short  period  2tj/o}  and  short  spatial 
scale  2Tr/k.  We  now  insist  that  (<l>20jj  +  0)  as  (jxj  i  »)  so  that 
0(x,) 

/  <|)2q  dx  =•  <|)2Q  =  0(1)  and  not  0(e  )  as  x  =■  0(l/e) 

0  X 

In  the  near  field,  however,  t20x  contain  a  constant  term.  To  evaluate 

this  term  we  write  the  governing  equations  for  ^20’  From  (2.2.5b),  (2.2.6b), 

(2.2.7b)  and  (3.1.3),  we  get: 

.2  2 

^20  “0  ("h  <  z  <  0)  (3.3.24) 

3x  3z 

9'I'20 

-r^  -  0  (z  =  -h)  (3.3.25) 

3z 


3i|) 


20 


3z 


X  +  *]x  H  gh  Vj^.U 


(z  =  0) 


(3.3.26) 
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where  U  “  (U,V)  Is  defined  by: 

ghU  =  -  +  *]^_o  (3.3.27) 

The  kinematic  condition  on  the  body  is  obtained  from  (3.3.2),  and  is: 

'*’20  “  “^^11  '*’11  ■*■*)+  ^lot  ^0  (3.3.28) 

X  XX  1 

Thus  il/jg  accounts  for  the  short  scale  variation  of  the  slow  potential.  We 
must  now  find  the  last  term  in  (3.3.22).  It  turns  out  that  it  is  not  neces¬ 
sary  to  solve  explicitly  for  <|)2o*  Because  of  Laplace  Equation  (3.3.24),  we 
may  apply  Gauss'  theorem  to  ij/jg  ®  control  volume  that  is  bound  by  the  sur¬ 
faces  (z  =  -h)  and  (z=0)  on  bottom  and  top,  the  surfaces  and  S„  on  the 
right  and  left,  and  the  rest  position  of  body  the  surface  Sq  (see  Figure 
3.2).  We  shall  write  the  flux  of  i\i2q  through  these  surfaces  in  terms  of  U 
and  add  them  up  to  find  the  value  of  the  last  term  in  (3.3.22). 

By  Gauss'  theorem. 


-B 


[/  -  /  +  /  "  /  ]  ‘iz  'P20X  ■^  /  ]  1': 


<  s;  sj  Sq 


20z 


dx 

z“0 


B 

+  /  i|> 

-B 


20z 


dx  -  f 
2  =  -D  - 


20z 


dx 

z  “  -h 


Now,  the  integrals  along  the  free  surface  are 

— B  00  — B  «»  — B  00 

[/  '*’  /  ]  ’l'9n  dz  *  h( J  +  J  )  U  dx  »  h  uj  +  huj  (z“0) 

B  -«  B  *  l-oo  |b 


(3.3.29) 


(3.3.30) 


after  using  (3.3.26). 
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The  integrals  along  the  vertical  walls  of  the  body  are 

(/  -  I  )•“  ♦20«  -  -1/  -  n  nix*  ^  Uf  -  /  )* 


s’*"  s~ 


s~ 


s" 


(3.3.31) 


Since  we  have: 

XX  zz 


/  (X  +  *)dz  =  I-Xjj  /  dz  +  *] 

-D  XX  -D  zz 


(-Xjj  (±B,  0)  +*)  =  (-  -^  ijijj  (±B,  0)  +  *) 


=  -hU(±B) 


(3.3.32) 


(making  use  of  \p  ^  —  {p  at  (2=0)  and  <(i 

11^  g  11  1 


-iwX.  at  (x  =  ±B)). 
11 


When  we  combine  (3.3.23)  and  (3.3.32),  we  get  an  expression  for  the  net 
flux  across  the  body  through  the  gap  (the  body  is  rigid  and  Impermeable): 

/  no  -f  no  b  “  not,  -  “  not,  • L 

„+  X  X  ■‘•B  1  1  -IB 


Because  of  Equations  (3.3.24)  and  (3.3.30)  we  get: 


TD  — D 

/  m  dz  -  J  4^20  dz  =  hU  I  -  hU  I  -  hU  I 

.+  X  X  ■‘•1 


(3.3.34) 


U(±<»>)  are  given  by  (3.3.27): 

gh  U(-»)  =  -  {lo)A*fQ(0)  [ikAQfQ(O)]  +  iaJl*A*fQ(0)  [-ikR  AQfQ(O)]}  +  * 


2mk|AQ|^(l-|R|^)fQ(0) 


(3.3.35a) 
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gh  U(»)  -  -  {iojA*  T*fQ(0)IlkAQ  T  f^CO)]  +  *} 


-  2a)k|AQ|^  |t|^  fj(0) 


(3.3.35b) 


Because  the  body  Is  freely  floating  as  far  as  the  high  frequencies  are  con¬ 
cerned,  there  is  no  energy  transfer  to  the  body;  hence, 

|r|^  +  |Tj^  -  1  (3.3.36) 

That  is  the  energy  carried  away  by  the  short  waves  is  equal  to  that  brought 
in  by  the  incident  waves.  Combining  (3.3.35)  and  (3.3.36)  we  get; 

U(+«)  -  U(-«)  (3.3.37) 


We  are  now  ready  to  match  (3.3.22)  to  (3.3.23).  The  quadratic  terms  in 
the  two  equations  are  equal,  since  is  continuous  in  the  passage  between 
the  near  field  and  the  far  field.  Matching  the  gradients  of  and  ^20* 
substituting  (3. 3. 34)  and  (3.3.36)  we  get: 


♦lOxj 


/  ’f’20  ~  I  ^20  *** 

,+  X  _+  X 


hU(«)  -hU(-«)  -  0 


(3.3.38) 


OD 


From  (3.1.9)  and  (3.3.36)  we  immediately  obtain: 

♦iQ  (3.3.39) 

Matching  the  gradients  of  <^jq  as  they  appear  in  (3.1.10),  and  using  (3.3.38), 
we  see  that  the  locked  wave  contribution  is  eliminated  by  virtue  of  (3.3.39). 
We  are  left  with: 

“♦jo(t)  "  *^10^^^  ^*1  “  (3.3.40) 
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Using  (3.3.39)  we  write  the  right  hand  side  of  (3.3.21)  in  terms  of  ♦iq* 
The  left  hand  size  is  written  in  terms  of  ®a*^fng  use  of  (3.3.40).  The 

result  is  simply: 


(3.3.41) 


Thus,  the  free  long  wave  propagating  to  the  right  is  equal  to  the  long  wave 
locked  to  the  reflected  wave  group.  Substituting  (3.3.40)  and  (3.3.41)  into 
(3.1.8)  and  (3.3.20)  we  get: 


'*’10^*^P  ”  ***10^*^P  “  (xj  =  0  ) 


(3.3.42) 


(|(jQ  at  (xj  “  0)  la  the  same  as  if  the  body  were  absent!  This  is  similar  to 
the  Froude-Krylov  approximation  in  the  linearized  theory  of  long  waves  past  a 
small  body. 

As  an  aside,  when  calculating  the  vertical  slow  drift  force,  there  will 
again  be  terms  due  to  which  are  formally  the  same  as  those  terms  for  a 

regular  wave  and  terms  due  to  the  slow  potential,  as  in  (3.3.8),  (3,3.9). 

The  vertical  drift  force  due  to  iji.g  for  a  general  hull  shape,  is: 


f 

-p  e  / 


4',n^  dx  dy  “  -G  P'1'1  A 


body  surface  1 


1 Ot ,  wp 


(3.3.43) 


where  is  the  water  plane  area.  As  pointed  out  in  the  introduction  to 
this  chapter,  when  the  water  plane  area  is  0(1),  the  heave  displacement  is 
O(e^).  The  contribution  to  the  heave-inducing  force  due  to  and  the 
correction  due  to  ij/^g  are  comparable,  yet  the  slow  heave  is  too  small  to  be 
interesting  in  practice.  When  the  body  is  shaped  as  a  bottle 


wp 


(3.3.44) 


46 


the  heave  amplitude  becomes  0(e),  yet,  from  (3.3.43)  we  see  that  the  effect 
of  ijiir)  on  the  vertical  force  Is  O(e^),  resulting  In  a  negligible,  O(e^), 
effect  on  the  heaving  displacement.  In  contrast  to  the  vertical  motion, 
there  are  cases  where  i|ijq  may  play  an  Important  role  In  the  sway  motion. 

This  happens  when  varies  appreciably  across  the  body  surface,  as  discus¬ 
sed  In  the  following  subsections. 

Numerical  results  will  be  presented  later. 

3.3.2  The  Small  Gap  Case 

When  the  gap  between  the  body  and  the  seabed  Is  narrow,  the  analysis  of 

Section  3.3.1  does  not  apply,  and  the  slow  potential  affects  the  slow  sway. 

This  can  happen  when  a  ship  Is  moored  In  very  shallow  water,  which  Is  only 

slightly  deeper  than  Its  draught  and  Is  encountering  beam  seas.  We  let: 

H  =  h  -  D  =  0(eh)  (3.3.45) 

First,  let  us  show  that  the  heave  and  roll  amplitudes,  both  fast  and  slow, 

are  O(e^)  and  can  be  neglected.  The  body  Is  forced  by  the  Incident  and 

diffracted  waves,  and  Is  restrained  by  buoyancy  and  the  radiated  wave  force. 

For  the  fast  motion,  the  forcing  Is  0(e): 

P  -  -p  -  0(e)  (3.3.46) 

t 

and  It  Is  O(e^)  for  the  slow  motion: 

-» ir  =‘io  ■  ‘lotj  ■ 

The  radiated  potential  for  the  fast  flow  acts  as  an  added  mass;  this  added 
mass  Is  usually  0(1).  In  the  present  case,  however,  the  narrowness  of  the 
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gap  makes  the  heave  added  mass  and  roll  added  moment  of  Inertia  much  larger 
-  0(e“^).  Their  dominant  parts  are  due  to  the  flow  Inside  the  gap.  This 
flow  Is  primarily  horizontal,  since  the  gap  is  narrow. 


Denote  the  fast  heave  and  roll  displacement  by  Z  and  0,  the  associated 

(Z)  (6) 

potentials  In  the  gap  by  and  (f>  ,  and  the  related  added  mass  and  moment 

of  inertia  by  ^2  ^9»  respectively. 

The  equations  of  continuity  for  are; 


-H 


^xx 


(-ioi)  z  =  0 


+  (~iu))  0x  =  0 

^xx 

Since  Z  is  constant  In  x  and  0x  is  linear  In  x. 

(Z  )  ( 0 ) 

(f)  Is  symmetric  and  (p  Is  antisymmetric 


so  that; 


4> 


(Z)  =  _ 


,  Z  x 

im 


(0) 


.  0  X 

^  “  hF 


(3.3.48) 

(3.3.49) 


(3.3.50) 

(3.3.51) 


from  which  we  find  and  l0  to  be: 

B  3 

- - 2~  /  ~  ^  fir  ^  (3.3.52) 

oj  Z 

B  5 

(1)  0  “B 

Since  0)  =  0(1),  and  the  forcing  Is  0(e),  the  fast  heave  and  roll  displacement 
are  0(e/e“M  “  ,  and  may  be  neglected. 

For  the  slow  heave  and  roll,  the  dominant  terms  of  the  restoring  force 
and  moment  are  due  to  buoyancy.  The  inertia  term  being  O(e^).  For  heave  the 
buoyancy  force  Is 
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-2BpgZ  =  0(Z) 


(3. 3. 54) 


while  the  restoring  moment  for  roll  is 

-g  M  m^  0  =  0(6)  (3.3.55) 

where  m^  the  raetacentrlc  height  Is  assumed  to  be  0(B).  Since  the  forcing 
for  slow  motion  Is  O(e^),  we  find  that  the  slow  heave  and  roll,  too,  are 
O(e^).  Therefore,  the  only  first  order  slow  motion  present  is  sway. 

Because  of  the  narrow  gap  under  the  body,  strong  blockage  occurs.  The 
slow  current  under  the  body  becomes  large  (0(c))  and  there  Is  a  large  pres¬ 
sure  gradient  due  to  the  potential  difference  in  between  the  two  ends  of 
the  gap.  Apart  from  the  difference  In  between  the  two  sides  of  the  gap, 
the  flow  In  the  gap  has  only  an  O(e^)  effect  on  mass  flux  which  does  not 
cause  any  additional  effect  on 

In  order  to  determine  ij/^g,  we  define  ipg  to  be  the  sum  of  and  eij»2o* 

*^0  =  ^'*’10  *  ^  ’•’20  (3.3.56) 

(pg  is  the  long  scale  slow  potential  In  the  near  field.  The  slow  potential 
has  a  component  with  short  length  scale:  2Tr/2k,  and  a  component  with  a  long 
scale:  0(e“^  2ir/2k).  The  kinematic  boundary  condition  on  the  body,  (3.3.2), 
combines  both  long  scale  and  short  scale  zeroth  harmonic  potential.  The  long 
spatial  scale  component  can  be  separated  by  integrating  (3.3.2)  along  Sg, 
keeping  only  the  propagating  modes  contribution,  and  not  that  of  the  evane¬ 
scent  modes  that  give  rise  to  short  scale  The  Integral  of  the  last  term 

on  the  left  hand  side  of  (3.3.2)  was  evaluated  in  (3.3.32).  The  long  scale 
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part  associated  with  the  propagating  modes  is  given  by  U(±oo)  at  the  far 
field,  which  was  given  in  (3.3.35)  and  (3.3.37).  We  note  that  U(±<»)  is  just 
Stokes'  drift. 

We  are  now  ready  to  write  the  resulting  kinematic  boundary  condition  on 
the  body.  Separating  the  long  scale  component  of  (3.3.32)  we  get 

j  (X,,  \|)  +  *)dz  =  -hl](±B)  =  -hU(±oo)  on  S“ 

11  ’^llxx  0 

When  substituting  this  result  in  place  of  the  last  term  of  (3.3.1)  averaged 
over  Sq  and  recalling  that  the  difference  between  h  and  D  is  0(e),  we  get  the 
boundary  condition  for  defined  in  (3.3.56),  in  the  form: 

1<0  =  ~  o"  S*  (3.3.57) 

X  1 

i(»Q  also  satisfies  the  rigid  lid  boundary  conditions: 

i|)q  =  0  (z  =  0,  -h)  and  on  bg  (3.3.58) 

z 

Recall  that  both  (J/^g  and  (|;2g  satisfy  (3.3.24).  i|>20z  equal  to  zero 

on  the  free  surface  while  the  short  scale  part  accounts  for  the  inhomogeneous 
term  in  (3.3.26).  With  (3.3.57)  and  (3.3.58)  ijjg  is  the  potential  for  a 
steady  flow  ev  past  a  body  moving  at  a  steady  velocity  Xiot2  ~U(±oo),  with 
a  parametric  dependence  on  t^. 

It  is  known  for  a  stationary  body  in  long  waves  that  the  outer  approxi¬ 
mation  of  the  near  field  in  terms  of  (x"  =  x  -  Xg)  is 

<|>0  ’  QCtj)  +  e(v  -  X^Q^^  +  U)x"  ±  e(v  -  X^^^  +  U)c 
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where  c  is  the  blockage  coefficient.  By  a  Gallllean  transformation,  we  go  to 
the  stationary  frame  of  reference  in  which  the  body  is  moving  at  the  velocity 
Xjotj  “  outer  approximation  of  the  near  field  potential  must  then 

be 

i|»0  ~  Q(tj)  +  [v  -  +  Ul(±ec)  +  vex  (x  -»•  ±«»)  (3.3.59) 

where  Q  is  an  unknown  integration  constant  and  ev  is  the  unknown  current, 

)j)Ox»  at  X  -►  ±"0.  c  is  the  blockage  coefficient  which  is  known  to  be  a  large 
quantity 

c  =  +  O(loge)  “  0(e'S  (3.3.60) 

(cf.  Flagg  and  Newman  (1971)). 

The  outer  expansion  of  the  near  field,  (3.3.59),  has  to  be  matched  to 
the  inner  expansion  of  the  far  field,  (3.1.10).  This  matching  gives: 

Q(tj)  +  [v  -  +  Ul  ec  =  (♦[q  +  ^  ^ 
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+  — 
*10  “  "*10 


(xj  =»  0) 


.J.3.62) 


and 


'10 


10 


I  R  — 

^*10  ^  *1^  ^  *in^ 


'10 


'10' 


T 

g 


+ 

*10 


/gh  *^1 


-  X 


lOt, 


+  U}  (3.3.63) 


where  Cj  =  ec. 

The  large  blocking  by  Che  narrow  gap,  causes  a  jump  in  =  ij;^g  between 
the  two  sides  of  the  block  given  by  (3.3.63),  although  the  gradients  of  the 
slow  potential,  on  either  side  of  the  body,  are  O(e^).  This  jump  is  impor¬ 
tant  to  the  dynamic  boundary  condition  on  the  body,  (3.3.14),  which  can  be 
written  in  terms  of  the  far  field  potentials.  Since  does  not  change 
appreciably  on  either  side  of  the  body.  This  gives: 


^  ^10  °°  P^^*10  ^  *10  "  ^* 


10 


*10 


+ 


pg 


7  C 
2  _& 


(x,  =  0) 


1 


(3.3.64) 


We  now  wish  to  eliminate  <J)^g  between  (3.3.62)  and  (3.3.63).  From  (3.3.62) 
and  (3.3.64)  we  get  after  using  (3.3.37): 


1 


KX 


10 


ph 


+  2<|» 


lOt, 


(3.3.65) 


Differentiating  (3.3.63)  with  respect  to  tj  we  can  substitute  its  right  hand 
side  in  the  left  hand  side  of  (3.3.64)  and  finally  substitute  for  <^tot  ^ 
using  (3.3.65).  This  yields: 

K  „  .  K 


^^*^1*^1  2ph/gh 


^lOtj  2phcj 


10 


[-(-^  |r|^  +  1  -  |r|^) 


It-  S  -f-  +  Ir' 


/gh 


gh  -  C 


2  2h 


at 


(3.3.66) 

2 
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where 


(3.3.67) 


Is  a  factor  appearing  in  the  radiation  stress.  Equations  (3.1.6),  (3.1.7), 
and  (3.1.9)  have  been  used  (details  are  given  In  Appendix  C).  From  (3.3.65) 
we  get  the  free  surface  elevation  associated  with  (J)io»  through  the  Bernoulli 
equation 


4o  '  '*  ♦lOt, 


(3.3.68) 


In  the  limit  as  H  0,  c^  -»■  <»,  the  last  term  on  both  sides  of  (3.3.66)  vanish 
and  the  resulting  equation  may  be  Integrated  with  respect  to  tj,  to  give: 


X 


lOt 


1 


+ 


K 

2ph/gh 


gh  -  C, 


^  S 
2h 


(3.3.69) 


Equation  (3.3.69)  is  Identical  to  the  equation  for  a  sliding  block  (Agnon  and 
Mel  1985)  to  be  discussed  in  Section  3.3.3.  At  the  other  limit,  as  the  gap 
becomes  large,  H  =  0(1),  c^  =  0(e)  and  the  last  term  on  each  side  of  (3.3.66) 
becomes  the  dominant  term,  we  get: 

Q 

K  Xio  -  pg  /  |r|^  |a|2  (3.3.70) 


Equation  (3.3.70)  Is  the  same  as  the  result  (3.3.19)  for  a  wide  gap.  There¬ 
fore,  (3.3.66)  may  be  regarded  as  practically  valid  for  all  gap  widths. 

We  shall  now  solve  (3.3.66)  explicitly  to  study  different  wave  envelope 
forms  and  their  effect  on  the  slow  sway,  as  listed  below: 

a)  A  Steady  Sinusoidal  Envelope 

b)  A  Sudden  Start  of  a  Sinusoidal  Envelope 

c)  A  Gradual  Start  of  a  Regular  Wave 

d)  A  Sinusoidal  Wave  Packet 
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a)  A  Steady  Sinusoidal  Envelope 


Let  the  incident  wave  amplitude  be: 

A  =  a  sin  fitj^  at  (Xj^  =  0)  (3.3.71) 

The  solution  to  (3.3.66)  is  then: 

-21fitj 

XjQ  =  Re[XjQ  e  +  X"^]  (3.3.72) 


where 


,  1  2i£jK  .  K  ,  ,, 

X'  =  D/(-4fi - ,  .  )/2 

*0  /_v  2pnc« 


2ph/gh 


X"  =  —  — S.  IrI^ 

^10  -  2  K  C  ®  '  ' 


1 


(3.3.73) 


(3.3.74) 


and 


„  ,  ^  |„|2  .  ,  .  |,|2, 


£_  s  a^  +  ^  —Si 


C  2hc 


1 


r|2  a2 


(3.3.75) 


In  Figure  3.3,  X^'^  ■  h|x|Qj/a2,  the  normalized  slow  sway  amplitude  is 
plotted  versus  fj  =  fi/u,  the  wave  modulation  ratio,  for  three  values  of  K. 

The  resonance  peak  occurs  near  2fi/a)  =  *^/2phcj.  2phcj  may  be  regarded  as  the 
apparent  mass.  X"q  will  be  discussed  in  Section  3.3.3  and  plotted  in  Figure 

3.5. 


In  all  the  computations  presented  in  this  section  kh  “  1,  M  =  M/ph^  “  1 
and  =  Cj/h  =  1.  The  nomallzed  elastic  constant  is  K  =  KQ/(epgh).  A 
corresponding  choice  of  physical  parameters  would  be,  depth  h  =  10  m  deep, 
the  vessel  width  B  =  10  m,  the  wavelength  X  =  2Tr.l0  s  63  m  which  corresponds 
to  a  wave  period  of  7.2  sec. 
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The  gap  width  is  5e/l  “  0.5  m  if  e  Is  0.1,  and  the  spring  constants  are 
0.1(0. 5,  1,  2)»10«10»1  =  (5,  10,  20)  10^  Newton/m  per  meter  length  of  the 
cylinder.  Each  unit  of  the  abscissa  (o/m)  translates  to  72  sec.  period  of 
the  modulation  envelope  and  each  unit  of  the  ordinate  XJq  translates  to  slow 
sway  amplitude  of  -  1  m  for  incident  wave  of  3  m  amplitudes  (or  6  m  wave- 
height)  and  is  proportional  to  the  square  of  the  waveheight. 

In  nondlmenslonal  form.  Equations  (3.3.73)  through  (3.3.75)  become: 

^I'o  “  ll  -  2lil  oj  K  +  ^)/2  (3.3.73’) 


1  c 
i_ZSL 
2K  C 


2 


D 


=  -ifi  0)  (C 
- -g 


(3.3.74') 

(3.3.75’) 


where 


and 


C  =  C  //gh  ,  C  =  C//gh 
-g  g  ~ 


(I) 


=  oj/h/g 


b)  A  Sudden  Start  of  a  Sinusoidal  Envelope 


Let 


A  =  H(t j  )  a  sin  jjt ^ 


(3.3.76) 


where  H  is  the  Heaviside  step  function.  The  solutions  of  the  homogeneous 
equation  corresponding  to  (3.3.66)  are; 
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,  “1*^1  “2*^1 
Xio  «  e  ,  e 


with 


^  -K  ^  f  K _ K  ,1/2 

4ph/ih  16p^h^gh 


(3.3.77) 


and  we  assume  that  the  square  root  does  not  vanish  so  that  oj  *  0.2  (when  = 

“1  ”  ®2  solutions  are  e®l^l  and  t^  e^l*"!  ,  the  results  are  not 
qualitatively  different  In  any  significant  manner).  In  nondlmenslonal  form. 
Equation  (3.3.77)  becomes: 

2  8K  1/2 
^1’2  “  ~ 

With  Xjg(O)  =  Xioti^O)  =  0  we  get: 


*in  -  - 


“2  ®  -  Oj  e 

^®2  ~ 


*  *l'o  I' 


Ojti  ®2^1 

-2lfitj  (02  +  2ln)  e  -  (oj  +  2ln)  e 

®2  ~  ®1 


(3.3.78) 


where  XJ^  and  X"q  have  been  defined  In  (3.3.73)  and  (3.3.74) 


c)  A  Gradual  Start  of  a  Regular  Wave  Train 


A  =  a  (H(tj)  H(|-  -  tj)  sin  n  t^  +  H(tj  -  |-)]  (3.3.79) 


Then  Xjg  Is  given  by  (3.3.78)  for  0  ^  t  j  Tr/2n  and 


ai(ti-ir/2n) 


R.  |2X  [1  - 


-  a,  e 


a2(ti— ir/2n) 
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+  [a2  Xj^Q  (Tr/2n)  -  X^q^  (ir/2n)] 


OiCti— ir/2n) 


®2  ~  “l 


(3.3.80) 


a2(ti-Ti/2n) 

-(<5  >io  -  Jjotj 

as  (tj  +  oo)  X.Q  tends  to  2X"  which  is  the  steady  drift  force  and  Is  also  the 


value  given  by  (3.3.19)  and  (3.3.70), 


d)  A  Sinusoidal  Wave  Packet 


A  =  H(tp  H(^  -  tj)  a  sin  n  t^ 

then  XjQ  is  given  by  (2.92)  for  0  ^  ir/o  and 

a^(t^— n/ O) 

V--~— 

1  2  1 


“  [oi  Xj^Q  (ir/n)  -  Xjg^  (ir/fi)] 


n) 


®2  "  “l 


(3.3.81) 


(3.3.82) 


In  nondimenslonal  form.  Equations  (3.3.78),  (3.3.80),  and  (3.3.82)  are  modi¬ 
fied  by  changing  (<*^2^  ^  to  Q,  and  Xj^  and  ^lo  and  2ijQ* 

In  Figure  3.4,  a  typical  example  of  the  displacement  X^q  from  (3.3.78, 
3.3.80,  and  3.3.82)  is  plotted  versus  nondimenslonal  slow  time,  fit^.  There 
Is  a  considerable  overshooting  before  a  steady  or  a  quasi-steady  state  Is 
reached.  For  the  gradual  start  of  a  regular  wavetraln  (Case  C),  the  over¬ 
shoot  can  be  an  order  of  magnitude  larger  than  the  steady  state  value,  where 
♦iQ  has  no  effect.  This  transient  drift  force  also  alternates  direction,  and 
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In  Plot  a)  K/pgh  =  I  and  in  Plot  b)  K/pgh 


60 


becomes  large  when  directed  both  to  the  lee  of  the  body  and  in  the  opposite 
direction.  The  modulation  ratio  Is  1  (£}  =  eoi)  so  that  the  envelope  period  Is 
72  sec  for  the  parameters  given  In  Case  a.  The  normalized  elasticity 
constant  Is  1. 


3.3.3  The  Case  of  No  Gap 


Finally  we  examine  the  case  of  complete  blocking.  There  Is  no 
connection  between  the  two  sides  of  the  body,  which  slides  on  the  seabed. 

This  case  was  studied  by  Agnon  &  Mel  (1985a).  The  dynamic  boundary 
condition,  (3.3.57),  can  now  be  directly  matched  to  (3.1.10),  the  far  field 
Inner  expansion,  since  there  Is  no  flow  underneath  the  body,  and  the  gradient 
of  4»jq  Is  uniform  throughout  the  near  field.  The  kinematic  boundary 
condition  Is  now; 


lOt, 


♦ 


lOx, 


(x. 


0-) 


(3.3.83) 


1  *'^"1 

Combining  (3.1.10),  (3.3.64),  and  (3.3.83)  we  can  eliminate  In  the  same 
way  that  (3.3.66)  Is  obtained  In  Appendix  B  to  get: 


lOt, 


10 


M^l' 

/— U  I 


1-  R^] 


ec 


S  A 


gh 


0 


(3.3.84) 


1  2ph/gh  /gh 

which  Is  (4.31)  of  Agnon  &  Mel  (1985a)  where  the  notation  for  S  Is  different. 

This  Is  the  differential  equation  for  the  slow  sway  of  the  body.  The  first 

term  signifies  the  effect  of  radiation  damping,  the  second  term  the  effect  of 

mooring,  and  the  third  term,  the  wave  momentum  flux.  Because  of  the  slow 

motion.  Inertia  Is  Ineffective. 

The  general  solution  to  (3.3.84)  Is  easily  found: 

C  ^  9  ^ 

X,n^t,)  -  -  +  (1-R  )]  - S 

✓gh  2h(gh 


S  “St,  t,  , 

-^e  ^  |Ap(-CgT)e^^  dx  (3.3.85) 
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where 


B  H  --  —  (3.3.86) 

2ph/gh 

We  can  see  Immediately  that  if  the  mooring  stiffness  K  increases,  it 
takes  less  time  for  the  slow  motion  to  reach  equilibrium  under  persistent 
forcing,  or  to  decay  when  forcing  is  removed.  Let  us  specify  the  incident 
envelope  to  be  such  that 


a  sin  fitj  (0  <  tj  <  T^) 

A  =  {  (3.3.87) 

0  otherwise 

where  a  is  the  maximum  amplitude  of  the  Incident  wave  potential;  then  the 
slow  displacement  is 


(3.3.88) 


^10(1  -  e  -  Reid  -  ^)  (e 

P 

-B(t,  -  T  ) 

Xio(Ti)  e  (tj  >  Tj) 


-  e  )]}  (0  <  tj  <  Tj) 


(3.3.89) 


where 


C  /gh  C  2 

XjQ  =  R  >  0  (3.3.90) 

/gh  gh  - 

o 

A 

Note  that  Xjg  is  Independent  of  fi,  is  inversely  proportional  to  the  spring 
constant  K,  and  is  positive. 

Let  us  consider  three  sub-cases: 


a)  A  Sinusoidal  Envelope: 

If  we  let  Tj  +  00  in  Equation  (3.3.87),  then  the  incident  wave  envelope 
is  sinusoidal  for  all  tj  <  oo.  The  limit  for  quasi-steady  state  ntj  >  1  is 

Xio(ti)  -  -XjQ  +  Xjq[1  +  cos  2n(tj  -  i)  Jltj  oo  (3.3.91) 
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where  t  is  the  phase  angle: 


tan  T 


(3.3.92) 


-  2n 


The  ratio  g/fi  measures  the  importance  of  mooring  relative  to  radiation  damp¬ 


ing.  Thus  XjQ  oscillates  about  the  mean  In  Figure  3.5,  the  dimen- 

*  ^  2 

sionless  2iio  “  ^10  plotted  for  three  values  of  M.  For  comparison 

the  first-order  high-frequency  displacement  Xjj  is  also  shown.  Relative  to 
the  mean,  the  amplitude  of  the  oscillatory  part  is 

*10  (3.3.93) 

which  decreases  as  g/ft  increases.  Thus  for  tighter  mooring,  the  body 

wanders  less.  A  typical  history  of  XjQ(t),  normalized  by  a^/h,  is  shown  by 

curve  a  in  Figure  3.6.  Throughout  Figure  3.6  we  have  chosen  the  modulation 

frequency  to  be  e  times  the  short-wave  frequency,  namely  £i  ■  w.  We  have 

also  set  M  *  1  and  Kg  =  epgh. 

“R  — T 

The  reflected  520  ®nd  transmitted  ^20  travel  outwards  at  the  group 
velocity.  Their  amplitudes  depend  only  on  R  and  are  plotted  in  Figure  3.7a 
for  three  different  values  of  M.  The  radiated  long  waves  which  travel  at 
the  speed  /gh  depend  further  on  K  through  XjotjJ  their  amplitudes  are 
equal  by  virtue  of  Equation  (3.3.65)  and  are  plotted  in  Figure  3.7b  only  for 
K  -  pgh. 

b)  A  Constant  Envelope: 

If  we  let  the  incoming  wave  have  the  form 

a  sin  nt.  (0  <  t.  <  Tr/2fi) 

A  -  (  (3.3.94) 

a  i-n/lQ.  <  tj) 
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Figure  3.5  Normalized  Displacement  Amplitudes  of  the 


Incident  Envelopes:  a)  Sinusoidal  Envelope  Started  From  Zero, 
b)  Uniform  Envelope  Started  From  Zero,  c)  Pulse  Envelope. 


Figure  3.7  Sliding  Block:  (a)  Normalized  Amplitudes  of  the  Long  Waves 

T 

Accompanying  the  Incident  Group  C20>  Reflected  Group 

R  T 

^20  and  the  Transmitted  Group  Z20’  Normalized  Amplitude 

of  the  Radiated  Long  Waves  Travelling  at  the  Speed  /gh  to 

the  Right  Left  ^2  0*  Normalization  Length  is 

a^/h,  K  =  1. 


we  get  in  Equation  (3.3.84)  a  constant  forcing  term  when  t^  >  ir/Zj}  and  the 


solution  tends  asymptotically  to  X^q  =  ~^^10  follows: 

-B(t.  -  n/2n) 


(3.3.95) 


For  0  <  tj^  <  Tr/2fi  the  solution  Xj^Q(tj)  is  still  given  by  Equation  (3.3.89). 
Note  that  X^g  is  negative,  implying  that  the  drift  displacement  is  opposite 
to  the  direction  of  the  incident  waves. 

When  a  cylinder  is  floating  or  immersed  so  that  fluid  is  allowed  to 
pass  from  one  side  to  the  other,  the  drift  force  due  to  normally  and  stead¬ 
ily  incident  waves  is  known  to  be 


2 

pga 


(3.3.96) 


and  is  in  the  same  direction  as  the  incident  waves  (see,  e.g.,  Longuet- 
Hlgglns,  1977).  The  corresponding  displacement  is  just 


2X"q  -  pg  a^|R^|  Cg/CK  (3.3.97) 

A 

It  is  different  both  in  sign  and  in  magnitude  (see  Figure  3.5)  from  -2Xjg, 
the  steady-state  limit  of  X^g  (see  (3.3.95)).  Does  this  contradict  our 
result? 

Observe  first  that  Equation  (3.3.97)  can  be  gotten  from  Equation 
(3.3.64)  if  one  sets 


0 


(3.3.98) 


In  the  long-time  limit  of  a  uniform  wavetraln,  the  long  wave  described  by 
(()jQ  becomes  a  steady  current.  For  a  two-dimensional  body  that  does  not 
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prohibit  steady  flow  of  water  from  one  side  to  the  other.  Equation  (3.3.98) 
holds  when  viscous  effects  are  ignored.  Now,  in  the  present  example, 

Ac^lOtj  does  not  vanish;  indeed,  it  also  contributes  a  new  part  in  the  jump 
of  mean  sea  level  across  the  body,  since 


gA^zo  =  (3.3.99) 

The  second  part  a[  ]  above  can  be  calculated 

tV-  (3.3.100) 

which  is  precisely  the  mean  sea  level  change  across  a  two-dimensional  body 
if  fluid  passage  is  allowed  (Longuet-Hlggins,  1977). 

To  confirm  our  result,  let  us  give  a  more  physical  derivation  of  the 
steady  drift  force  on  a  fixed  plate  which  seals  off  all  communications 
between  two  fluid  regions  (x  >  0)  and  (x  <  0).  We  assume  that  the  envelope 
of  the  Incident  wave  train  grows  steadily  from  zero  to  a  constant  value. 

The  fluid  on  the  right  is  never  disturbed.  In  the  fluid  on  the  left  the 
steady-state  spatial  average  of  the  normal  radiation  stress  in  the  standing 
wave  is  known  to  be 


S 

XX 


2E(-^  -  2  )  =  2SE 


where  E  =  y  pg  a 


2 


(3.3.101) 


The  steady-state  mean  set-down  is  also  twice  that  of  the  Incident  progres¬ 
sive  wave,  so  that  the  corresponding  hydrostatic  pressure  is 

A 

l>PgC,n  - - -  (3.3.102) 

1  -  C^/gh 

The  sum  of  the  two  gives  the  drift  force  on  the  plate  as  tj 
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1 


hpgC 


+  s  =  -2SE  r- 

20  XX  '■ 


.  -  C^/gh 


-  1)  -  -2SE(- 


(3.3.103) 


gh  -  c; 


On  the  other  hand.  If  we  take  t^  -*■««.  In  Equation  (3.3.95)  and  R  *  1  In  Equa¬ 
tion  (3.3.90),  we  get  the  steady-state  displacement  whose  product  with  K  is 
precisely  equal  to  the  steady  drift  force  given  by  Equation  (3.3.103).  Thus 
the  mean  set-down  on  the  up-wave  side  is  responsible  for  the  negative  drift 
force. 

In  Figure  3.6  the  typical  transient  motion  ^^^(tj)  is  plotted  as  curve 
b  for  the  case  where  the  Incident  envelope  becomes  uniform  after  the  first 
peak  at  fit^  =  Tr/2. 


c)  A  Pulse  Envelope: 


Let  the  pulse  envelope  have  the  total  duration  =  x/fi.  After  the 
pulse  expires  (jlt^  >  ir),  slow  sway  of  the  body  gradually  attenuates  as  shown 
by  curve  c  in  Figure  3.6.  The  rate  of  attenuation  Increases  with  K  through 
6;  see  (3.3.89).  The  maximum  of  ^^^(t^)  lags  behind  the  peak  of  the  inci¬ 
dent  wave  envelope  at  JJt^  *  Tr/2.  This  is  due  to  the  time  constant 
1/6  «  1/K. 

To  see  the  effect  of  the  mooring  force  we  plot  in  Figure  3.8  the  effect 
of  6/£2  for  the  same  wave  packet.  For  smaller  6/n  (or  K),  the  maximum  dis¬ 
placement  is  greater  but  is  realized  later. 

The  long  waves  are  particularly  interesting  for  this  case.  As  shown  in 
Figure  3.9,  C20*  '^hlch  accompanies  the  Incident  wave  packet,  is  a  solitary 

depression  (set-down)  travelling  at  the  group  velocity.  The  reflected  and 

"*R  “T 

transmitted  wave  packets  are  also  accompanied  by  set-downs  ^20  ?20 

amplitudes  are  simply  related  to  the  reflection  coefficient  R.  Ahead  of 
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Figure  3.8  Effects  of  Elastic  Mooring  Constant  on  the  Transient  Slow 

Displacement  X^q  of  a  Sliding  Block  Due  to  a  Wave  Packet  for 
M  =  ph^  and  kh  =  1.25. 
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Figure  3.9  Scattering  and  Radiation  of  Long  Waves  Due  to  an  Incident 
Wave  Packet;  kh  =  1.25  M  =  1 ,K  =  pgh. 


them  are  the  two  long  waves  travelling  at  speed  /gh.  Because  X^^Ct)  is 
pulse-llke,  X^otj  C20  have  points  of  inflection  (cf.  (3.3.65). 

3.3.4  Comments  on  the  Steady  State  "Set  Down"  and  Drift  Current 

We  wish  to  Interpret  the  slow  potential  in  the  case  of  a  uniform 

wave  train.  It  is  known  that  for  a  uniform  wave  train  there  can  be  an 
arbitrary  second  order  mean  free  surface  displacement  and  mean  current  (Mei 
(1983)).  This  can  be  written  as: 

♦  10  “ 

where  a  and  b  are  undetermined.  If  we  let  the  waves  propagate  onto 
otherwise  calm  water,  with  the  envelope  of  the  waves  given,  say,  by 
(3.3.79): 

A  *  a[H(tp  (3.3.104) 

or  any  other  transient  start  of  a  uniform  wave  train,  we  find  from  (3.1.9) 
and  (3.3.13)  that  ({ijq  will  evolve  to  the  form: 

♦10  “  (3.3.105) 

gh/Cg  -  1  g  g 

once  the  waves  have  reached  their  full  amplitude. 

The  usual  practice  in  treating  uniform  wave  train  problems  is  to 
redefine  the  origin  so  that  the  mean  set  down  is  zero.  We  have  seen  that 
the  initial  value  approach  resulted  in  different  mean  water  levels  on  the 
two  sides  of  the  sliding  block  contributing  to  a  negative  drift  force. 
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On  the  rotating  earth,  Coriolis  effect  will  affect  the  steady  state  drift 
current  as  shown  by  Ursell  (1950).  The  Lagranglan  drift  Ujq  of  (2.4.16)  will 
approach  zero. 

However,  reaching  this  steady  state  lasts  on  a  time  scale  comparable  to 
the  Coriolis  period,  which  is  many  hours  long,  and  is  much  longer  than  the 
time  scale  considered  by  us,  hence  we  may  neglect  rotation  in  the  present 
study. 


3.4  Large  Aaplltude  Sway 

In  many  practical  situations,  the  mooring  is  very  weak  and  the  hydrody¬ 
namic  blocking  is  small.  This  results  In  a  slow  sway  displacement  much 
larger  than  the  wave  amplitude.  It  is  no  longer  feasible  to  carry  out  the 
Taylor  exoansion  about  the  rest  position  of  the  body;  a  different  approach  is 
required. 

Trial tafyllou  (1982)  has  used  multiple  time  scales  and  a  moving  coordi¬ 
nate  system  to  treat  this  problem,  arriving  at  decoupled  equations  for  the 
fast  and  slow  motions.  He  did  not  account  for  the  long  spatial  scales  asso¬ 
ciated  w  th  the  long  period  motion.  As  a  result  the  diffraction  of  the 
locked  Icng  waves  and  the  damping  by  radiation  of  long  waves  were  not 
considered. 

We  carry  out  the  multiple  scale  expansion  in  both  time  and  space  for  the 
present  problem,  with  a  view  to  examining  forces  due  to  long  wave  diffraction 
and  radiation  damping. 

We  assume, 

H  -  0(h)  (3.4.1) 

then  there  is  no  large  blocking,  and 

Kq  -  6%  -  0(e2) 
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so  that 
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0 


(-h  <  z  <  0) 


(3.4.11) 


(z  »  0)  where  o  ■  —  (3.4.12) 

g 


(z  “  -h)  and  on 


(3.4.13) 


-0)  M  Xjj  -  loip  [  /  dz  -  /  dz] 

O  I  “  o  i"*" 

^0  ^0 


(3.4.14) 


where  Sq“  are  the  mean  body  surface  moving  In  only.  These  equations  are 
formally  Identical  to  those  for  In  Se^'tlon  3.2  and  are  also  formally  the 
same  as  Che  equations  for  In  a  regular  wave  train  with  no  slow  motion. 
Therefore,  the  solution  Is  the  same  as  that  obtained  In  Section  3.2  for  a 
body  In  regular  waves.  This  means  that  the  drift  force,  to  the  leading  order 
(0(6^))  la  given  by  the  "steady  drift  force"  formula.  Equation  (3.3.12).  We 
point  out,  however,  that  at  the  next  order  (0(e^))  there  Is  a  contribution 
associated  with  Che  body  velocity  This  contribution  can  be  computed 

using  the  theory  of  a  moving  ship.  Only  a  part  of  the  effect  of  the  current 
can  be  attributed  to  Che  0(e)  Doppler  shift  of  Equation  (3.4.8).  The 
computation  becomes  more  Involved.  We  refer  to  Huljsmans  and  Hermans  (1985) 
who  solved  a  similar  asymptotic  problem  for  a  slowly  moving  ship  In  very  deep 
water.  They  found  that  the  leading  order  correction,  for  a  small  forward 
speed.  Is  proportional  to  that  speed.  We  do  not  calculate  Chat  effect  and 


merely  write  It  symbolically  In  the  form: 


pg  X'  e^|A|^  eXqj.  -  O(e^)  X’ 
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for  the  leading  order  change  in  the  drift  force.  This  change  is  expected  to 
be  proportional  to  the  forward  speed  (or  equivalently  the  current).  Since 
the  changes  occur  in  the  second  order  potential,  numerical  value 

of  X'  is  expected  to  be  0(1).  It  might  be  larger  at  specific  frequencies, 
for  some  bodies.  We  shall  see  that  this  damping-type  term  can  be  of  the  same 
order  as  the  damping  due  to  long-wave  radiation  and  should  thus  be 
considered.  The  physical  origin  of  the  present  term,  however,  is  due  to 
interaction  of  the  body  motion  with  the  short  waves. 

Because  the  forward  speed  is  small,  0(£),  the  Doppler  shift  effect  on  i|;^^ 
is  negligible.  The  slow  motion  is  related  to  t|;Q,  which  is  a  function  of  the 
long  spatial  scale  x^.  Relative  to  that  scale,  Xq  Is  a  small  quantity  and 
the  change  of  mean  position  does  not  affect  the  slow  motion  hydrodynamics  at 
the  leading  order. 

We  may  therefore  study  the  slow  potential  in  the  rest  coordinate  system 
(x,  z,  t).  This  is  different  from  the  approach  taken  by  Triantafyllou 
(1982).  The  equations  for  the  near  field  potential  ij»Q  are: 

=0  (-h  <  z  <  0)  (3.4.15) 

ipQ  =  0  (z  =  0,  -h)  and  on  b^  (3.4.16) 

z 

The  velocity  of  the  body  is  larger  than  in  the  small  amplitude  slow 
motion  case,  and  is 

e  X  =  0(e)  (3.4.17) 

1 

We  discard  from  (3.3.2)  the  second  order  terms 

e^u  =  e^['p20x  ^  ^*^11  ^11*  *  *  O(e^)  (3.4.18) 
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To  the  first  order  0(e)  the  slow  kinematic  boundary  condition  is 

i(»10  =  o..  S"  (x  =  ±B  +  Xq)  (3.4.19) 

The  dynamic  boundary  condition  is  obtained  in  the  same  way  as  (3.3.14) 
was  derived  for  the  small  amplitude  slow  sway.  Including  the  body  Inertia, 
we  write: 

2  2  ^ 

MXq^  +  K  Xq  =  -p  [  /  dz  -  /  dzlj.  +  pg  |a|  |r|  ^  (3.4.20) 

^  ^  s“  s'*"  ^ 

The  slow  potential  i|/jj  can  be  separated  into  two  parts:  The  first  part 

is  a  function  of  t ^  only  and  is  the  near  field  of  the  slow  potential  of 

Section  3.3.1  given  by  (3.3.38)  through  (3.3.42).  Its  corresponding  far 

field  ^10^  consists  of  the  locked  long  waves  which  travel  with  the  short  wave 

groups  (I,  R,  and  T)  and  of  the  free  diffracted  wave  In  particular, 

( 3  ) 

satisfies 

'I'lOx  “  °  ^0  (3.4.21) 

Just  as  in  (3.3.16b). 

(2) 

The  remaining  part  of  the  slow  potential,  i|;q  ,  is  related  to  the  radi¬ 
ated  wave,  satisfying  (3.4.19)  and  radiation  conditions  that  state  that  the 

(2) 

waves  are  outgoing  from  the  body.  i|)g  is  simply  the  solution  to  the  linear 

(2) 

problem  for  wave  radiation  by  the  slow  sway  motion.  The  order  of  t|iq  and 

(2) 

the  corresponding  ^g  will  be  determined  later. 

The  present  case  is  different  from  the  case  studied  in  Section  3.3.1  in 
that  the  radiation  problem  is  Important,  due  to  the  larger  slow  sway  velo¬ 
city.  In  addition  to  <|>ig\  we  have  that  is  the  potential  of  radiation 

due  to  the  slow  sway.  This  potential  will  have  the  effect  of  an  added  mass  y 
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and  a  damping  coefficient  eX.  These  are  derived  in  Appendix  E  by  using 
matched  asymptotics. 

(2) 

In  the  far  field  (t)Q  is  given  from  (E.14)  by: 


±  e  Re 


-2in(t^ 


} 


(xj  >  0) 
(xj  <  0) 


(3.4.22) 


{ 1 )  3 

The  force  due  to  ijjjQ  is  0(e  )  as  shown  for  the  identical  in  Section 

(2) 

3.3.1.  The  force  due  to  iJiq  is  given  in  terms  of  the  added  mass  and  damping 
coefficient  at  the  low  frequency  2j2: 

+  X  (3.4.23) 


The  introduction  of  the  damping  coef f licient ,  eX,  is  made  possible  by 
accounting  for  a  free  wave  in  the  far  field.  This  wave  is  radiated  by  the 
slow  sway,  carrying  away  energy,  which  results  in  damping.  When  c  =  0(1),  it 
is  seen  from  (E.17)  that  eX  is  in  fact  O(e^).  Multiplied  by  Xot^  it  will 
yield  an  O(e^)  contribution,  while  the  rest  of  the  terms  in  the  equation  of 
motion  are  O(e^).  However,  near  resonance,  where  the  inertia  and  elastic 
force  are  nearly  balanced,  damping  is  important  and  serves  to  make  the  reso¬ 
nance  finite.  This  effect  was  not  accounted  for  in  previous  theories,  where 
a  rigid  lid  approximation  results  In  using  a  zero  frequency  limit  for  the 
slow  sway  which  excludes  radiation  damping  (Trlantafyllou  (1982)).  The  equa¬ 
tion  of  motion  for  the  slow  sway  is 

<"  *  ")  *  *  *0t,  *  "  *0  ■  '>*1*1^  <1“!"  f 

where  the  last  term  on  the  right  hand  side  indicates  the  third  order  effects 
for  which  the  details  have  not  been  worked  out.  Note  that  terms  in  (3.4.24) 
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can  be  of  different  orders  of  magnitude.  In  particular  y  and  X  depend  on  the 
gap  width  H  which  Influences  the  blockage  coefficient  c.  If  H/h  Is  0(1), 
then  c  =  2BD/h  =  0(1).  It  follows  from  (E.17)  that  X  is  0(e)  and  y  is  0(1). 
Damping  Is  small.  If  the  mooring  constant  Is  assumed  to  be  O(e^)  then  Equa¬ 
tion  (3.4.24)  predicts  unbounded  resonance  at  leading  order.  To  obtain 
finite  response  at  resonance  it  is  necessary  to  calculate  all  the  O(e^)  terms 
which  are  tedious  to  work  out.  It  Is  seen  from  (3.4.22)  that  when  c  is  0(1) 
the  radiated  potential  is  0(e)  and  is  comparable  to  . 

There  is  a  situation  where  (3.4.24)  can  be  used  without  adding  higher 
order  corrections.  This  is  when  c  is  0(1, Ve)  which  corresponds  to  an  inter¬ 
mediate  gap.  We  also  choose  intermediate  mooring  =  0(e^/^).  From 
(E.17)  we  see  that  y  is  0(l//e)  and  eX  is  0(e).  Equation  (3.4.24)  is  still 
at  order  and  Xq  decreases  to  0(/e).  The  damping  term  eXXo^^,  increases 
to  0(e^/^).  The  damping  due  to  the  modification  of  the  short  waves, 
eX'Xot^,  becomes  0(e^/^)  (since  Xg  decreased  to  0(/e)). 

In  the  results  presented  we  shall  use  c  =  0(1 //e)  and  neglect  the  higher 

(2)  (2)  — 

order  effect  of  X'.  <fro  ^  '^'O  0(/e)  as  seen  from  Equation 

(3.4.22)  and  are  larger  than 

In  Figure  3.10  the  value  of  the  normalized  =  hXg/a^,  the  slow  sway 
amplitude  for  a  periodic  envelope  is  plotted  versus  the  modulation  ratio  il/ai 
for  various  choices  of  the  parameters.  It  is  clear  from  (3.4.24)  and  (E.17) 
that  increasing  c  as  well  as  increasing  K  will  reduce  ^g.  In  Figure  3.11 
transient  responses  corresponding  to  the  inputs  of  Cases  b,  c,  and  d  of 
Section  3.3  are  plotted.  Equations  (3.3.71)  through  (3.3.82)  can  be  carried 
over  directly  to  the  present  case,  with  Xjg  replaced  by  Xg.  The  only 
modifications  are: 
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Figure  3.11 


887 


1)  defined  by  (3.3.73)  and  (3.3.75)  is  replaced  by 

1  ^  9 

X(^  =  -  j^|S./[-(M  +  p)4fi  -  2lfiX  +  K)  (3.4.25) 

2)  aif2  (3.3.77)  is  now  given  by: 


-  ~  4K(M  +  u) 

“l’2  ~  2(M  +  y) 


(3.4.26) 


The  blockage  coefficient  corresponds  to  a  cylinder  that  is  10  m  wide  and 
has  a  7 .5  m  draft  in  water  of  10  m  depth.  The  spring  constant  is  e  (0.1) 
times  that  of  Section  3.3  and  a  unit  sway  motion  is  1/e  (10)  fold  larger:  10 
m  per  unit  for  3  m  Incident  wave  amplitude.  We  turn  to  the  conclusions  of 
the  present  chapter. 


3.5  Conclusion 


The  response  of  a  floating  offshore  structure  in  slowly  varying  waves 
depends  on  body  geometry,  mooring  stiffness  and  wave  characteristics.  We 
have  directed  most  of  our  attention  to  the  slow  sway  motion,  which  is 
expected  to  be  the  largest,  since  it  is  not  affected  by  the  buoyancy  force 
and  moments.  By  a  multiple  scale  analysis  we  have  been  able  to  separate  the 
low  frequency  part  of  the  second-order  fluid  motion  from  the  high  frequency 
part.  By  the  notion  of  near  and  far  fields,  we  have  shown  (i)  that  the  slow 
drift  problem  can  be  solved  without  solving  explicitly  for  any  second  order 
potential  and  (ii)  that  the  slow  motion  in  the  near  field  is  accompanied  by 
the  propagation  of  long  waves  in  the  far  field. 
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The  sway  response  of  the  two  dimensional  body  depends  on  the  so  called 
blockage  coefficient,  which  measures  the  obstruction  by  the  body  to  the  long 
period  current  around  it,  and  on  the  mooring  stiffness. 

When  the  mooring  is  not  too  weak,  the  slow  sway  is  comparable  to  fast 
sway,  this  was  discussed  in  Section  3.3.  In  particular,  when  the  blockage 
was  not  large,  the  body's  response  was  passive,  in  agreement  with  Newman's 
(1974)  generalization  of  the  steady  drift  formula.  We  found  that  the  long 
wave  field  is  Independent  of  the  slow  motion.  These  results  are  expected  to 
be  valid  for  three  dimensional  bodies  that  have  0(1)  dimensions,  since  the 
blockage  of  a  uniform  flow  by  such  a  body  is  small,  regardless  of  the  gap 
under  it;  this  is  because  the  fluid  can  flow  around  the  three  dimensional 
body. 

When  the  blockage  coefficient  is  large  it  has  the  effect  of  large  added 
mass  and  damping  coefficients.  The  differential  equation  of  slow  motion 
becomes  second  order  and  resonance  occurs.  The  response  for  a  transient  may 
exhibit  considerable  overshooting.  This  case  was  solved  analytically.  For  a 
long  vessel  with  a  narrow  clearance  underneath,  we  expect  a  similar  blockage 
effect  to  make  the  slow  potential  dynamically  important.  Finally,  we 
presented  results  for  a  sliding  block  where  the  blockage  is  infinite.  No 
resonance  occurs  in  this  case. 

When  the  mooring  is  very  weak  and  the  blockage  is  not  large,  the  slow 
sway  becomes  0(1).  We  found  two  types  of  damping:  one  is  due  to  the  radia¬ 
tion  of  long  waves  and  the  other  is  due  to  the  "forward  speed”  effect  of  the 
slow  sway  on  the  short  waves.  Both  these  effects  are  small.  The  long  wave 
radiation  was  computed  and  the  slow  sway  was  found  for  Intermediate  gap 
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width.  When  the  gap  Is  0(e)  or  smaller  the  slow  sway  becomes  0(e)  even  if 
the  mooring  is  weak,  since  the  added  mass  and  damping  coefficient  become  very 
large. 

The  heave  motion  too,  was  discussed.  We  have  distinguished  between  a 
bottle  shaped  body  and  a  regular  shaped  body.  The  slow  potential's 
contribution  to  the  vertical  drift  force  is  at  leading  order  for  a  regular 
shaped  body,  but  the  resulting  motion  is  0(e2).  For  a  bottle  shaped  body 
where  the  heave  is  large,  the  leading  order  slow  vertical  force  is  given  by 
the  steady  drift  formula  and  the  slow  potential  effect  is  negligible. 

We  now  turn  to  other  applications  of  the  theory. 
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4.  TSAPPING  AND  RESONANCE  OF  LONG  HAVES  ON  A  SHELF 

4.1  Introduction  and  Definition 

It  Is  well  known  that  shallow  water  waves  can  be  trapped  on  a  shelf  or  on 
a  ridge  (or  outside  a  canyon).  These  waves  must  be  propagating  on  the  ridge 
In  a  direction  oblique  to  Its  axis.  It  Is  Impossible  by  any  linear  mechanism 
for  such  waves  to  enter  the  ridge  from  the  deeper  sea;  they  must  be  generated 
by  local  forcing  on  the  ridge.  The  wave  length  and  period  of  the  trapped 
modes  are  much  longer  than  those  of  swell  with  typical  periods  of  0(20  sec), 
but  can  be  close  to  the  modulatlonal  frequency  of  swell  groups,  typically  a 
few  mlntues. 

We  shall  show  that  swell  groups  Incident  from  the  sea  can  excite  trapped 
long  waves  on  a  shelf.  If  the  excitation  Is  resonant,  these  waves  can  be 
greater  than  second  order,  and  cause  large  slow  oscillations.  Hence  the  long 
waves  may  present  a  hazard  to  off~’shore  platforms. 

The  subject  of  Interaction  of  slowly  varying  waves  with  topography,  and 
trapping  on  a  shelf  was  studied  by  Foda  and  Hel  (1981)  and  by  Mel  and 
Benmoussa  (1984)  who  gave  references  to  other  work.  The  latter  authors 
studied  a  slowly  varying  topography  In  which  reflection  by  the  sloping  bottom 
Is  negligible.  In  the  present  work  we  study  the  case  when  reflection  of  the 
short  waves  Is  significant. 

For  demonstration  of  the  physics,  we  choose  a  relatively  simple  bathy¬ 
metry  In  the  form  of  a  rectangular  step.  The  width  of  the  shelf  Is  comparable 
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to  the  wave  group  length,  which  is  much  greater  than  the  short  wave  wave¬ 
length.  The  water  depths  are  comparable  to  the  short  wave  length.  To  solve 
for  the  flow  potential,  we  perform  a  perturbation  expansion  in  terms  of  the 
short  wave  steepness,  and  utilize  multiple  time  and  (horizontal)  space  scales. 

At  the  leading  order,  the  short  waves  are  found  by  the  solutions  of  a 
linear  diffraction  problem  studied  by  Mel  and  Black  (1969)  using  a  variational 
method,  or,  equivalently,  by  Galerkin  method.  For  convenience,  the  problem  Is 
solved  as  a  superposition  of  a  symmetric  and  an  antisymmetric  problem.  In 
contrast  to  the  procedure  employed  for  a  floating  body  (Agnon  and  Mel  1985a), 
the  reflection  by  the  shelf  Is  very  sensitive  to  the  short  wave  wavelength. 

This  Is  because  the  shelf  is  long;  even  small  variations  within  a  narrow  banded 
spectrum  give  rise  to  large  phase  differences  that  affect  reflection.  An 
alternative  and  equivalent  approach  is  to  break  down  the  shelf  into  two  steps, 
matching  their  far  fields  (cf.  Newman  (1965))  while  retaining  the  dependence  of 
the  wave  amplitude  on  the  slow  (stretched)  time. 

After  finding  the  short  wave  field,  we  examine  separately  the  near  field 
within  a  few  short  waves  of  the  depth  steps,  and  the  far  field,  a  few  wave 
groups  away.  We  find  the  forms  of  the  potentials  for  the  slow  waves  in  each 
field,  and  match  their  asymptotic  values  to  get  the  solution.  It  is  found 
that  two  types  of  long  waves  exist.  The  first  one  is  locked  to  the  shortwave 
groups  and  propagates  at  the  group  velocity.  The  other  is  a  shallow  water 
wave  generated  at  the  shelf  edges  and  radiated  away  from  them  in  directions 
different  from  those  of  the  short  waves  and  the  locked  long  waves.  When  the 
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depths  and  angle  of  Incidence  are  such  Chat  the  second  long  wave  Is  Crapped  on 
the  shelf,  resonance  can  occur.  The  resonance  obtained  from  the  present 
approximation  Is  unbounded,  since  there  Is  no  radiation  damping  due  to  com¬ 
plete  trapping  of  the  long  waves.  The  onset  of  resonance,  however,  takes  a 
very  long  time. 

We  present  results  for  the  amplitude  of  the  free  long  waves  due  to  peri¬ 
odic  and  transient  Inputs.  The  short  waves  are  Incident  normally  or  obliquely. 

The  Cartesian  coordinates  (x,  y,  z)  are  the  same  as  those  of  Section 
2.1.  The  bathymetry  Is  given  by; 

depth  ”h  (|x|>L,  -oo<y<<«>) 

'  '  (4.1.1) 

depth  =•  h'  (jxj  <  L,  -<»  <  y  <  ») 

where 

h  >  h*  (4.1.2) 

that  Is,  we  are  studying  a  ridge  that  runs  parallel  to  Che  y  axis.  The  water 
Is  of  Intermediate  depth: 

ah  -  0(1)  ;  ah’  •  0(1)  (4.1.3) 

and  the  waves  are  Incident  from  both  left  and  right  for  the  symmetric  and 
antisymmetric  problem.  The  width  of  the  ridge,  2L,  Is  assumed  to  be  comparable 
to  the  length  of  the  wave  groups: 

kL  -  0(e~^)  (4.1.4) 

The  case  of  a  shelf  ending  with  a  vertical  coast  Is  equivalent  to  a  special 
case  of  the  present  problem  where  Che  Incident  waves  from  the  left  and  the 
right  are  symmetric. 
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The  regions  within  a  few  short  waves  from  the  ridge  edges  are  called  the 
near  fields,  and  the  potential  there  is  denoted  by  tl* 

ij)  =  ♦  k(x  ±  L)  =  0(1)  (4.1.5) 

The  regions  a  few  wave  groups  away  from  the  ridge  edges  are  called  the  far 
fields,  and  the  potential  there  is  denoted  by  (j»: 

(J(  =  *  k(Xj^  -  =  0(1)  ;  k(x^  +  =  0(1)  (4.1.6) 

where 

(xj,  y^,  t^;  Lj)  =  e(x,  y,  t;  L)  (4.1.7) 

are  stretched  quantities  used  to  study  the  wave  groups  and  the  long  waves. 

In  each  region  the  potential  is  expanded  into  orders  and  harmonics,  as  in 
Section  2.2.  Let  us  first  examine  the  short  wave  potentials, 

4.2  The  Short-wave  Potential 

In  the  near  fields,  the  equations  for  are  formally  the  same  as  the 
equations  for  the  diffraction  of  regular  waves  by  a  step.  Two  approaches  can 
be  taken.  One  is  to  spectrally  decompose  the  incident  wave  into  its  short 
wave  components  and  to  superpose  the  respective  solutions  to  get  the  total 
diffraction  wave  field.  It  is  well  known  that  the  diffraction  by  a  long  shelf 
is  very  sensitive  to  small  changes  in  the  wavenumber.  Along  the  shelf,  small 
phase  differences  due  to  dispersion  have  a  large  accumulating  effect  on  the 
Interference  of  waves  reflected  from  the  two  edges,  as  seen  in  the  results  of 
Mel  &  Black  (1969).  An  alternative  way  is  to  use  the  multiple  scale  analysis, 
and  to  follow  the  local  wave  amplitude  of  the  wave  group,  which  carries  the 
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required  phase  information.  Since  the  shelf  is  large,  a  wide  spacing  approxi¬ 
mation  can  be  employed,  taking  into  account  only  the  propagating  modes  for  the 
interaction  between  the  two  edges  (cf.  Newman  (1965)),  as  explained  below. 

Since  the  shelf  is  much  longer  than  the  wave  length  of  the  short  waves, 
we  may  assume  that  only  the  propagating  modes  from  each  edge  reach  the  oppo¬ 
site  edge  and  that  the  evanescent  modes  may  be  neglected. 

Consider  an  infinitely  wide  depth  step  at  (x  =  -L).  Let  the  potential  of 
a  wave  Incident  from  the  left  of  the  step  (from  deep  to  shallow  water)  have  a 
unit  amplitude: 

*^11  ""  ^0  <  “L)  (4.2.1a) 

with  a  ”  k  cos  0;  y  =  k  sin  0 

where  0  is  the  angle  of  Incidence,  then  the  propagating  modes  of  the  reflected 
and  transmitted  waves  are  given  by 

(|)^J  *  ^  ^0  ^  (4.2.1b) 

-  T  fg  expdax  +  iyy)  (x  >  -L)  (4.2.1c) 

respectively,  with  k'^  ™  a'^  +  R'  and  T*  are  used  to  denote  the 
reflection  and  transmission  coefficients  for  a  wave  incident  from  the  shallow 
side  of  the  depth  step.  The  values  of  R,  R',  T,  and  T'  are  determined  using 
the  Galerkln  method  used  in  Section  3.2  for  the  floating  rectangular 
cylinder.  is  replaced  by  {f^}  which  are  the  eigenfunctions  for  the 
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onshelf  potential.  given  by  (2.3.5a,  b)  with  h  replaced  by  h'. 

Simple  adjustments  are  made  for  oblique  incidence  as  in  Black  (1970).  Values 
of  R,  R',  T,  and  T'  are  plotted  in  Figure  4.3.  Energy  conservation  was 
verified  by  checking  the  numerical  results  against  the  identity: 

(1  -  |r|^)  Cg  cos  e  =  |t|  cos  6 
where  tan  9'  “  y/a' 

Now  consider  the  shelf  of  finite  width  2L.  We  denote  the  short  wave 
potential  in  the  far  field  off  the  shelf  by: 

=  (A  +  B 

(4.2.2) 

*11  ”  fQ  (x  >  L) 


We  study  the  case  of  waves  incident  from  the  left  only.  The  case  of  waves 
incident  from  both  sides  and  in  particular  symmetric  incidence,  which 
corresponds  to  a  shelf  of  width  L  terminated  by  a  vertical  coast,  can  be 
treated  in  the  same  way. 

The  propagating  modes  on  the  shelf  are  given  by: 

(D  +  E  (jj^l  <  L)  (4.2.3) 


(see  Figure  4.2).  Matching  the  propagating  modes  in  (4.2.2)  and  (4.2.3)  and 
making  use  of  the  reflection  and  transmission  coefficients  for  each  step,  we 
immediately  get: 


RA  +  T'E  e“^®'^“^^  »  B 

TA  e^®^“^^  +  R'E  e“^®'^“^^  =  D 


(4.2.4) 


at  (x  =  “!.),  and 
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Ridge  Geometry 

Matching  Boundaries  for  a  Ridge 


9 


Figure  4.3  The  Reflection  Coefficients  ((rI  »  )r'|')  and  Transmission 

Coefficients  (T  and  T')  for  a  Single  Step  versus  kh';  h'/h  *0.5 
Angles  of  Incidence  are:  a)  0®(normal),  b)  40°  and  c)  55°. 


(4.2.5) 


T-D  =  C 

R-D  -  E 

at  (x  “  L) 

In  general.  A,  B,  C,  D,  and  E  should  be  treated  as  functions  of 
scale.  However,  for  a  sinusoidal  envelope  It  Is  easier  to  decompose 
dent  wave  as  follows: 


long  time 
the  Incl- 


^  ^Kax+yy-oit)  ^  ^Ko^x+Y^y-m^t) 


+  gl(o  x+Y  y-tu  t) 


(4.2.6) 


where  A"*"  =  A“  Is  a  real  constant. 


+  _*  -in(tj  -  yx  -  vy  ) 

A  =■  Re  {(A  +  A  )  e  } 


kC  ’ 


kC 


g 


and 


+  —  +  _ 
a  =*a+en;a  =  a“  en»  u  =  u+eu  ;  u)“  =  u  •"  eft  5  eft/w 

The  rest  of  the  amplitude  functions  are  written; 


(4.2.7) 

(4.2.8) 

0(e)  (4.2.9) 


B  -  Re  {(B  +  B  )  e 


_*  -lft(tj  +  yXj  -  vy^) 


} 


+  _*  -in(t,  -  yx  -  vy,) 

C  =  Re  {(C  +  C  )  e  } 

+  _*  -ln(ti  -  y’x  -  vy^) 

D  -  Re  {(D  +  D  )  e  } 

.  _*  -lft(tj  +  y'x,  -  vy,) 

E  -  Re  {(E  +  E  )  e  } 


(4.2.10) 


where 
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(4.2.11) 


or,  equivalently  In  the  form  of  (4.2.6). 

+  +  +  +  + 

Substituting  A~,  B  ,  C  ,  D“,  and  E"  in  place  of  A,  B,C,  D,  and  E  in 
(4.2.4)  and  (4.2.5)  results  in  a  matrix  equation: 


la*L 

e 

0 

0 

-T'e^“'  ^ 

B- 

0 

0 

-la'L 

e 

-R'e^“'  ^ 

• 

C* 

TA*e-^“^^ 

0 

la~l* 

-e 

0 

D- 

0 

0 

0 

-la^'b 

“e 

E- 

0 

which  can  be  Inverted  to  give: 


•  m 

R'  exp(-lo^)  +  R’TT'  exp[l(4o"'-o")L] 

c- 

A“ 

T'T  exp(2l(o~’-a')L] 

D“ 

1-R’^exp(4la*  L) 

T  exp[l(a“  “a*)b] 

E* 

•  m 

R'T  expli(3a“’-a*)L] 

(4.2.12) 


cf •  Newman  (1965). 

Values  of  B-/A-  are  plotted  in  Figure  4.4.  Since  the  small 
differences  between  a"*"  and  a”  have  an  0(1)  effect  on  the  phase  for  waves 
travelling  across  the  whole  shelf,  the  amplitudes  for  the  two  components,  ()-, 
differ.  The  solution  given  by  (4.2,12)  could  be  found  without  using  the  far 
spacing  approximation.  An  exact  solution  for  B  and  C  was  given  by  Mel  and  Black 
(1969)  using  a  variational  method.  Their  results  for  h'/h  ■  0,5  and  L/h  “  3 
given  in  their  Figure  2  are  in  excellent  agreement  with  the  present  results 
where  we  assume  a  long  shelf  -  in  the  range  of  frequencies  corresponding 
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0.4 


Figure  4.4  The  Reflection  Coefficient  |b-/A±|  for  a  Ridge  versus  kh'.  h'/h 
0.5;  a)  L/h  =  10,  b)  L/h  =  3.  Results  of  Mei  &  Black  (1969)  are 
Shown  in  a  Dotted  Line  for  Case  b. 
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to  (4.1.3))  as  shown  In  Figure  4.4.  For  a  transient  problemA,  B,  C,  D,  and  E 
are  more  general  functions  of  Xj,  y^,  tj.  Results  of  using  (4.25)  (which  relies 
on  the  multiple  scale  approach)  for  transient  envelopes  will  be  presented 
later.  For  a  sinusoidal  envelope  we  can  substitute  (4.2.8)  and  (4.2.10)  Into 
(4.2.3).  We  end  up  with  the  same  solutions  for  B"*",  B“,  C"*",  C~,... 

An  alternative  to  (4.2.12)  Is  to  decompose  the  short  wave  potential  at, 
say,  m'*’.  Into  a  component  that  Is  symmetric  with  respect  to  (x  “0)  and  an 
antisymmetric  component.  This  results  In  four  sets  of  equations,  for  two 
unknowns  each.  The  results  are,  of  course,  the  same. 

This  completes  the  determination  of  the  propagating  short  waves.  We  may 
turn  to  the  long-wave  potential. 

4.3  The  Long-wave  Potential 


The  governing  wave  equation  for  the  slow  potential  in  the  far  field  is 
(2.4.18).  We  are  still  considering  the  monochromatic  envelope  Input  (4.2.6) 
and  wish  to  study  the  oscillatory  motion  that  results  from  It.  There  are  two 
types  of  long  waves.  One  Is  locked  to  the  short  wave  groups,  and  Is  similar 
to  the  waves  given  by  (3.1.9).  The  only  difference  is  in  the  direction  in 
which  these  waves  propagate.  Off  the  shelf,  where  thr  short  waves  are  unidi¬ 
rectional,  the  locked  long  waves  are  given  by  the  solution  to  (2.4.20). 


Re  ^e 


-2in(tj  ±  u  Xj  -  vyj) 


2  2 


gh  -  c. 


I  (k^  -  o^  +  (4.3.1) 


where 
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X  *  A,B,C 


(4.3.2) 


A,B,C,D,  and  E  refer  to  the  four  wave  trains  given  by  (4.2.7)  and  (4.2.10); 
(-)  Is  taken  for  A  and  C;  (+)  Is  taken  for  B. 

On  the  shelf,  the  dispersion  makes  the  wave  envelopes  propagate  at  a 
different  direction  from  the  short  waves.  From  (2.4.18)  we  get: 


2 


gh*  -  c- 


(k*^  -  +  2(0  (a'  u'  +  YV))  (4.3.3) 


where 


(-)  is  taken  for  D  and  (+)  for  E. 


X  =  D.E 


(4.3.4) 


There  will  also  be  free  waves,  propagating  away  from  the  edges.  The  free 


waves  are  solutions  of  the  D'Alembert  equation: 

♦lOtjtj  ■  ’Mo  ■  ° 
where  h  is  replaced  by  h'  on  the  shelf. 


(4.3.5) 


Because  the  potential  is  periodic  in  y^  for  the  periodic  envelope 
(4.2.6),  the  free  waves,  too,  must  be  periodic  with  the  same  wavenumber  as 
jAj^,  namely  2nv. 

The  Xj  (cross-ridge)  wavenumbers  of  the  free  waves  2in8,  2jib',  must 


satisfy: 


2  2  -1 
+  V  -  (gh) 

2  2  -1 
g'^  +  -  (gh') 


(4.3.6) 


(4.3.7) 
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We  shall 


3,  3*  will  be  real  or  imaginary,  depending  on  the  signs  of  3^,  3'^. 
study  the  interesting  case  where  the  free  long  waves  are  trapped  as  propagat¬ 
ing  modes  on  the  ridge,  and  are  evanescent  off  the  ridge.  In  this  case  reso¬ 
nance  may  occur.  The  free  waves  will  have  the  form: 

-2ift(tj  -  V  yj)  +  2n3xj^ 


♦lO  ■  '*’10 
with 


♦lo  “  *10  ■  ^  ® 


10 

H 

♦lO 

♦lO 


-2lfi(t,  -  3'x,  -  V  y, ) 


Re  H  e 


Re  I  e 


-2ln(tj  +  3'xj  -  V  y^) 


(xj  <  -Lj) 


(|xi|  <  Lj) 


(4.3.8) 


and 


♦lO  ” 


Re  J  e 


-2ln(tj  -  V  yj)  -  2n3Xj 


(xj  >  Lj) 


with 


o2  1  2  .  „  .  „,2  1  2 

-6  -  -  V  <  0  <  6-  ■  Jip-  -  V 


(4.3.9) 


G,  H,  and  J[  will  be  determined  by  matching  to  the  near  field  potential. 
The  total  far  field  is: 


A  B  G 

♦lO  “  *10  *10  *10 


(xj  <  -Lj) 


*10  "  *10  *10  *10  *10  ^|*l|  ^  ^1^  (4.3.10) 

”  ♦n  ♦in  ^*i  ^  ) 


from  (4.3.1)  and  (4.3.8) 
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In  the  near  field,  the  geometry  Is  the  same  as  in  the  wide  gap  case 
Section  3.3.1.  Due  to  the  arguements  of  Section  3.3.1  (Equation  3.3.16),  the 
slow  potential  i|>jq  is  Independent  of  short  scale  (x,  y,  z).  It  does,  however, 
depend  on  y^  due  to  variation  along  the  shelf: 

"^10  “  '^'lO  (4.3.11) 

where  tj;^Q  is  the  potential  at  the  right  edge  and  is  the  potential  at  the 
left  edge  of  the  shelf. 

Just  as  in  Section  3.3.1,  ij/^g  does  not  depend  on  and  the  second  order 
slow  drift  velocity  in  the  x  direction  is  accounted  for  by  i(»20x*  Integrat¬ 
ing  the  normal  derivative  of  i|)2o  around  the  'boundaries'  of  the  near  field 
S'*",  S“  ,  S„  and  S„  (see  Figure  4.1)  at  say  (x,  “  -L,  =  -  eL),  we  find  as  in 

Section  3.3.1  a  result  similar  to  (3.3.34).  It  is: 

^  ’^20x  /  ’*’20x  "  -[h'U(-L+o.)  -  h  U(-L-»)]  (4.3.12a) 

S‘*'(-L)  S‘(-L) 

00  00 

where  S-(-L)  and  (x  +  -  L  ±  «)  both  stand  for  the  matching  boundaries  to  the 
00 

right  and  to  the  left  of  the  edge  at  (x^  =  -Lj).  The  derivation  of  (4.3.12) 
is  even  simpler  than  that  of  (3.3.34)  in  the  floating  body  case,  since  there 
are  no  moving  boundaries.  At  (xj^  *  L^),  we  have,  similarly, 

/  iP20x  “  -[hU(L+«.)  -h'(U(L-o.)]  (4.3.12b) 

sV)  S“(L) 

00  00 

To  match  with  the  near  field  near  the  depth  discontinuity  at  (x  =  -!>),  we 
may  write  the  inner  expansion  of  (see  3.1.11)  as  follows: 
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A  B  Q  A  B  G 

♦lO  t*10  *10  ♦lol  .  *10  ^  ♦lOx) 

Xj  -  -L  1  1  _ 

(xj  -  L^) 

'^’lO  '*’  *10  '''  *10  *10^  r  *10  *10  *10  *10^ 


*1  "  "^1 


(Xj  ■►  -  Lj^) 


*10  ‘  t*?o  *10  *10  ^  *io]  ^  *lt"‘'’  *10  *10  *10  ■  *io] 


xi  -  Lj 


-  r  C  J  1  .  r  C 

*10  l•*lo  *ioJ  ,  *iot,  ■  nox, 

'l  ■  ‘■1  ‘  ‘ 


J  1 

*10,.  ] 


(Xj  +  Lj) 


(xj  ♦  Lj)  (4.3.13) 


Use  has  been  made  of  the  simple  relations  between  t^  and  x^  derivatives  of 
these  terms. 

Matching  (4.3.13)  to  the  outer  expansion  of  the  slow  potential  from 
(4.3.11)  and  (4.3.12)  we  get: 

*10  *10  *10  “  *10  “  *10  *  *10  *10  *10  ^*1  "  "^1^ 


D  E  H  I  +  C  J 

*10  *  *10  *10  *10  “  *10  “  *10  *  *10 


(xi  -  Lj) 


^  *iotj  *iotj  *i0xj  ■  *10  *10  "  *10  *  *io]j 


(xi  -  -L^) 


-  h'  U(-L  +  «)  -  h  U(-L  -  ») 
h'  [-w’  tio  *10  “  *10  *lo]  ■  *10t,  ~  *10x  ] 


h  U(L  +  »)  -  h'  U(L  -  «) 


(xi  -  Lp 


(4.3.14) 
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On  the  other  hand,  the  Stokes  drift  termb  nil  and  h'U  are  determined  from 


(3.3.27)  as  in  (3.3.33).  They  are: 


u(i) 

-2ifitj 

e 

_  2(1)0 

g 

u(2) 

-2intj 

e 

2(1)0' 

g 

„(3) 

-2lfitj 

e 

2(1)0 ' 

g 

u(4) 

-2l0tj 

e 

2(1)0 

fj(0)  -  |i|^) 

■f>)  C|D|^-  |E|^) 
■f>)  (H^-  |E|^) 


<1.1  .  -Lj) 


(Xj  .  -Li) 


(Xi  =  I-j) 


(x,  -  Lj) 


(4.3.15) 


where  j  j  stands  for  the  oscillatory  part  of  |  j  ^  and  has  an  exponential 
factor  as  in  (4.3.1). 

We  note  that  in  contrast  to  the  cases  of  Section  3.3,  where  the  Stokes' 
Drift  terms  eventually  cancelled,  they  have  a  non  zero  net  contribution  in  the 
present  problem.  The  reason  is  that  in  the  floating  body  problem,  conserva¬ 
tion  of  energy  was  equivalent  to  conservation  of  mass  in  the  drift  terms, 
since  the  depth  and  wavenumber  were  the  same  on  both  sides  of  the  body.  At 
the  present  case  the  depth  difference  results  in  a  net  mass  flux. 

When  (4.3.1),  (4.3.3)  and  (4.3.8)  are  substituted  into  (4.3,14)  we  arrive 
at  four  linear  equations  for  £,  H,  and  These  can  be  solved  to  yield  a 
full  solution  for  the  long  waves,  without  having  to  explicitly  solve  for  the 
full  second  order  potential. 

We  get: 
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where  the  matrix  b  Is  given  by: 


“2fi3L2 

e 

-2lnB'Lj 

2lft3'Lj 

-e 

-e 

0 

0 

2lfi6'Lj 

e 

-2in3*Lj 

e 

-2n3Lj 

-e 

-2nBL. 

2nh0e 

-2lfiB'L, 

-Zlfih'B'e 

21£JB*L 

Zlnh'B'e 

0 

(4.3.17) 

0 

2lflB'L 

-2ltih'8'e 

-2ln3’L 

Zlfth'B'e 

-2nBL 

-2fthBe 

and  ^^9  ^2*  ^3* 

and  a^  are  the  locked 

wave  components  In 

(4.3.14): 

-ZlftyLj^  2iftyLj  “2lSlyLj  2iftyLj 

a,“-Ae  “Be  +De  +Ee 


2l{jyL.  “21nyL.  21nuh. 

“C©  "E©  “D© 

-2lnyL.  2lnyL. 

aj  -  -h[(-yA  e  +  y B  e  )(-2ln)  +  u''*''] 

-Zlny'L  21ny'L  f  . 

+  h'  [(-y'De  +  y'i  e  )(-2ln)  +  ] 

-2lny’L.  Zlny'Lj 

a^  -  -h  [(-y'E  e  +  y'D  e  )  (-2lfi)  + 

2lnyL 


+  h'[“  y'C  e 


^  (-2ln)  + 


The  matrix  (4.3.17)  becomes  singular  when: 


4lng'L 

Im((Bh  +  IB'h')^  e  =  0  (4.3.18) 

(see  Mei  (1983,  p.  141))  which  means  that  there  is  infinite  resonance  of  the 
trapped  free  waves  at  the  frequencies 

2n^B’Li  =  tg"^  Y  nil  (4.3.19) 

We  note  that  the  second  and  third  columns  of  (4.3.17)  are  complex  conju¬ 
gates  of  each  other.  This  means  that 

i(2fi)  =  -H*(2fi)  (4.3.20) 

Hence  the  free  wave  on  the  shelf  is  of  the  form 

2lft(t^  -  B'xj  -  vy^)  *  2ln(tj^  +  3'  Xj^  -  vy.) 

Re[H  e  -He  ] 

2in(tj-vyp 

=  RetjHj  e  ]  sly(2ji3'x  -  argH) 

which  is  a  standing  wave  in  the  x^  direction,  propagating  in  the  y^  direction. 
II  is  found  from  (4.3.16)  by  Cramer's  method; 


(4.3.21) 
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where  the  left  hand  side  of  (4.3.18)  was  written  as  a  sine. 

Figure  4.5  shows  the  variation  of  the  normalized  h|llj/A^  with  n/u.  The 
resonant  frequencies  correspond  to  the  roots  of  (4.3.19).  The  spacing  between 
them  is  simply  ir/(4g'Lj). 

The  location  of  the  first  resonance  Oq  is  given  by 

-  FT7 

The  leading  behaviour  of  Ug  is  due  to  the  factor  l/(B'Lj).  In  Figure  4.6,  fig 
is  plotted  versus  the  angle  of  incidence  of  the  waves. 

In  the  special  case  of  normal  incidence  (y  =  v  “  0),  Equation  (4.3.8)  is 
replaced  by: 


p  -2lfi(t,  +  x.//gh) 

♦lO  “  *10  “  ®  ^*1  ^ 

“^lO  "  *10  *10  ^|*l|  ^ 

with 

-2ifi(t,  -  x.//gh') 

♦lO  -  ReH  e  (4.3.23) 


j  -2ifi(tj  -  x^//gh’) 

♦  lO  “  ® 

and 


y  T  -2ifi(t.  -  X  //gh) 

*10  "  *10  "  ®  ^*1  ^  "^1^ 

The  free  wave  off  the  shelf  is  now  propagating.  Equations  (4.3.16)  and 

(4.3.17)  are  still  valid  (with  6  imaginary).  The  normalized  h|H|/A^  is 

plotted  in  Figure  4.5.  No  infinite  resonance  occurs,  and  the  variation  with  fi 

or  L  is  minute. 
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h/h'=  .25  .5  .75 


109 


Anplc  o'  incidence  (degrees' 

Figure  4.6  The  Modulation  Frequency  Hq  at  Which  the  First  Pole  Occurs  and  the 
Value  of  Hg,  the  Residue  at  fig  for  Three  Depth  Ratios: 
a)  h'/h  =  0.25,  b)  h*/h  =  0.5  and  c)  h'/h  =  0,75; 

(L/h  =  10.,  ah  =  4.). 


4.4  Transient  Response 


We  have  found  that  the  trapped  long  waves  due  to  obliquely  Incident  short 
waves  can  have  an  Infinite  resonance  at  certain  modulation  frequencies.  This, 
of  course.  Is  physically  impossible.  Note,  however,  that  the  present  theory 
relates  to  times  on  the  scale  of  t^,  whereas  the  frequency  response  Is  calcu¬ 
lated  for  a  periodic  Incident  wave  envelope.  Longer  time  scales,  and  hence, 

higher  order  effects  need  to  be  considered.  Foda  &  Mel  (1981)  have  derived  a 

•• 

nonlinear  Schrodinger  Equation  for  the  development  of  the  resonated  long  waves 
on  a  mildly  sloping  ridge.  They  found  that  a  higher  order  radiation  damping 
limits  the  magnitude  of  resonance.  We  shall  examine  instead  the  transient 
response  of  the  ridge  to  a  narrow  banded  incident  wave  of  finite  duration  on 
the  time  scale  of  tp  In  particular,  for  oblique  incidence  we  shall  focus  on 
the  case  where  the  spectrum  peak  coincides  with  a  resonant  frequency. 

We  start  by  looking  at  the  case  of  normal  incidence,  which  does  not  lead 
to  wave  trapping,  but  Is  the  easiest  to  perform  experiments  for. 

When  Is  positive,  the  governing  equation  for  can  be  easily  solved 
using  finite  differences  in  time.  The  response  is  passive,  and  the  free  waves 
and  the  forced  waves  'echo'  on  both  edges  of  the  shelf. 

The  method  of  solution  is  as  follows:  first  we  solve  for  the  short 
waves.  Given  A(Lj,tj)  we  note  that: 

D(Lj,tj)  »  D  (-Lj,  tj  -Li/C^);  E(-Lj,  tj)  =  E(L^,  t^  "Lj/Cp 

making  use  of  these  relations  to  account  for  the  propagation  across  the  ridge, 
we  solve  (4.2.4)  and  (4.2.5)  at  each  time  step,  to  find  the  short-wave 
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amplitudes:  A,  B,  D  and  E  at  (xj^  L^)  and  C,  D  and  E  at  (x^  *■  L^),  see 
Figures  4.7a.  The  locked  long-waves  are  computed  from  these  amplitudes 
through  (4.3.1)  and  (4.3.3)  and  the  Stokes*  drift  terms  are  found  from 
(4.3.15). 

We  can  now  eliminate  the  Xj  derivative  from  (4.3.14)  via: 

G  1  I  G  J  I  .  I  d 

*10  "  "|®|  *10  ’  *10  "  |®|  *10^ 

*1  ^1  *1  ^1 

to  get  coupled  equations  for  (o  “  G,  H,  I,  J)  that  involve  only  tj 

derivatives.  Making  use  of 


*10^^!’ 


H  /  T 
*10^"^1’ 


ti  -L^/ 


✓gh ' ) ; 


J  /T 
*10^^1 


.tl) 


J  /  , 
*10^"^1» 


/gh*  ) 


We  Integrate  (4.3.14)  numerically  with  respect  to  time  from  (tj*0)  on  and 
determine  the  free  long  waves  at  -Lj  and  at  L^. 

In  Figure  4.7b  the  wave  characteristics  of  the  long  waves  are  sketched  In 
an  (x^,t^)  diagram. 

The  resulting  transient  response  to  a  slnosoldal  packet  at  the  left  edge 

(xj  -  -Lj)  and  at  the  right  edge  (xj  -  Lj)  is  plotted  in  Figure  4.8. 

Interference  of  the  free  waves  and  the  locked  waves  can  be  seen  after  one 

reflection.  At  later  times  the  response  Is  negligible. 

F 

For  oblique  Incidence  with  wave  trapping,  becomes  evanescent  off  the 
shelf  and  we  shall  use  a  Fourier  Transform  of  the  Incident  wave  In  conjunction 
with  the  frequency  response  to  determine  the  transient  response.  Let  the 
Incident  wave  envelope  be  a  Gaussian: 
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(a)  The  Short-wave 
Amplitudes. 


Groups  and  their 


Figure  4.7 
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o 

o 


Figure  4.^  The  Long-wave  Transient  Response  of  a  Ridge  to  Normally  Incident  Waves 

a)  at  the  Left  Edge  and  b)  at  the  Right  Edge  of  the  Shelf, 
versus  nt,  Ch’/h  =  0.5,  L/h  »  10,  oh  =  1,  n/u  =1)  1)  the  locked 
Waves.  2)  the  Free  Waves,  and  3)  the  Total  Slow-wave. 
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(4.4.1) 


,  1  2  2 
“in«p  -  w  p 

A  =•  Re  e 

where  p  Is  the  phase  function: 

p  =  ti  -  Pi  Xi  -  V  yi 

and 

w  «  Oq 

Is  the  characteristic  length  of  the  packets. 

The  Fourier  Transform  of  |a|^  with  respect  to  2j2  is 


^2inp 


„  ^  ,1  ^  I  2illoO-»%^  1  -21!1oP-w%^ 


2,2 


.2,  2 


2,  2  /w  -(fi+JV,)  /w 

[2e‘“  ^  ♦  P  ] 


Sw/ir 


(4.4.2) 

(4.4.3) 


(4.4.4) 


Let  us  examine  the  transient  response  for  H,  the  right  going  free  wave  on 
the  shelf  at  a  point  along  the  axis  (xj  =  0,  yj  »  0). 

“  I  1 2 

H(0,  0,  tj)  -  2Re  /  d2j2  H(2n)  e  F.T.jAp)  (4.4.5) 

—00 

since  II  (2n)  is  the  complex  amplitude  of  the  frequency  response  to  a  component 
of  |Aj^  in  the  form 

jRee"^^”P  (4.4.6) 
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To  accentuate  the  resonance  effect,  we  choose  2ng  that  lies  on  a  pole  of 
Il(2jj).  The  leading  term  of  (4,4.5)  will  be  the  contribution  from  that  pole. 
The  Gaussian  dependence  will  make  all  other  contributions  negligible  because 
of  (4.4.3),  and  only  the  second  term  of  (4.4.4)  should  be  kept.  The  result  Is 

,  «  -2int,  -  (jl-iU)^/w^ 

H(0,  0,  tj)  =  — ^Re  /  d  (4.4.7) 

where  Hg  Is  the  residue  of  11  at  the  pole  of  H  as  a  function  of  2(i  at  2(2 
2(2g.  The  values  of  Hg  will  be  discussed  later. 

The  Integration  path  Is  slightly  above  the  real  axis  In  order  to  make  H 
tend  to  zero  as  tj  goes  to  -<«».  This  Integral  can  be  analyzed  by  the  method  of 
steepest  descent.  In  the  complex  plane  2(2,  the  phase  function 

g(2(2)  -  -2i(2tj - 2 —  (4.4.8) 

w 

has  a  saddle  point  s  at: 

2(2  “  2(2q  -  2iw^tj  (4.4.9) 


(see  Figure  4.9).  We  transform  the  Integration  path  to  a  path  through  s: 


H(0,  0,  t,  )  »  Re  / 


■2iw^t 


1  Hq  -2i(2tj  -  ((2-(2q)^/w^ 

d2(2  ---■  e 


4w/x 


-00-21W  t, 


Hq  ~2i(^t. 

+  U(t,  )  Re[-2Trl  — ^  e  ] 

4w/x 


(4.4.10) 


where  U  Is  the  unit  step  function  and  the  last  term  Is  due  to  the  residue  and 
Is  added  by  Cauchy's  theorem  when  the  pole  at  2(2g  Is  traversed  for  positive 
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t^.  The  contributions  from  the  vertical  segments  connecting  the  present  path 
with  the  real  axis  at  both  infinities  vanish  because  of  the  Gaussian  factor. 

The  Integral  in  (4.4.10)  can  be  evaluated  by  the  following  change  of 
variable: 

a  *  2n  -  +  2lw^  tj  (4.4.11) 

It  becomes 


4wi/ir 

This  last  integral  has  been  evaluated  (Carrier  1970)  and  is  given  by 

2  2 
w  ^ 

sgn(tj^)  in  e  ^  [1  -  erf  wjtjj]  (4,4.13) 

(4.4,10),  (4,4,12)  and  (4,4,13)  can  be  combined  to  give: 


-w^tj  -  2tnQtj 


/ 


2/,  2 
-a  /4w 


a  -  21w  t, 


(4.4.12) 


H(0,  0,  tj)  a  Re — {sgn(tj)  (1  -  erf  w|tj|]  -  2  U(tj)}  e  (4.4,14) 


Equation  (4.3.21)  describes  the  spatial  structure  of  the  total  free  wave  on 
the  shelf . 

The  strength  of  the  transient  response  (4.4.14)  depends  on  the  energy  in 
the  input  packet  via  1/w  which  represents  its  length,  and  on  Hq. 

To  find  Hq,  we  replace  the  denominator  of  (4,3.21)  by  23'Lp  hence  Hg  is 
proportional  to  1/(3'L^)  and  is  an  almost  linear  functions  of  Oq  cf. 

(4.3.22).  In  Figure  4,6  Hg  is  plotted  versus  the  angle  of  incidence  of  the 
short  waves . 
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Figure  4.10  shows  the  response  of  the  trapped  waves  to  the  transient 
Input  (4.4.1).  As  seen  In  (4.4.14)  there  Is  a  transient  part  that  builds  up 
before  the  arrival  of  the  peak  of  the  envelope  and  decays  after  Its  arrlv'il, 
and  a  steady  ’wake’  due  to  the  contribution  from  the  pole,  that  follows  the 
short-wave  groups  passage  along  the  ridge  and  keeps  oscillating  at  the 
resonance  frequency  on  the  ridge,  long  after  all  other  waves  are  gone.  The 
steady  reverberation  In  the  wake  Is,  of  course  valid  only  for  time  scales  that 
are  within  the  scope  of  the  theory,  namely  on  the  t^  scale.  Over  much  longer 
time  scales,  higher  order  effects  will  diminish  the  wake.  This  behaviour  Is 
typical  of  various  wave  guides. 

We  now  turn  to  the  diffraction  of  slowly  varying  waves  In  a  two  layer 
fluid. 
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Figure  4.10 


O 

O 

• 

o 

UfSJ 


122 


5.  EXCITATION  OF  INTERNAL  HAVES  IN  THE  LEE  OF  A  BREAKHATER  BT  SLOHLT 

VARYING  SURFACE  HAVES 

5.1  Introduction  and  Definitions 

It  Is  well  known  that  surface  waves  can  Interact  nonlinear ly  to  produce 
Internal  waves.  Past  research  In  this  field  has  been  concentrated  on  the 
resonant  Interaction  of  two  surface  wave  components  with  an  Internal  wave, 
in  an  open  sea.  This  type  of  work  was  started  by  Ball  (1964)  and  is  still  being 
continued  (c.f.  Dysthe  and  Das  (1981)  and  references  therein).  In  this  chapter 
we  shall  study  the  Interaction  away  from  resonance  combined  with  diffraction  by 
a  long  breakwater.  The  short-wave  groups  and  the  "locked"  internal  waves 
propagating  with  them  will  be  shown  to  exist  in  the  lit  zone  but  not  in  the 
shadow  zone.  Free  internal  waves  will  be  required  to  maintain  continuity  of  the 
long  wave  field  across  the  shadow  boundary.  We  shall  see  that  the  free  long 
waves  are  slower  than  the  locked  long  waves  and  hence  the  free  long  waves  will 
radiate  away  from  the  shadow  boundary  into  the  shadow  zone  (and  the  lit  zone). 
This  mechanism  can  be  Important  for  resonant  excitation  in  bays  and  behind 
straits.  The  method  of  solution  will  again  Involve  the  use  of  multiple  scales 
and  asymptotic  matching  of  near  fields  and  far  fields. 

The  sea  is  modeled  as  a  two  layer  fluid  with  the  lower  layer  of  density  p 
and  depth  h,  and  the  upper  layer  of  density  p'  and  depth  h'.  The  origin  of 
the  coordinate  system  is  chosen  at  the  interface  between  the  two  layers  with 
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the  z  axis  pointing  upwards.  The  breakwater  is  lying  on  the  positive  y  axis 
and  the  waves  are  normally  incident  from  (x  -►  -«)  (see  Figure  5.1). 

Let  us  first  determine  the  forced  internal  wave  in  an  open  sea.  As  in 
Lamb  (1932,  p.  372)  the  short  waves  have  a  potential 

lk.5t  ^  lev 

(^11  =  AgfQ  e^  s  Aq  cosh  k(z+h)  e  (-h  <  z  <  0)  (5.1.1) 

in  the  lower  layer,  and 

lltv  ^  le  V 

(frjl  »  AqI^  e  h  Aq  (C  cosh  kz  +  B  slnh  kz)e^^''  (0  <  z  <  h’)  (5.1.2) 


in  the  upper  layer. 

Matching  the  vertical  velocities  at  (z=0)  we  get: 

-k  sinh  kh  =  -kB  =  lub  (5.1.3) 

where  Agb  is  the  Interface  displacement  amplitude  for  the  short  wave. 
Matching  the  pressures  at  the  interface  gives: 

pdoj  cosh  kh  -  gb)  =  p'dojC  -  gb)  (5.1.4) 

At  the  free  surface  (z  “  h')  the  usual  boundary  condition  leads  to 
2 

u  (C  cosh  kh'  +  B  slnh  kh')  gk(C  slnh  kh'  +  B  cosh  kh')  (5.1.5) 
b,  B  and  C  may  be  determined  from  (5.1.3)  and  (5.1.4): 
b  ”  Ik/u)  slnh  kh 


B  “  slnh  kh 

C  “  p/p'  cosh  kh  -  gk/a)2  (p-p')/p'  sinh  kh 

Substituting  these  in  (5.1.5)  yields  the  dispersion  relation 

4  2 

0)  (p  coth  kh  coth  kh'  +  p')  -  u  p(coth  kh'  +  coth  kh)  gk 

+  (p  -  P’)  gV  -  0 


(5.1.6) 


(5.1.7) 
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far  field  . 


TOP  VIEW 


SIDE  VIEW 


Figure  5-.1  Definition  Sketch. 
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The  surface  wave  mode 


For  a  fixed  k  this  Is  a  quadratic  equation  for 
corresponds  to  the  larger  frequency  of  the  two  positive  roots  of  (5.1.6). 

The  Internal  wave  mode  corresponds  to  the  smaller  uj.  When  one  mode  Is 

t 

selecred.  Equation  (5.1.7)  can  be  substituted  In  (5.1.2)  to  define  fp.  For 
reference  numerical  solutions  for  oi/k/  /g(h  +  h'  versus  k(h  +  h')  are  shown  in 
Figure  5.2  for  both  surface  and  internal  modes,  with  (p-p')/p  “  0.01; 
hV(h  +  h')  =  0.05,  .01  and  0.2. 
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Figure  5.2  The  Normalized  Phase  Velocity,  u)/k//g(h+h ' ) ,  for  Surface  and 
Internal  Waves  versus  k(h+h’):  p'/p  =•  0.99;  h'/(h+h')  is  a) 


=  p[g  C20  *iotj "  r 

Equations  (5.2.1)  -  (5. 2. A)  can  be  solved  for  ^jg,  ^20  ^20’ 

t 

Let  us  assume  that  ^  monochromatic  envelope: 

X.  . 

-  [Re  a  exp  -in  “  ^)]  e  *  (5.2.5) 

g 

where  Cg  Is  the  group  velocity  of  the  short  surface  wave  mode. 

The  slow  potential  is  forced  by  the  square  of  the  short  wave.  It  will  have 

a  component  that  is  time  harmonic  and  proportional  to  exp(-2i£itj)  and  a  steady 

component.  From  here  on  we  shall  study  only  the  first  one.  We  proceed 

•  —  _  » 

to  determine  the  amplitudes  of  the  harmonic  components  of  ^^g,  ♦lo*  C20 
520*  •^he  locked  long  waves  that  propagate  at  the  velocities  ±Cg,  the 

homogeneous  part  of  Equations  (5.2.1)  -  (5.2.4)  has  the  following  coefficient 


^  +  *)] 


(5.2.4) 

(z=0) 


matrix: 


The  determinant  of  (5.2.6)  is 


(5.2.6) 
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0  and  kh'  +0.  This  means  that  the  short-wave  groups  will  resonate  the  long 
Internal  waves  when  they  propagate  at  the  same  velocity,  as  Is  well  known. 
Equation  (5.2.7)  can  be  written  as 

4 

\  p  (C^  -  cj)(c^  -  c:^)  (5.2.8) 

c  K  A  S  i 

g 

where  and  are  the  phase  velocities  of  the  long  surface  and  Internal 

2 

modes  respectlvley .  This  factor  is  similar  to  the  factor  (gh  -  Cg)  In 
Section  2. 

Equation  (5.2.8)  is  equivalent  to  (5.1.6)  in  the  limit  of  shallow  water 
waves.  and  are  given  by: 

cj  -  g/2  {h+h*  -  [(h-h*)^  +  4  (p'/p)h  h*]^^^}  (5.2.9) 

-  g/2  {h+h'  +  ((h-h')^  +  4  (p’/p)h  h']^^^}  (5.2.10) 

The  long-waves  solutions  of  (5.2,1)  -  (5.2.4)  will  Include  locked  and  free 
long-waves.  The  locked  long  waves  are  found  by  solving  the  matrix  equation: 

[b] 

where  (a^,  aj,  Sj,  a^}  are  the  forcing  terms  in  (5.2.1)  -  (5.2.4): 


»  m 

♦lO 

*1 

L' 

♦lO 

®2 

-L 

^20 

*3 

-L’ 

^20 

«4 

IS 
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^2*  •"  *^x^"  ■"  *^x^ 

The  values  of  a^,  a.2,  and  a,^  are  given  by: 


(z  =  0) 


(z  -  0) 
(z  =  h') 


k^Uo  ^ 


-2ln 


wgC, 


slnh  2k.h 


[(B^  +  C^)  slnh  2kh'  +  4BC  slnh^  kh’] 


ojgC^ 


p  -  p'  (C^  -  B^) 


2  2 
C  -  B 


(5.2.13) 


**  2  2 

where  jAg|  Is  the  harmonic  part  of  Aq|  .  It  Is  given  (from  (5.2.5))  by: 


A  ^  «  Jl  _ 8^  I  ”■  Re  e 

^0  2  2la)f;(h')  '  ® 


-2l«(t^  --) 


(5.2.14) 


as  In  (3.3.13). 

We  denote  the  complex  amplitudes  of  and  by  E  and  E',  respectively: 
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-  Re  E  exp  -2in(tj  -x^/C^)  ;  <j»jQ  -  Re  E'  exp  -2i{2(tj  -Xj/C^)  (5.2.15) 

—  * 

and  those  of  ?2  0  ^20  ^  respectively: 

^20  =  Re  Z  exp  -2in(tj  -Xj/C^)  ;  ^20  ”  Z'  exp  -2lfi(tj  -Xj/C^)  (5.2.16) 

After  eliminating  the  factor  exp  -2lft(tj  -  x^/Cg),  the  coefficients  E,  E', 
and  Z'  are  found  from  the  matrix  equation  (5.2.11)  numerically. 

Figure  5.3  shows  the  locked  long  waves  amplitudes  Z  and  Z'  versus 
k(h  +  h’).  The  parameters  are  h/(h  +h')  “  0.4,  0.5,  0.6  and  (p  -  p')/p  =  0.1. 
The  resonance  at  the  roots  of  (5.2.8)  is  apparent. 

The  free  long-wave  will  be  determined  in  Section  5.4. 

5.3  •  The  Short-wave  Diffraction 


For  simplicity  we  shall  only  consider  normal  incidence  of  a  surface  wave 

train  from  x  -  -«.  For  a  uniform  wavetraln  the  diffraction  can  be  written 

in  terms  of  the  well  known  diffraction  factor.  We  write  as  the  sum  of  a 

right  going  wave  and  a  left  going  (reflected)  wave  (see  Figure  5.1): 

AC/.*  -ikx.  .  ^  A  fi/.* 

*^11  ”  ^0^0^*^  ®  +<fr  e  )  »  ♦jj  ^0^0^^  ®  +♦  e  ) 

where  in  f^  and  f^,  must  be  associated  with  the  large  (surface  mode) 

solution  of  (5.1.7). 

The  diffraction  factors  and  are: 
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Figure  5.3  b)  The  Locked  Surface  Long-wave  Amplitudes  Iz'l; 

((p'-p)/p  *  O.l),  h’/(h+h')  is  a)  0.4,  b;  0.5  and  c)  0.6. 


133 


{ 


D(-y,  /x) 


(x  >  0) 
(x  <  0) 


(5.3.1) 


^  n(y,  /-x 


(x  >  0) 
(x  <  0) 


where 


D(y,  |xj^/2)  +0  (x  >  0  ,  y  -►  -co) 

D(yt  jxj^/2)  (x  >  0  ,  y  oo) 


(5.3.2) 


and  D  can  be  found  by  a  parabolic  approximation  (cf.  Mel  1983,  pp.  486-490): 


D(y,  /x)  =  {4  +  C(a)  +  i[4  +  S(o)]} 

n 


(5.3.3) 


/irkx 

and  C  and  S  are  Fresnel  cosine  and  sine  integrals. 

The  derivatives  of  D  with  respect  to  x  and  y  are: 

1  3D  .-1 

kD  -57  - 

kD  3y 


(5.3.4) 


(5.3.5) 
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For  a  slowly  varying  wave: 


£  r*'*'  A  j.  a“  A  j.  “ikx, 

“  ^0^'^  V*"!  ■  ♦  V*"!  ^  > 

g  g 


♦il  ■  fflf*''  *0<'l  -  ^  ♦'  *0<'^1  c^>  "'“‘’‘j 

g  g 


(5.3.6) 


If  the  incident  wave  is  given  by  (5.2.5)  then 


AQ(tj)  =  a  cos  Jitj  g/(2ia)fQ(h')) 


as  in  (3.3.13). 

The  far  field  is  defined  by: 

a  >  0(e"h  (5.3.7) 

since  there  the  leading  order  derivative  is  that  of  the  phase  factors 

e  and  e  (5.3.8) 

and  D  can  be  approximated  in  the  far  field  by  its  asymptotic  value  0  or  1 ,  see 
(5.3.2). 

When  we  substitute  the  form  of  ♦u*  ♦'n  for  diffraction  problem  into 
(5.2.1)  -  (5.2.4),  we  find  right  away  that  the  locked  long-wave  in  the  far 
field  is  given  by  the  solution  of  (5.2.11)  multiplied  by  the  square  of  the 
diffraction  factor  from  (5.3.1),  which  is  summarized  in  the  following  table: 

TABLE  5.1 

Quadrant:  I  II  III  IV 

Right  Going:  #+0111 
Left  Going:  *-0100 
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The  near  field  is  defined  as  the  region  where  the  diffraction  factor 
variation  is  significant.  It  is  the  inside  of  the  "parabola" 

a  =  0(e"S  (5.3.9) 

On  the  scale  of  Xj,  x  is  0(e~^)  and  y  is  0(e~^/2)^  go  that  yj=  0.  At  the 
leading  order  4»io»  near  field  slow  potentials,  are  functions  of  ti* 

but  not  of  yj.  Hence  we  have, 

'*'10°  '*'10  “  ’*'10^*1  '}  (5.3.10) 

This  near  field  slow  potential  should  be  matched  to  the  far  field  poten¬ 
tial: 


♦  lO^yi  =  0  ,  Xj,  tj)  »  t|)jQ  (xj,  t^)  =  “  0  »  *1*  tj)  (5.3.11) 

The  normal  gradient  of  the  slow  potential  should  also  be  matched  as  in  the 
previous  chapters: 


♦loy/yi  =  0  .  *1.  tj)  =  ^20  ^y 


(5.3.12) 


*1’  “  ’*'20  ’^l’ 

1  y 

In  the  present  problem  there  are  no  evanescent  modes,  and  hence  4/20 
represents  a  uniform  current.  To  see  that,  we  note  that  the  forcing  term  for 
ij)20z  at  (z  =  0)  as  given  in  (2.2.6)  and  (3.3.26)  can  be  evaluated  as 
follows  (using  a  well  known  vector  identity): 

®  ^h  •  \  '"11  +  *1  =  \  '''11  "  ’''11  \  ’''iJ 


=  la),|;jj  7^  t|»*j  -  imtl.*!  7^  - 

2  2 

=•  [~  k  'i'jj]  “  1“  ^  *  0 


(5.3.13) 
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since  solves  a  single  Helmholz  equation  and  has  no  evanescent  modes 
This  means  that  \l/20  is  divergence-free  and  we  may  write 


*10v  ^^1  “  °  “  '*'20  “  ♦lOv  “  0  »  *i»  (5.3.14) 

where  'l>2Qy  ^'ss  the  same  value  at  both  matching  boundaries 


yji  “  0  >  y  +  +00  and  yj 


y  -»■  -« 


(5.3.15) 


for  any  given  x^,  t^.  This  Is  again  a  simple  special  case  of  (3.3.34)  where 
the  net  mass  flux  through  the  near  field  must  be  zero.  At  the  outer  limits  of 
the  near  field,  the  y  component  of  Stokes'  drift  Is  zero  since  the  waves  are 
propagating  in  the  x  direction.  Hence,  i|>20y  accounts  for  the  total  drift 
current  In  the  y  direction.  equations  identical 

to  (5.3.11)  and  (5.3.14).  The  near  field  Is,  in  effect,  transparent  to  the 
long-wave  and  the  matching  is  carried  out  between  the  far  field  potentials,  on 
both  sides  of  the  shadow  and  reflection  boundaries.  The  near  field  merely 
smooths  the  transition.  The  locked  waves  are  already  known  and  we  now 
determine  the  free  waves. 


5.4  The  Free  Long-wave 

Let  us  write  the  slow  far  field  potentials,  and  as  sums  of 

L  L '  F  F ' 

locked  long-waves:  ^jq,  (^jg;  and  free  long-waves:  <fjQ,  (J^g: 

'*’10  “  '*’10  '*'10  ’  '*’lo  ’  '*'10  '''  '*’10  (5.4.1) 

F  F ' 

The  governing  equations  for  ^jg,  ^jg  are  then  obtained  from  the  homogen¬ 
eous  counterparts  of  (5.2.1)  through  (5.2.4)  with  3/3xi  replaced  by  Vh» 

1 .  e . , 
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'Ltj  '"'Mo  ■  “ 

-F '  -F  2  F ' 

^20tj  ■  ^20tj  “  ° 

(5.4.2) 

— .p  *  F  ^ 

-gC2o  “  'J'lOtj 

p'(g52o 


Eliminating  ^20  ^20 


3t 


^2  2. ,3 

=1 


^,2  2. 
c*  vp 


6^' 

♦10 

F 

♦10 


From  (5.3.11)  and  (5.3.12)  we  get  the  boundary  conditions: 


(5.4.3) 


♦10 


0^ 

o" 


■♦10 


yj  =  0’ 


(5.4.4a) 


b'" 

^10 


yi  = 


0 

o’ 


10 


'1 


0 

o' 


(5.4.4b) 


'lOy, 


0^ 

o' 


-f 


lOy, 


yi 

yi 


0^ 

o' 


(5.4.5a) 


^lOy, 


y,  -0^ 


Yl  -0 


.L' 

hOy 


=  0 


y,  “  0 


(5.4.5b) 
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L  L ' 

and  are  not  functions  of  yj)*  From  Equations  (5.4.4),  (5.4.5)  and 

L  L '  F  F ' 

the  form  of  ^^gt  Immediately  see  that  ^^g,  <^^g  are  antisymmetric  In 

and  In  y^.  It  Is  therefore  sufficient  to  solve  for  the  free  potentials  In, 

say,  the  first  quadrant.  The  boundary  conditions  are; 

along  the  breakwater: 


♦lOxj  “  "^ICxj 


on  (x,  “  0)  (5.4.6a) 

(yj  >  0) 


and  along  the  line  yj  =  0: 


♦ 

♦ 


F' 

10 

F 

10 


-2ln(t 

e 


1_ 

C 

g 


X 


(y^  “  0) 

(xi  >  0) 


(5.4.6b) 


Because  of  the  form  of  the  boundary  condition  (5.4.6b),  we  expect  the 
solutions  of  <^^Q,  i^jg  to  have  a  wavenumber  2{i/Cg  In  the  x^  direction. 
Equation  (5.4.3)  describes  waves  of  two  types.  One  mode  Is  the  barotroplc 

t 

surface  mode  with  horizontal  wavenumber  2ft/Ci  and  the  other,  the  barocllnlc 
Internal  mode  with  horizontal  wave  number  Zn/C^.  For  a  wide  range  of 
parameters,  2n/Cg  lies  between  these  two,  l.e., 

20 
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This  means  that  the  Internal  mode  of  the  free  waves  will  be  propagating  and 
have  real  wavenumbers  In  both  the  Xj  and  y^  directions,  while  the  free  surface 
mode  will  be  evanescent  In  the  y^  direction  (since  the  Xj  wavenumber  Is 
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greater  than  the  total  horizontal  wave  number). 

Away  from  the  shadow  boundary  for  the  free  waves,  which  will  be  defined 
shortly,  we  have 
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where  the  subscripts  bt  and  be  stand  for  the  surface  (barotroplc)  and  the 
Internal  (barocllnlc)  modes,  respectively,  and 
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The  direction  of  the  Internal  mode  forms  an  angle 
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with  the  positive  axis.  There  will  be  a  parabolic  transition  zone  between 
the  area  Into  which  the  free  Internal  waves  are  emitted  and  the  area  near  the 
y^  axis  where  there  are  no  free  waves.  This  transition  can  be  described  In 
terms  of  a  diffraction  factor,  similar  to  D  of  (5.3.3) 
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and  (x|,  y^)  is  a  Cartesian  coordinate-system  with  x'^lncllned  at  ^  with 
respect  to  the  xj  axis: 
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The  diffraction  factor  Dj  multiplied  by  the  periodic  waves  of  Equation  (5.4.7) 

describes  the  'diffracted'  free  long-wave  field  In  the  first  quadrant  (xj,  yj 

>  0).  The  significance  of  the  existence  of  free  waves  Is  greatest  In  this 

region.  The  free  waves  elsewhere  are  simply  given  by  the  free  wave  field 

being  antisymmetric  with  respect  to  both  x^  and  y^. 

In  Figure  5.4  the  transition  zones  for  the  emitted  free  waves  are 

sketched.  We  see  that  free  Internal  waves  are  emitted  Into  the  "shadow  zone" 

behind  the  breakwater,  where  no  short-waves  or  locked  long-waves  are  present. 

To  determine  the  amplitude  of  the  free  Internal  waves,  we  first  write 

down  the  ratio  of  the  amplitudes  of  the  potentials  for  the  two  modes,  the 

barotroplc  (bt)  and  the  barocllnlc  (be).  For  the  surface  (bt)  mode  of  the 

2 

free  long-waves,  we  replace  all  operators  In  the  equations  (5.4.2)  by 
'222  -'  - 

l/C^  3  /9t^.  We  then  eliminate  ^20  integrate  with  respect  to  ti  and 

determine  the  ratio  ♦bt^^bt’ 
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Replacing  by  C^,  we  get  for  the  Internal  wave  (be)  mode: 
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shadow  n  shadow 
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Figure  5.4  The  Diffraction  Transition  Zones  for  the  Free  Long-waves. 
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Note  that  the  signs  of  and  ^,g  for  the  Internal  wave  mode  are  opposite. 


From  (5.4.6b)  and  (5.4.7)  we  have: 
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Equations  (5.4.12)  -  (5.4.14)  ean  be  Immediately  solved  to  give  ^be» 
and  In  partleular; 
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The  baroellnle  (propagating)  eomponent  of  the  Interfaee  dlsplaeement  is  found 
directly  from  (5.2.1); 
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In  Figure  5.3a,  Z,  the  normalized  amplitude  of  the  locked  long  Internal 
waves,  C20»  plotted  versus  k(h  +  h').  The  resonance  at  (kh  -►  0)  is  due  to 
Cg  coinciding  with  C^,  whereas  the  higher  frequency  resonance  Is  due  to 
coinciding  with  Cj^.  In  Figure  5.5,  the  normalized  ^^c  Equation  (5.4.16) 
is  plotted  versus  k(h  +  h').  When  the  depths  of  the  two  layers  are  equal  the 
high  wave  number  resonance  is  absent  in  the  baroellnle  mode.  This  is  because 
the  forcing  is  barotroplc,  E  and  E'  are  both  large,  but  the  barolinlc  mode  (see 
5.4.15)  remains  finite.  When  the  layer  depths  are  different  the  resonance  of 
the  locked  wave  is  reflected  in  the  free  waves.  Near  the  resonance,  however,  0 
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Figure  5.5  a)  The  Free  Internal  Long-wave  Amplitude  Uhrl: 

(p'-p)/p  -  0.1,  h'/(h+h')  is  a)  0.4,  b;  0.5  and  c)  O.b. 
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Figure  5.5  b) 


The  Free  Internal  Long-wave  Amplitude  Ubcl- 

(o'-p)/o  =  0.01,  h'/(h+h')  is  a)  0.4,  b)  0.5  and  c)  0.6 
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of  Equation  (5.4.8)  Is  nearly  zero.  This  means  that  the  free  waves  penetration 
angle  Into  the  shadow  zone  Is  very  small,  and  the  free  waves  are  almost  parallel 
to  the  already  existing  locked  long-waves.  For  waves  with  larger  wave-numbers 
than  the  resonant  wave-number,  the  Internal  waves  are  evanescent.  The  values  at 
which  becomes  zero  represent  a  change  in  phase  of  the  free  internal 
waves.  We  expect  the  interface  displacement  to  grow  with  decreasing  reduced 
gravity  as  seen  when  comparing  Figure  5.5  (a)  and  (b). 

This  concludes  our  treatment  of  slowly  varying  waves  in  a  stratified 
fluid.  The  significance  of  our  findings  is  in  the  propagation  of  energy  into 
the  sheltered  region  in  the  lee  of  the  breakwater.  It  is  of  interest  to  study 
the  possibility  of  resonance  in  a  basin  sheltered  by  a  peninsula  or  a  strait. 
Additional  aspects  include  the  effect  of  continuous  stratification  and  bottom 
topography.  The  related  problem  of  harbour  resonance  in  a  homogenous  fluid 
will  be  commented  on  in  the  following,  conclusion  chapter. 
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6.  COHCLDSION 


6.1  Siunary 

In  this  work  we  have  strived  to  gain  Insight  Into  the  Interplay  between 
wave  scattering  and  diffraction  by  floating  bodies  and  topographical  features 
-  and  the  wave  nonlinearity.  The  waves  were  assumed  to  be  slowly  varying. 
The  water  depth  and  the  short  wave-length  are  comparable.  We  carried  out  a 
multiple  scale  perturbation  expansion  in  both  time  and  the  horizontal  space 
coordinates.  The  governing  equation  in  open  water  was  formulated,  purpose¬ 
fully  allowing  for  waves  propagating  in  a  few  directions  at  the  same  loca¬ 
tion.  Locked  long-waves  were  found  that  are  forced  by,  and  travel  with  the 
short-wave  groups. 

We  found  that  when  the  short-waves  were  modified  by  a  floating  body, 
depth  variations  or  a  breakwater,  so  were  the  locked  long-waves.  The  changes 
in  the  mean  mass  flux  and  free  surface  elevation  associated  with  these 
modifications  are  balanced  by  the  generation  of  free  shallow  water  waves  that 
are  radiated  away  from  the  region  of  short-wave  modification.  Specific 
problems  were  then  solved  by  studying  the  near  field  of  the  wave  modification 
regions  and  asymptotically  matching  to  the  far  field. 

The  lion  share  of  the  thesis  treats  the  drift  motion  of  a  floating  body. 
Specifically,  the  leading  order  slow  drift  forces  on  a  swaying  horizontal 
rectangular  cylinder  were  studied  in  detail.  Various  mooring  stiffnesses  and 
depth-draught  ratios  were  examined.  The  method  demonstrated  can  be  readily 
applied  to  more  general  body  shapes  and  modes  of  motion.  We  have  studied  the 
practically  important  case  in  which  the  slow  drift  amplitude  is  large  and 
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comparable  to  the  short-wave  length.  The  present  work  justifies  the 
results  of  the  more  intuitive  engineering  solutions  that  have  been  in  use. 
It  also  defines  the  limitations  on  its  validity. 

We  have  found  that  there  is  a  first  order  slow  potential  that  corresponds 
to  a  second  order  pressure  and  second  order  slow  surface  elevation.  This 
pressure  is  at  the  leading  order  slow  pressure  and  has  not  been  considered 
before.  The  force  due  to  that  pressure  was  determined.  It  was  found  to  be 
small  for  quite  a  general  class  of  geometries.  The  slow  first  order  poten¬ 
tial  is  coupled  to  the  slow  second  order  potential.  In  particular,  the  ef¬ 
fect  of  the  second  order  potential  on  the  leading  order  slow  motion  and  waves 
was  derived  without  explicitly  solving  for  any  second  order  potential.  We 
have  not  computed  the  third  order  drift  force  due  to  the  second  order  slow 
potential,  nor  the  damping  due  to  the  interaction  of  the  slow  sway  with  the 
short  waves.  However,  the  simplicity  of  the  formulation  facilitates  ex¬ 
tension  of  the  treatment  to  determine  higher  order  corrections.  In  all  the 
cases  studied  both  periodic  and  transient  envelopes  were  Investigated.  The 
sway  response  and  the  long-wave  field  were  determined.  For  several  config¬ 
urations  resonance  was  Identified. 

We  have  also  studied  the  trapping  of  long  waves  on  a  shelf /ridge  by  obli¬ 
quely  incident  slowly  varying  waves.  Here,  too,  emphasis  was  placed  on  the 
new  aspects:  those  of  partial  reflection  of  the  modulated  short-waves  -  and 
the  resulting  local  second  order  slow  potential  which  is  coupled  to  the  long 
waves.  Due  to  the  depth  difference  between  the  two  sides  of  the  shelf  edge 
there  is  a  net  slow  mass  flux  gradient  associated  with  the  short-waves  that 
contributes  to  the  long-wave  generation.  Free  shallow  water  waves  are  emlt- 
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ted  away  from  the  shelf  edges.  These  free  waves  can  be  trapped  on  the  shelf 
and  resonance  may  occur. 

At  the  leading  order  the  resonance  appears  infinite  in  the  frequency  res¬ 
ponse  formulation.  We  have  chosen  to  study  the  excitation  of  an  edge  wave  by 
a  transient  envelope.  After  the  input  short -waves  and  locked  long-waves  have 
passed,  the  trapped  oscillation  (which  is  a  standing  wave  in  the  cross-shelf 
direction  and  propagating  along  the  shelf)  remains  behind  as  a  wake,  which 
can  be  described  as  a  'ringing  effect'.  This  effect  can  be  of  quiet  a  large 
amplitude.  At  longer  times,  higher  order  effects  that  have  not  been  included 
will  diminish  the  wake.  Results  for  the  response  to  normally  incident  waves 
were  also  included.  No  trapping  occurs  and  the  response  is  passive. 
Possible  extensions  of  the  theory  include  application  to  standing  internal 
waves,  which  have  been  observed  on  seamounts,  and  combining  the  effects  of 
shear  currents  (that  can  also  trap  waves)  and  rotation. 

Finally  we  have  considered  the  diffraction  of  slowly  varying  waves  in  a 
two  layer  fluid.  In  addition  to  the  surface  waves,  interfacial  waves  are 
also  present.  When  the  short  waves  are  diffracted  by  a  long  breakwater,  the 
locked  long  surface  and  internal  waves  are  modified.  We  had  asymptotically 
matched  the  far  field  in  the  lee  of  the  breakwater,  through  the  shadow  bound¬ 
ary  to  the  'lit  zone'.  We  found  that  free  baroclinic  long-waves  are  emitted 
from  the  shadow  boundary  into  both  regions.  Long  waves  are  emitted  from  the 
reflection  boundary  as  well. 

The  free  internal  waves  propagate  obliquely  into  the  shadow  zone  at  an 
angle  that  depends  on  the  ratio  of  the  internal  wave  velocity  to  the  larger 
surface-waves  group-velocity.  These  internal  waves  can  resonate  in  a  basin. 
Even  small  internal  waves  can  be  important  to  mixing  processes  and  to  bio- 
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logical  phenomenon.  Space  shuttle  pictures  of  the  straits  of  Gibraltar  show 
internal  waves  propagating  obliquely  into  the  lee  of  the  straits. 

All  of  the  problems  studied  were  designed  to  bring  out  the  physics  of 
short  waves  -  long  waves  Interaction.  the  formal  assumptions  of  the  asym¬ 
ptotic  theory  were:  1)  small  waves  steepness  and  li)  narrow  wave  spectrum. 
For  steeper  waves  and  wider  spectra,  higher  order  corrections  need  to  be  con¬ 
sidered.  The  simplicity  of  the  present  formulation  facilitates  such  exten¬ 
sion.  The  asymptotolc  results  should  be  compared  with  more  accurate  results 
to  assess  their  range  of  validity.  An  additional  difficulty  is  associated 
with  a  directional  spread  of  the  waves.  When  the  range  of  wave  directions  is 
continuous,  the  difference  wave-numbers  form  a  continuum  and  are  not  strictly 
multiple  scales.  It  is  likely  that  this  obstacle  can  be  overcome  by  observ¬ 
ing  that  the  physical  quantity  of  the  gradient  of  the  slow  potential  is  se¬ 
cond  order.  More  work  along  these  lines  may  allow  easier  treatment  of  three 
dimensional  problems.  The  various  features  described  can  be  combined  with 
additional  ingredients  to  study  realistic,  complex  situations.  Finally,  let 
us  sketch  our  approach  to  the  interesting  problem  of  harbour  reasonance. 
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6.2  Epilogue:  Harbour  Resonance 

Harbours  are  known  to  undergo  seiches  at  their  natural  frequencies  which 
are  of  the  order  of  several  minutes.  Several  mechanisms  have  been  suggested 
to  be  associated  with  the  excitation  of  such  oscillations.  These  include 
surf-beats  and  edge  waves.  Bowers  (1977)  has  studied  the  problem  of  waves 
incident  in  a  narrow  channel  terminated  by  a  narrower  basin.  He  did  not  con¬ 
sider  the  local  field  near  the  basin  entrance.  The  present  theory  shows  that 
since  in  uniform  depth  energy  conservation  translates  to  mass  conservation  in 
the  drift  terms,  the  final  result  is  not  affected.  We  choose  to  study  a  more 
realistic  configuration,  where  the  ocean  is  three  dimensional.  Consider,  for 
definiteness  the  configuration  given  in  Figure  6.1.  The  short-wave  groups 
are  incident  bn  the  breakwater  and  diffracted  by  it.  As  in  Chapter  5  there 
is  no  slow  locked  waves  inside  the  harbour,  while  there  is  a  slow  standing 
wave  outside.  Matching  the  two  regions  through  the  shadow  boundary  (that  is 
studied  by  conformal  mapping)  the  free  slow  oscillations  can  be  determined. 
Since  the  eigenf requencies  of  the  harbour  can  fall  in  the  ranre  of  the  mod¬ 
ulation  frequency  of  the  short  waves,  harbour  resonance  will  be  excited.  In 
contrast  to  the  internal  waves  studied  previously,  the  free  waves  outside  the 
harbour  will  be  evanescent,  since  they  are  faster  than  the  locked  long- 
waves.  A  full  solution  to  that  problem  will  be  completed  in  the  future. 


Figure  6.1  Sketch  of  Harbour. 
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APPENDIX  A 


On  The  Scattering  and  Radiation  of  Regular  Waves  by  a  Seal-lnmersed 

Rectangular  Cylinder 

In  this  appendix  we  supply  some  details  on  the  determination  of  the 
reflection  coefficient,  the  hydrodynamic  force  for  the  scattering  of  regular 
waves,  the  radiated  waves,  and  the  added  mass  and  damping  coefficients  for 
harmonic  sway  motion  of  a  semi -immersed  horizontal  rectangular  cylinder. 

Multiplying  (3.2.14)  by  {Fq,...,F(j}  and  integrating  with  respect  to  z 
from  -h  to  -D  we  obtain,  from  (3.2.16)  and  from  the  orthonormality  of  {Ffj} 
the  following  N+1  linear  algebraic  equations. 

“^‘*0  ^00  °  ®00  I  ®n  ’^n  *^n  ®  ®0n 

n=l 

•  •  • 

•  •  •  (A. 1 ) 


where 


A  A 

-2bn  Pr.M  *  8^  J  B  k  cos  k  B  g„ 

0  ON  “ON  0  0  n  n  n  ^Nn 

n=l 


-D 


P  =  f  dz  f  F 
mn  rat 

— h 


(A. 2) 


[^1  -D  +D  [^] 

g  =  /  /  P  (2)  F  (z’)  G  (z,  z')  dz  dz ' 

ran  ^ ,  m  n 

— n  *"11 


(A. 3) 


From  3.2.16  we  get 


[J]  .  «  P  P 

"nm  ran  k  Im  Ira  k 

1=1  £ 

coth  k^B  ^£n^)tm 

[tanh  k  bJ  K 
£=1  £  £ 


(A. 4) 


From  (3.2.15)  multiplied  by  {Fj,...F(^} 
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(A. 5) 


-2a^P.,  -  y  -A  K  sinh  K  B  g, 
^01  n  n  n  ®ln 

n=l 


-2a^P.„  =  y  -A  K  sinh  K  B  g„ 

0  ON  n  n  n  '’Nn 

n=l 

We  need  to  find  Rg  and  which  are  given  by  (3.2.13a), 


i 

Rc  =  r  f  y  “A  K  sinh  K  B  P.  ]  +  1 

0  S  k  I-  n  n  n  On-l 

n=l 


(A. 6) 


“■o  \  t  [“o  •'oo  ^  I ,  -  \  ^  '  ‘’onl  *  ‘  <*•'> 

n=l 

(A.l)  can  be  rearranged  to  be  a  system  of  equations  in  terms  of  {R;^, 
Bg,...B[^}  and  (A. 5)  can  be  written  as  a  system  of  equations  for  {Rg, 
Ap...A()},  then  using  Cramer's  method  Rg  and  R^  are  computed. 

The  hydrodynamic  force  can  be  evaluated  directly.  The  symmetric 
component  yields  no  net  horizontal  force,  so  that  for  an  incident  wave  of 
unit  amplitude: 


F  =  2  ^A  *  /  *11  “  iwpl(l+R^)Mp  Sq  +  ^  a^]  (A. 8) 

-  n-1 

^0 


where 


M  =  f  f  dz  n“0,l. 


and  from  (3.3.13),  (3.2.7) 


a  a  H  I  I  p 

®0  2la>fQ(0)  ® 


(A. 9) 


Alternatively,  F  can  be  found  from  Haskind  relationship  as  will  be  discussed 
later.  In  Figure  A.l  the  values  of 
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are  plotted.  They  agree  with  those  of  Mei  and  Black  (1969  (Figure  6))  and 
Black  (1970  (Figure  9.9))  within  graphic  accuracy.  The  order  of  the 
determinant  was  10  for  and  9  for  Rs.  ^0  terms  are  used  in  the  series  of 

(A. 4).  The  resulting  errors  are  less  than  one  percent. 

For  the  radiation  problem,  the  potential  is  antisymmetric,  (3.3.9)  is 
replaced  by 

(x+B) 

4l  “  ^O'^O  ®  I  ^  ^n  ® 

n*l 


the  boundary  condition  on  the  sides  of  the  swaying  cylinder  is  (3.2.20). 
Instead  of  (3.2.14),  the  following  integral  equation  holds: 


lM_f»  0.  M  f  -D 

-(-F^+  I  =  /  dz'  G^z,z-)  (z*) 

n**l  n  -h 


(A. 12) 


Mutliplying  (A. 12)  by  {Fq,...Fn}  and  integrating  with  respect  to  z  from  -h 
to  -D,  we  get  for  a  unit  sway  velocity: 


^0  ^00  0  On  n  n  n 

n=2 


(A. 13) 


where 


-q,,  =  q„„  Bn  J  gxi  ®  cosh  K  B 
^ON  0  Nn  n  n  n 

n*2 


OP 

q  =  f  M„  +  y  M  P  /k 
n  k  0  On  n  mn  n 

m*l 


B  =  {Bq,...Bi]}  Is  found  by  solving  the  set  of  equations  (A. 13) 


(A. 14) 


-ikb„  =  +  y  g  B  K  cosh  K  B 

0  0  0  00  Nn  n  n  n 

n=l 


which  is  used  to  determine  bg.  From  (3.3.13)  we  get  the  amplitude  of  the 
radiated  wave  in  terms  of  the  sway  velocity: 


hence 


ia>fo(0) 

^11  \ - i —  ’^O  ^1  ® 

V^l  =  -i.>X^^ 


^  fo(0)bo  e 


(A. 15) 


(A. 16) 


(A. 17) 


Haskind  relation  is: 


-  „  2ikB 
F  =•  - -  R 

fQ(0)  a 


(A. 18) 


In  Figure  A. I  the  values  of  jRj.j  are  plotted.  They  also  agree  with  the 
values  of  Black  (1970,  Figure  3.20)  within  graphic  accuracy  as  do  the  values 
of  F  computed  from  (A. 8)  or  from  (A. 18). 

The  hydrodynamic  force  on  the  swaying  cylinder  can  be  written  in  terms  of 
the  real  valued  added  mass  m  and  damping  coefficient  \  defined  by 


-(-u  m  -  IuA)  “  2la)p  J  dz 


(A. 19) 


From  (A. 11)  and  (A. 16)  we  get: 

»  N  (A. 20) 

-i-J  m  -  1.x)  Xj,  =.  21.P  (i  +  j  +  q„Bg  +  J  B  K  co.h  K  B 

n=»l  n“l 

Since  B  is  known,  (“u^m  -  iuX)  is  determined  from  (A. 20).  The  damping 
coefficient  can  be  related  to  Rj.  via  energy  balance  as  follows:  1 


1  cf .  Newman  (1977) 


X  “  2ajpk 


(A. 21) 


where  bQ  Is  given  by  (A.IA). 

The  value  of  X  computed  from  bg  (or  Rj.  In  Figure  A.l)  according  to 
(A. 20)  and  (A. 21)  agree  each  with  the  ocher. 
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APPENDIX  B 


Derivation  of  R,  the  Reflection  Coefficient  for  the  Short  Waves, 
for  a  Swaying  Block  with  a  Narrow  Gap,  or  no  Gap  Underneath 

As  discussed  In  Section  3.3.2,  the  short-wave  potential  Is  not  affected 
by  the  existence  of  a  narrow  (h/L  =  0(e))  gap  underneath  the  rectangular 
horizontal  cylinder.  The  solution  for  a  sliding  block  was  given  In  Agnon  and 
Mel  (1985a)  and  Is  repeated  here. 

(2.3.4)  Is  written  In  the  form: 


-k  X 


f  +  . +v  ,  Ikx  .  r  .  +  f  n 

<^0  ■  “o*  *0®  I,  V  « 

I  =  I 

ni  I  Ik  «  k  X 

f.  +  y  -  b"  f  e  " 

0  n  n 

n=l 


(ao  e 


+  bo  e 


(x  >  0) 
(x  <  0) 


(B.l) 


from  (3.2.4)  get 


“  ^0  "  ^0 


^0  ”  ^0  =  ®0 


(B.2) 


b  =  b  =  B 
n  n  n 


When  (B.l)  Is  substituted  Into  (3.2.6)  and  the  scalar  product  of  the  resulting 
equation  Is  taken  with  fj^  for  each  n,  we  find 


Ao  "  =■  -  (-B„  F„  +  y  B  F  )F„ 

0  0  k„M  0  0  m  m  0 

0  m=l 


and 


-2 


^  F  (-B„  F„  +  y  B  F  ) 
M  n  U  U  m  m 


(B.3a) 


(B.3b) 


m=l 


where 


F  = 
n 


0 

/  f  dz 
•'  n 
-n 


(B.4) 


Define 


(B.5) 


(B.6) 


(B.7) 


(B.8) 


(B.9) 


(B.IO) 


(B.ll) 
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APPENDIX  C 


Derivation  of  the  Ordinary  Differential  Equation  for  the  Slow  Sway  In  Terws 

of  the  Radiation  Stress 

The  derivation  of  the  Ordinary  Differential  Equation  for  the  slow  sway 
was  outlined  In  Section  3.3.1.  When  (3.3.63)  Is  differentiated  with  respect 
to  tj  and  substituted  Into  (3.3.64),  we  get: 

T  + 

KXio-»h2o^  [(-^-^)  -X  1  <C.I) 

g  /gh  11  11  1 


+  pg  |a2  r2  f 


we  now  substitute  for  ,  using  (3.3.65),  to  get: 

t, 


♦lot  ti 

^10  “  *^1  C  “  ^t,  5 

g  111 


(C.2) 


Rearranging,  we  get,  using  (3.1.6),  (3.1.8),  (3.1.9)  and  (3.3.33). 


^°*^1*^1  2ph/gh 


^lOtj  2phcj  ^10  “  2h^j  hi  hi  ^ 


2(i)k  „ 


(C.3) 


✓gh  1  g 


from  (3.3.13)  It  Is  seen  that 
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(C.4) 


From  (3.1.8,  3.1.9)  we  get 


T  9  9  9 


(C.5) 


hence 


9  2  2  2  2  2 

k  -  a  =  k  (1  -  tanh  kh)  =  k  /ch  kh 


(C.6) 


^  _  g  sh  kh  ch  kh  +  kh 

e  “  ch^  kh 


(C.7) 


C  =  -^  =  ^  th  kh 
k  0) 


(C.8) 


-  <T^  =  -^  (2C  -  C) 

gh  g 


(C.9) 


♦fot,  -  if  <=8  fo<«>  ho  I 


»k  fg(0)  2  C  _  2(gh  -  'i),  I.  |2  ,2 


gh  - 


■]  hor 


(C.IO) 


Defining 


(ok  f.(0)  C  „  ,  .  , 

g 


s  =^-i 

^  -  C  2 


(C.ll) 
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we  may  write  (C.3)  (using  (3.3.13))  as: 


(3.3.66) 
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APPENDIX  D 


Evaluation  of  the  Steady  Drift  Force  on  a  Floating  Rectangular  Cylinder 

In  Regular  Waves 

In  this  appendix,  we  evaluate  the  steady  drift  force  on  a  sliding  block 
in  regular  waves  via  direct  integration  of  the  mean  pressure  on  the  body,  and 
demonstrate  that  it  is  equal  to  that  given  by  the  general  formula  derived  by 
Maruo  (1960)  and  by  Newman  (1967)  via  momentum  flux  consideration,  and  used  in 
(3.3.19).  We  shall  then  generalize  the  result  to  a  floating  cylinder. 

In  the  right  hand  side  of  (3.3.8)  we  let  P20  be  given  by  (3.3.9)  with 
’I'lOtj  bo  be  zero,  for  regular  waves.  We  get  for  the  drift  force: 

* 

F  =  -p[  /  P20  dz  -  /  p^Q  dz  +  6(ia))|»jj  +  *) 


+  6(g  qp] 

P20  "  h'*'ll|  ^11  +  *)  +  gz 


(D.l) 


(D.2) 


The  contributions  of  the  terms 


2  ,  * 
j  and  (ia)i{)^j^  +  *) 


(D.3) 


to  the  integral  are  zero  because  of  the  kinematic  boundary  condition  (3.2.4). 
We  shall,  however,  retain  the  first  of  these  terms  for  completeness,  and 
because  it  facilitates  the  generalization  to  a  floating  cylinder.  The 


hydrostatic  terra  can  clearly  be  discarded  of  since 

f  z  dz  -  f  z  dz  =*  0 

Sq  ^0 


(D.4) 


Integrating  ipHz  ^  by  parts,  we  find 
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/  ♦  ^  dz  -  -  /  +  1|»11  i|»Ji  I 

.±  XX  oi  “  2 


-0 

h 


*  ^2 

/  i|»n  'I'll  i"  hill  z-0 

.±  XX 


‘'0 

where  use  has  been  made  of  (2.3»1),  (2.3«2)  and  (2»3.3)( 
From  3.2.4  we  get 


/  X  *)  dz  -  /  (lu4(  X  *)  dz  “  0 

<;“  *  s'*"  * 

Substituting  (D.5)  and  (D.6)  Into  (D.l)  and  using 

lot  . 

5ii  g  hi 

we  get 

?“  /  ^’•'11  ’•'11  hii  "  / ^’*'11  *11  1*11 

P  XX  '  x'  _+  XX  X 


(D.5) 


(D.6) 


(D.7) 


(D.8) 


If  we  write  In  the  form: 

*11  ”  *11  ■"  *11  ^0^®  ‘^A  ®  ^^0  k  n  ®  " 


n*! 
k  (x+B) 


(D.9) 


f  (x  <  -B) 
n 


where  for  (x  >  B)  i|»jj  Is  written  in  the  same  form,  with  i|»jj  symmetric  with 

A 

respect  to  (X*0)  and  antisymmetric. 

S  A 

1|»JJ  Is  decomposed  Into  and  \|ijj  so  that  the  result  can  be  generalized 
to  the  floating  rectangular  cylinder  case. 

Making  use  of 
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(x  =  B)  “  =•  -B) 


(D.IO) 


=  B)  =  =  -B) 


(D.ll) 


we  get  from  (0.8)  that: 


*  *  ^k 

^  “  2  /  i|)^j  +  *)]dz 

^  XX  XX  XX 


(0.12) 


substituting  (0.3)  into  (0.12)  and  making  use  of  the  orthoromallty  of  {fg. 


f , . . .}  we  get: 


^  =  2  (-k^  A  B*  (1+R  )(1-R*)  -  ^  k^  A  B* 

0  uu  A  S  nnn 

^  n“l 


-  k^  A_B*  (1-R.)(1-R*)  -  I  k^  A^  B*]  +  * 
00  A  S  nnn 

n=l 


i  -  4  AqB*  (1  4  R*)  4  * 


(0.13) 


(0.14) 


For  a  wave  incident  from  (x  =  -•)  with  amplitude  A,  we  get,  from  (3.3.13)  and 
(0.9) 


A  =  B  »  _ _ 

*0  “O  4lo)fQ(0) 


(0.15) 


The  reflection  and  transmission  coefficients  are  given  by 


(0.16) 


Substituting  (0.16)  into  (0.14)  we  get: 


-  -  8  k  (1  +  R  -  T  )  Re  A_B 


Substituting  (0.15)  into  (0.17)  we  get 
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(D.18) 


but 


2C 


g/(ii  [sh  kh  ch  kh  +  khl/ch  kh 
C  g/ti)  th  kh 

making  use  of  (3.3.34)  and  (D.19),  (0.18)  can  be  written  as: 

|2 


(D.19) 


F  -  pg  A  R 


2  _£ 


(D.20) 


This  result,  (3.3.12),  can  be  generalized  to  a  floating  rectangular  cylinder. 
If  D  Is  smaller  than  h,  we  need  to  subtract  from  the  right  hand  side  of  (0.8) 
the  following  expression 


-0 


/  'I'll  Uii  I  )  dz 

-h  XX  x' 


X  =•  -B 


X  “  B 


(0.21) 


Integrating  the  first  term  by  parts,  we  get,  as  In  (0.5)  that 

-0  *  -0  , 

/  dz  -  / 

-h  XX  -h  2 


dz 


Since 


'11 


0  (x  =  ±B;  z  -  -h,  -0) 


(D.22) 

(D.23) 


(cf.  (2.3.3),  (3.2.5));  Thus  (0.10)  may  be  written  as: 

x  -  -B 

-h  '  **'  X  =  +B 


dz 


(D.24) 


which  Is  the  net  mean  horizontal  momentum  flux  Into  the  volume  under  the 
floating  cylinder.  This  quantity  must  be  zero  since  there  Is  no  momentum 
transfer  from  the  swaying  cylinder  to  the  fluid  below  It.  We  conclude  that 
(3.3.12)  Is  valid  for  swaying  floating  rectangular  cylinders. 
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APPENDIX  E 


Derivation  of  the  Added  Mass  and  Damping  Coefficient  for  Large  Aiq>lltude 

Slow  Sway 

The  presentation  in  this  appendix  is  based  on  the  work  of  Professor  H. 

S.  Choi  from  Korea  who  was  visiting  Parsons  Laboratory  during  1985, 

The  method  of  a  solution  involves  asymptotically  matching  a  near  field 
solution  and  a  far  field  solution.  The  near  field  is  within  a  few  short 
wavelengths  from  the  body.  In  the  near  field  a  rigid  lid  upper  surface  can 
be  assumed  for  a  slow  potential.  We  now  apply  the  same  reasoning  which  led 
to  Equation  (3.3.59).  Let  the  ambient  flow  velocity  be  ev,  where  we  assume 
that  the  sway  motion  has  a  unit  amplitude.  The  sway  velocity  is  thus 

-2int. 

Re  -2ine  (E.l) 

and  its  acceleration  is 

2  -21nt, 

Re  -4fi  e  (E.2) 

Let  the  near  field  potential  be 

-2lnt. 

t)/  =  Re  i|;'  e  (E.3) 

\p'  will  have  the  following  outer  approximation  in  the  body  coordinates  (l.e. 
x”  =  X  -  Re  exp(-2lntj)); 

i|»'  -*■  e(v'  +  2in)  (x”  ±  c)  (kx"  -►  ±  «<>)  (E.4) 

where 

-2l£2t, 

V  =  Re  v'  e  (E.5) 

In  the  fixed  frame  of  reference,  we  have  instead 

(j)'  e[v'x  ±  (v*  +  2in)c]  (kx  ±  a>)  (E.6) 

c,  the  blockage  coefficient  was  given  by  Flagg  and  Newman  (1971).  In 


normalized  form,  we  have: 
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C  B.2_  «/Hv2 


(E.7) 


In  the  far  field,  the  radiated  long  waves  are 


^  “  Re  4)*  e 


-2int, 


(E.8) 


where 


,  ‘“i  *1 

=  ±  a  e 


(xj  >  0) 
(xj  <  0) 


(E.9) 


2  2 

a  Is  an  amplitude  function  and  ghk^  =  (2n)  ,  Use  has  been  made  of  the  fact 
that  the  slow  velocity  field  due  to  sway  Is  even  In  hence  Is  odd  In 
x^.  Note  that  ^  Is  a  solution  to  the  long-wave  equation;  we  do  not  Impose  a 
rigid  lid  condition  In  the  far  field* 

The  Inner  expansion  of  Is: 


±  a  [1  +  Ik^  Xji  +  0(kxp  ] 


2,  (x,  >  0) 


(xj  <  0) 


(E.IO) 


Matching  the  asymptotic  expansions  (E.6)  and  (E.IO)  Immediately  yields: 


a  “  v' 


(E.ll) 


a  *  e(v'  +  2ln)c 


(E.12) 


from  which  we  find; 


~e2lnclkj 

elkjC-1 


(E.13) 


The  radiation  coefficient,  which  Is  defined  as  the  amplitude  of  the  potential 


4iq  per  unit  amplitude  of  sway  displacement.  Is  given  by 

„  .  -e21nc 

elkjC-1 


(E.14) 
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To  find  the  added  mass  y,  and  damping  coefficient  A,  we  compute  the 
horizontal  dynamic  force  on  the  body: 

-2lfit 

F  “  pi  /  '1'^  dz  -  J  dz]  =  p  Re  -2lne  [  J  4»'  dz  -  J  ip'  dz]  = 

s“  ^  ^  s“ 

-2lftt. 

=  p  Re[-2lite  (ha  -  2lfi  2BD)]  (E.15) 

(cf.  Beck  &  Tuck  (1972)). 

The  added  mass  and  damping  coefficient  are  defined  by: 

,  -2lJit, 

-F  =  V  '  2ifiA)e  (E.16) 

From  (E.15)  and  (E.16)  we  find: 

-4fi^(p  +  M)  -  2lftA  =  4ft^ph  (E.17) 

since  M  =  2pBD. 

When  the  gap  H  Is  relatively  large,  the  blockage  coefficient  tends  2BD/h 
(this  Is  because  the  Increase  of  the  current  In  the  gap  Is  approximately: 
h/(h-D)  -  1  -»■  D/h  as  D  0). 

The  normalized  potential  jump,  2a/(-21fi),  is  approximated  by  -2c.  This 
Implies  that  the  force  Is  mainly  In  phase  with  the  acceleration: 

p  -  p(-2BD  +  2hc)  -►  2pBD  (E.18) 

which  are  the  values  for  slow  sway  in  the  absence  of  free  surface  effects. 
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